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Mass Spectra of the Deuteromethanes 
By Vernon H. Dibeler and Fred L. Mohler 


The four deuteromethanes have been synthesized and their mass spectra obtained with 


a Consolidated mass spectrometer. 
are equal within the estimated uncertainty. 
and CH,D, but is significantly low for CD,. 


not predictable from CH, by statistical considerations alone. 


The sensitivities of the molecule ions of CH, and of CD, 
The total ionization is the same for CH,, CH;D, 
The mass spectra of the deuteromethanes are 


In agreement with previously 


published data on CH,;D, probability factors of 1.18 and 0.45 were observed for the dissocia- 
tion of an H atom or a D atom compared with the dissociation of an H atom from CH,. 


Whereas these factors remain 


same molecule, they differ markedly for the three protiodeuteromethanes. 


approximately constant for subsequent dissociations in the 


Wide divergence 


is observed for some probability factors obtained in dissociations of two or more atoms. 


I. Introduction 


The use of deuterium and deuterated molecules in 
tracer experiments has given much information on 
reaction mechanisms, kinetics, and catalysis. One 
of the most rewarding methods of observing isotope 
exchanges and analyzing the reaction products of 
tracer experiments is that of mass spectroscopy. In 
all but the simplest experiments, however, this 
method requires calibration patterns of the substances 
encountered in the investigation, and it is obvious 
that eventually a complete catalog of mass spectra of 
deuterated and other isotopic molecules will be re- 
quired. In addition to their practical value, mass 
spectra of isotopic molecules are of considerable 
theorectical interest. Studies of a few simple deuter- 
ated hydrocarbons have been reported by Delfosse 
and Hipple [1], ' Evans, Bauer, and Beach [2], Mohler 
and Dibeler [3], and Turkevich, Friedman, Solomon, 
and Wrightson [4]. These studies indicate that al- 
though the sensitivity (ion current per unit of pres- 
sure) of the molecule ion and the total ion current for 
all dissociation the same for isotopic 
molecules, the effect of deuterium substitution on the 
relative probabilities of the various dissociation proc- 
esses Is not predictable from statistical considerations 
alone. Appropriate corrections must be made for 
differences in dissociation probability of C—H and 
(—D bonds, as well as the difference in behavior of 
H atoms in normal hydrocarbons and H atoms in 
deuterohydrocarbons. 

\lthough information has been published on the 
Raman [5] and infrared [6] spectra of all of the 
deuteromethanes, mass spectra of only CH,;D |[2, 4] 
and CD, [4], and a partial spectrum of C HD. [4] have 
been reported. Accordingly, we report the complete 
inass spectra of all of the ‘deuteromethanes to extend 
the published data on this group of molecules. Also, 
we wish to reexamine the statistical relationship of 
Cll, to the deuteromethanes for the complete series 
now available. 


processes 1s 


II. Experimental Details 


‘he mass spectra were obtained with a Consoli- 
dated mass spectrometer by using conventional 
techniques. The energy of ‘the ionizing electrons 


s in brackets indicate the literature references at the end of this paper 


was nominally 50 v. The temperature of the ioni- 
zation chamber was about 245° C. Possible traces 
of hydrogen were pumped off of all samples at liquid 
hvdrogen temperature, and the sample storage bulbs 
were cooled to liquid nitrogen temperature before 
samples were admitted to the mass spectrometer. 
The pressure of the sample in the inlet reservoir was 
measured with a diaphragm-t ype micromanometer [7]. 

Methane was “‘research grade”’ material obtained 
from Phillips Petroleum Co. and used without 
further purification. The purity, as stated by the 
manufacturer, was 99.7 mole percent. The most 
probable impurities were nitrogen, ethane, and 
carbon dioxide. 

Mono-, di-, and trideuteromethane were prepared 
by the reaction of lithium aluminum deuteride on 


the appropriate bromo-methane, according to the 
method of Finholt, Bond, and Schlesinger [8]. The 
apparatus and techniques were similar to those 


previously used in the reduction of halomethanes 
with lithium aluminum hydride [9]. 

In addition to the methanes, small amounts of 
heavier hydrocarbons and hydrogen were observed 
in the products. These were effectively removed as 
described above before the sample was admitted to 
the mass spectrometer. As the attainment of high 
isotopic purity became more difficult with increasing 
numbers of D atoms in the molecule, no attempt 
was made to prepare tetradeuteromethane by this 


method. The CD, was very kindly given by Leo A. 
Wall. It was prepared by the reduction of carbon 
dioxide with deuterium over a supported nickel 
catalyst. The product contained 1.8 mole percent of 
CHD,;. Hydrocarbons heavier than CD, were not 


detected by the mass spectrometer. 

The observed sensitivity of CH, (31.2 div/micron) 
and of CD, (31.9 div/micron) were considered equal 
within the estimated experimental uncertainty. The 
average value (31.5) was assumed for the sensitivities 
of CH,D, CH,.D., and CHD,, and on this basis the 
mass spectra of the deuteromethanes were corrected 
for isotopic impurities. After correcting the mass 
spectrum of “CH,D” for the contribution of ions 
containing carbon-13, a residual peak on mass 18 
above the normal H,O background indicated the 
presence of 1 percent of CH, SDs. The sum of the 
partial pressures of CH,;D and CH,D, (obtained from 
the mass 17 and 18 pe ak heights and the sensitivities) 
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subtracted from the total sample pressure indicated 
1 percent of CH, also present. An exact correction 
of the mass spectrum was made for the latter, and 
an approximate correction was made for the CH,D,. 

Similarly, the mass spectrum of the sample of di- 
deuteromethane showed a trace of CHD;. The total 
partial pressures of the CHD, and CH,D, subtracted 
from the sample pressure indicated the presence of 10 
percent of CH,D plus CH,. Approximate corrections 
were made to the mass spectrum assuming the amount 
of CH, to be negligible. 

The isotopic purity of the trideuteromethane was 
the least satisfactory of the group. The product ob- 
tained from the reaction of lithium aluminum deu- 
teride with bromoform was only 60 percent of CHDs. 
The partial pressures of the CD, impurity (obvious 
from the mass spectrum) and the CHD, subtracted 
from the total pressure gave a partial pressure as- 
sumed to be CH,D,. The corrected mass spectrum 
was in good agreement with the spectrum derived 
from a second sample of CD, obtained from Wall, 
containing 30 percent of CHD,. This sample con- 
tained only a small amount of CH,Dz, as indicated 
by the partial pressure computations. Finally, the 
mass spectra summarized in the following section were 
obtained by successive approximations, as indicated 
above. With the possible exception of CHDs,, this 
process should give satisfactory results. Fortunately, 
CHD, was not an appreciable impurity in any of the 
preparations. The estimated uncertainty in the 
principal ions of the spectra of CH,, CH, D, and CD, 
is about 1 percent. It is probably some ‘what higher 
for CH,D, and is estimated as 5 percent for CHD,. 


III. Results 


The relative ion intensities obtained for methane, 
and the four deuteromethanes are given in Table 1 
Correction has been made for ions containing car- 
bon-13. The sum of the carbon-group ions for each 
molecule is given near the bottom of the appropriate 
column. The measured sensitivity (molecule-ion cur- 
rent per unit of pressure) is also given for CH, and 
CD,. No doubly charged ions were observed. One 
ion due to a metastable transition was observed in the 
spectrum of CD,. This ion observed at m/e= 16.2 was 
only 0.05 percent of the abundance of the molecule 
ion and was attributed to the delayed dissociation 
CD,*-—-CD;*+D._ No independent check on the iso- 
topic composition of the methanes was employed. 
However, the mass spectrum of CH,D is in excellent 
agreement with previously published data from two 
different laboratories [2, 4]. Also, the mass spectrum 
of CD, and the odd-mass ions of CHD, are in good 
agreement with the spectrum of CD, and the partial 
spectrum of CHD, previously reported by Turkevich 
et al. [4]. 

In the dissociation patterns of CH;D and CH,D,, 
some of the mass peaks come from two different 
ions; for example, mass 14 is CH,* and CD*. In 
order to obtain information on the relative probabili- 
ties of dissociating H atoms and D atoms from 


various isotopic configurations, some assumption 


TABLE 1. Mass spectra of methane and the deuteron 
for 50-v electrons 
m/e CH,  CH;D |} CHyD, CHD, CI 
1 7.05 5.92 4.40 2.05 
2 0.34 0. 68 1.24 1. 86 aT 
; il 0.18 0.11 
4 . 04 09 0.2 
12 2. 57 2 46 2.39 2. 30 2.19 
13 8. 21 4.90 2.80 1. 42 
4 16.3 &. 80 6. 39 6. 41 7.2 
1h 86.1 20.9 9.79 6.75 
16 100.0 77.2 30.7 13.2 12.5 
17 100. 0 62.4 51.1 
18 100.0 27.7 &3. 1 
19 100.0 
a» 100 
otal C-group ions 213. 2 214.3 214 208.9 N44 
Sensitivity 1.2 4] 


must be made to derive the relative amounts of each 
contributing ion obtained in the dissociation proe- 
In the case of CHDs,, this problem does no 
arise, as each peak comes from a single type of ion 
[t is unfortunate that experimental errors are greatest 
in this spectrum. In the computations for C H;D 
and CH,D,, it will be assumed that the probability 
of removing one atom (either H or D) 1s equal to 
the probability of removing one H from CH,; the 
probability of removing two atoms is equal to that 
of removing two H atoms from CH,, and so forth 
This assumption is consistent with the fact that th 
sums of the ions in CH,, CH;D, and CH,D, an 
nearly equal, as shown in the following computations 
Starting from the low-mass end of the CH;D speec- 
trum, the contribution of each type of ion is com- 
puted in the following manner: Mass 13 results from 
CH* only and is equal to 4.90 (table 1). As CH 
and CD* are both formed by the dissociation of 
three atoms, the above assumption permits us to 
calculate the CD* abundance as the difference: CH 
in CH, (8.21) minus CH* in CH;D (4.90) equals 
3.31. Further, the total abundance of ions of mass 
14 in CH;D is 8.80; consequently, the CH,* abun- 
dance is given by 8.80—3.31=5.49. Continuing ir 
this manner, one computes the abundance of CH,D 
as 76.0 instead of the observed value 77.2. Theres 
a residual error, because the sum of masses 15, | 
and 15 in CH, is not exactly equal to the sum ol 
masses 13 to 16 in CH;D. If the values for CH, 
are multiplied by 1.01, the sums are exactly equal 
and no residual error occurs. This device mere 
distributes any error in the experimental data or 0 
the assumptions. Similarly, for CH,D, a factor o! 
1.013 eliminates the residual errer. In this manne! 
the ion abundances listed in table 2 were derived 
There are included in table 2, the @ priori probabil 
ties of removing H or D atoms from the three mol 
cules to give the various ions. The abundances 
the various ions are not in the ratios of the a pr 
probabilities. In CH;D, for example, the obser: 
abundance of the CH.D* ion is not 4X87.0=65~ 
but 1.183 times 65.2 (77.2), where 87.0 is the eo 
rected abundance of CH} in CH,. On the oth 
hand, the CH# ion in CH;D is 0.45 4X 87.0 = 
The factors 1.18 and 0.45 are a measure of the pre! 


esses. 
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Taste 2. Computed ion abundance in isotopic methanes 
CHsD CHD, CHD; 

\ pri- \ pri- A pri- 

Abun- ori A bun- ori A bun- ori 

Ton dance | proba- lon dance proba- Jon dance | proba- 
bility bility bility 
H:D | 77.2 84 CHD:| 62.4 ls CDs 27.7 \% 
CH 9.8 4 CH,D 24.9 le CHD, 51.1 A 
HD 11.1 lo CDs 5. 84 ly CDy 13.2 ly 
Hy 5.40 le CHD 9.79 24 CHD 6.75 ly 
CD 3.40 4 CH, 0.87 le cD 6. 41 ry 
CH 4.90 % cD 5. 52 ly CH 1.42 yy 

CH 2. 80 le 


ability of dissociating an H atom or a D atom from 
CH,;D compared with the dissociation of an H atom 
from CH, These probability factors are denoted 
by the symbols H (H;D) and D (H;D) in table 3. 
The letter outside the parenthesis indicates the atom 
dissociated, while the letters inside the parenthesis 
indicate the configuration before the dissociation, 
the presence of a carbon atom being understood. 
One can lasinty paveninoreed H (H.D,) and D (H,D,). 
In CHDs;, one uses the observed sum of CD? and 
CHD} without reference to CH, to compute the 
corresponding factors H (HD;) and D (HD;). 
Removal of two H atoms from CH,D has an a 
priori probability of one-half compared with CH,, and 
the probability factor can be expressed as the product 
of the probabilities H (H,;D) and H (H,D). There is 
an assumption involved in expressing the probability 
of removing several atoms as the product of the in- 
dependent probabilities, and it is of interest to see if 
the assumption leads to consistent results. There 
are two ways of removing H+-D atoms from these 
molecules; removing H first or removing D first. 
For CH,D, the probability for the observed CH; 
compared with the @ priori value can be expressed as: 


CH; from CH,D 


CH; from CH, 4 H(H,D)x 


D(H,D)+ 1/4 D(H,D). 


All terms are known except D(H,D), which can be 
computed from the above. There are three ways of 
removing 2H+D atoms to give CH*, and all con- 


stants are known except D(HD), which can be 
derived in a manner similar to the above. 


lante 3. Probabilities of dissociating H or D atoms from 
deuteromethane tons 
CH sD CHD, CHDs 
H(HyD) =1.18 H(H 2D») =1.43 H( HDs) =1.41 
D(HyD) =0.45 D( HD») =0.57 D( HDs) =0.87 
H(H,D)=1.14 H( HD») =1.49 D( HD») = .78 
H(HD) =1.21 (DCH2D) =0.54) (DCHD) =1.07) 
D(H 2D) =0.73 * H(H,D)=1.16 H(HD»)) =1.48 
D(HD) = .78 +b D(H Dy») =0.79 (HCH D) =0.95 


(HCH D) =—1.05) 
(DCH D) =0.90 


rived by using value of D(H D,) from CHD, spectrum 
ved by using H(H,D) from CHsD spectrum 





Computations for the dissociation of CHD, can 
be carried out in the same manner by using only the 
data from the CHD, spectrum. In CH,D, more 
constants are involved, and data from CHD, or 
CH,D must be used to derive one of the constants. 


IV. Discussion 


The constants derived from the three spectra are 
not entirely consistent. It should be noted that 
errors are necessarily cumulative in computing 
probabilities of removing successive atoms from the 
methanes, consequently computational errors are 
greatly magnified. H(H,D) can be derived either 
from CH,D or CH,D,, but in the latter case the 
value of D(HD,) from CHD, must be used. The 
pairs of values derived for H(H,D) or for D(HD,) 
are entirely consistent in the three spectra. Values 
for H(HD,) from CH,D, and CHD, are equal within 
experimental error. Other values do not check well. 
The values for D(H,D) are 0.73 from CH,D and 
0.54 from CH,D,. The former is probably more 
reliable, as the CH; ion in CH,D, (table 2) is very 
small and any uncertainty in computing table 2 
makes more uncertain the value derived from CH,D,. 
Values for H(HD) and D(HD) can be derived from 
all three spectra, and the values are quite divergent. 
Computations involve products of three probabili- 
ties, and thus errors are magnified. The values de- 
rived from the CHD, spectrum, D(HD)=1.07 and 
H(HD)=.95, are definitely contrary to expecta- 
tions, as it seems to be quite generally true that the 
probability of removing H from a radical containing 
D is greater than 1, and the probability of removing 
D from a radical containing H is less than 1. In 
table 3 values that are considered less reliable are 
inclosed in parentheses. 

It is concluded that the assumptions used to derive 
the probability constants are at least approximately 
true. The inconsistencies shown in table 3 are no 
greater than are to be expected if there are experi- 
mental uncertainties amounting to several percent of 
the maximum peak, and this is very likely true in 
the CH,D, and CHD, spectra. 

In criticism it can be said that it is probably not 
accurately true that the probability of removing one, 
two, or three atoms from these molecules is equal to 
the probability in CH,. If one uses CD, as the basis 
of comparison, somew hat different values are derived. 
There 1s also reason to doubt whether the probability 
of removing two or three atoms is accurately equal to 
the product of the probabilities of removing single 
atoms, for appearance potential measurements show 
that H atoms can be removed in pairs to give H, and 
not 2H. However, L. G. Smith has shown [10] that 
most ions in methane are formed by the dissociation 
of H atoms. In any event, the assumption seems to 
be approximately true, and without this assumption 
there is little possibility of extending data on prob- 
ability constants observed in one molecule to 
another. It has been pointed out by E. Wells 
that there are enough interrelations between con- 
stants to derive all the constants without the 
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assumptions that probabilities of removing atoms are 


Evans, Bauer, and Beach lead to nearly the 





equal to those found in methane. Unfortunately, results in CH;D, their assumptions did not sev to 
one must base such a derivation on the values derived be generally applicable to the three spectra (|/,D 
from CHD, that are considered least reliable experi- CH,D,, and CHDs, necessitating the dit 
mentally. approach used in this work. 
V. Conclusions 
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First Dissociation Constant of Succinic Acid from O” to 
50° C and Related Thermodynamic Quantities 


By Gladys D. Pinching and Roger G. Bates 


of succinic acid was determined at intervals of 5 





The first dissociation constant, A,, 
degrees from 0° to 50° C from measurements of the electromotive force between hydrogen ‘1 
and silver—silver-chloride electrodes in cells without liquid junction. In order to reduce f 
the correction for “overlapping” of the two acid groups, which are of nearly equal strengths, 
the 10 solutions studied contained twice as much succinic acid as acid succinate ion The 
residual correction of about 0.02 in log A, was made with the aid of values of the second ~ 
dissociation constant from an earlier determination n 
The first dissociation constant between 0° and 50° can be expressed by a 
log K,= 1206.25/ T— 3.3266 + 1.011697 7, ' 
where 7 is the absolute temperature. The changes of free energy, heat content, entropy, 
and heat capacity for the first dissociation step of succinic acid in the standard state were S 
calculated from the dissociation constant and its temperature coefficient g 


I. Introduction 


A determination of the second dissociation con- 
stant of succinic acid from 0° to 50° was recently 
reported by the authors [|1].'. The two acid groups 
of this dicarboxylic acid are of nearly equal strengths, 
and there is considerable overlapping, or simultaneous 
dissociation. Hence, it was necessary to estimate 
the first dissociation constant, A,, in order to com- 
pute the ionic concentrations required for an accu- 
rate determination of the second dissociation con- 
stant, A». Values of A, were therefore calculated 
from electromotive-force measurements of solutions 
of sodium acid succinate with added sodium chloride 


Figures in brackets indicate the literature references at the end of this paper 


from 0° to 50 
product of the two constants {2}. 


Although adequate for correcting the deter! 


tion of A, these values of log AK. were estimat: 


have an uncertainty of +0.01. Hence, it is do 
if the earlier results for A, are of sufficient acc 


to warrant 


the calculation of the heat-conten! 


by means of a determination of U 


entropy changes that are of considerable theoreti 


interest. It seemed desirable to undertal 
accurate redetermination of the first 
constant in order to obtain reliable informatio 
cerning the thermodynamic quantities asso 
with the first dissociation step. 

Ten mixtures of hydrochloric acid and s 
acid succinate were prepared, each contain 


disso¢ 
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Paste 1. Electromotive force of cell 1 from 0° to 50 
Pemperature, ° C 
m 
0 5 10 15 20 2 ‘) ) 45 “) 
Li t t r t 7 r t 
0 OvesY 0. 50441 0. SOST2 0. 50698 0. SORIS 0. Sog29 0. 51085 Oo. 51145 0. 51250 0. 51351 0. 51448 0. 515 
09062 SOONG 50812 M40 5LOST S172 51290 51399 51503 51620 ‘ATS 
OS029 MO84 51107 51242 51375 51500 51622 51740 51853 51958 52073 
O7054 51200 51414 515 4 51697 5ISSI 51954 52078 52182 52302 52408 
050385 42112 52272 52424 52574 52725 52861 52000 53133 53262 AS392 
04983 2735 S28S8 33022 M41 4 
02977 4086) 4248 “M410 4575 ‘Mr24 
O1QS7 131 M20 549s MiOOT MAT 
O15114 es wid 238 ‘ONO iG6Y 
OOGS00) iOS 7174 7383 Tuyl SSlvO 





atter at one and one-half times the molality, m, of 
the former, and the electromotive force of the cell, 


Pt; Ho(g, 1 atm), NaH succinate (1.5 m), HCI (m), 
AgCl; Ag, (I) 


was measured at intervals of 5 degrees from 0° to 50 
The first dissociation constant of succinic acid was 
caleulated from the experimental data. Corrections 
were applied for the amount of bivalent succinate 
ion formed in the solutions by dissociation of the 
acid succinate ion. The changes of free energy, 
heat content, entropy, and heat capacity for the 
dissociation of | mole of succinic acid in the standard 
state were evaluated from the dissociation constant 
and its temperature coefficient 


II. Experimental Procedures and Results 


The stock solution of hydrochloric acid was pre- 
pared by dilution of a portion of redistilled acid. 
The mean deviation of five gravimetric analyses from 
the mean molality was less than 0.02 percent. 

Sodium acid succinate, prepared from reagent- 
grade succinic acid and a standard solution of sodium 
hydroxide, was recrystallized from water with the 
addition of sufficient sodium hydroxide to bring the 
salt to the proper composition. After samples of 
the salt were dried at 110°, they were found to neu- 
tralize 100.03 percent of the theoretical amount of 
standard alkali. The titration was carried to pH 
8.7 with the use of phenolphthalein and a color 
standard. 

Nitrogen was passed through each buffer solution 
for 2 hr to remove dissolved oxygen. The cells 
were filled and placed in the water thermostat where 
they were kept at 25° for about 6 hr to insure initial 
quillbrium. The emf measurements from near 0 

5° were made on the second day, and from 25‘ 
to 50° on the third, after which the temperature of 

thermostat was lowered to 25° for a final measure- 
ment The average difference between the initial 
and final emf of the 10 cells was less than 0.08 mv, 
or 0.0013 in log Ay. The electrodes were prepared 
i (he same manner as for the determination of the 
second dissociation constant [1]. 

wh emf value, &, given in table 1 is the mean of 

tentials between two pairs of electrodes in the 





same cell. The moialityv, m, of hydrochloric acid 


appears in the first column 


III. Calculation of the Dissociation Constant 


Withdrawal of electrical energy from cell | causes 
the following electrochemical reaction to take place: 
 H,(g, 1 atm) + AgCl(s) = Ag(s) + H*+ Cl (1) 
The initial, reversible potential of the cell depends 
upon the extent to which the system is displaced 
from the equilibrium condition, that is, upon the 
difference between the reaction quotient, ©, and 
equilibrium constant, ©, for reaction 1. If the 
usual conventions regarding the standard state with 
respect to which activity, @, is measured are observed, 
we may write 
() ay,a.\, (2) 
and 
AF=RT In 0/0,—AF RT In Q. 3) 
in eq 3, AF is the molal free-energy change in reac- 
tion | for any particular state of the system and 
AF® is the corresponding quantity for unit activities 
of reactants and products 
Inasmuch as the solution is well buffered, the 
hydrogen ion formed by reaction 1 is largely utilized 
in converting acid succinate ion (HSuc~) into 
molecular succinic acid (H,Suc). The equilibrium 
between these constituents of the solution is the first 
dissociation of succinic acid: 
H.Suc=H* + HSuc”; A (4) 
The activity of hydrogen ion is fixed by the activities 
of the two succinate species and by the first dissocia- 
tion constant, A,, of succinic acid Hence, () can be 
written 


(/ A Ay suc 1 Busuc - (2a) 
As only electrical work is of concern here, 
AF LF, (5) 
and 
AF E F. 5a) 
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In eq 5a, E£° is the standard potential of reaction 1 
[3, 4] ?, and F is the faraday. Substitution of eq 2a, 
5, and 5a in eq 3 and rearranging gives 


> 
log kK, E , 3026) p eb lex Me |+ 


MH, Suc log Iu,suet cr ; 
log = T OR ’ () 
M usu I suc 


where / is an activity coefficient on the molal scale. 
In cell I, Moyo m 

If eq 6 is compared with the mass-law expression 
for the first dissociation (eq 4) in logarithmic form, 
it is seen that the portion of eq 6 enclosed in brackets 
actually represents the value of —log( fy fey) in the 
particular solution for which the emf, /, was meas- 
ured. For convenience, this experimental quantity 
will be called pwH. 

The last term of eq 6 cannot be evaluated exactly 
without a knowledge of A, However, in agreement 
with the predictions of the Hiickel equation [5], 
similar activity-coefficient terms for a number of 
weak monobasic acids have been found to be small 
and to vary linearly with ionic strength. By defini- 
tion, this term is zero at infinite dilution. Hence, 
log AK, at zero ionic strength was evaluated by 
extrapolation of log Aj, the value of log A, caleu- 
lated by eq 6 with the last term set equal to zero. 

The computation of the next to the last term on 
the right of eq 6 May be quite difficult for systems 
with two or more overlapping equilibria. The diffi- 
culty is aggravated if one group is moderately strong 
and the concentration of hydrogen ion is conse- 
quently an appreciable age tion of the ionic stre *ngth. 
Although the ratio of K, to K, for succinic acid is 
only 26 and pane re lg is considerable, neither 
dissociation constant exceeds 10~*, and the hydrogen- 
ion correction was readily evaluated. 

When sodium acid succinate is dissolved in a solu- 
tion of hydrochloric acid, as in the preparation of 
the buffer solutions used in this study, the reaction 
between hydrogen ion and acid succinate ion does 
not go to completion because of the small (but 
appreciable) dissociation of the molecular succinic 
acid formed (eq 4). Furthermore, the second disso- 
ciation constant is not small enough, nor the hydro- 
gen-ion concentration large enough, to prevent the 
dissociation of a part of the acid succinate ion, 


HSuc H*+Sue™; Ao. (7) 


It can be shown that the molalities of succinic acid 
and acid succinate ion remaining at equilibrium are 
given by the following equations: 


Mi suc m— My Meue™, (S) 
and 


Muysue 0.5 m+ MmMy—2 Meu” (9) 


The concentrations of the constituents of the cell 
solutions were so chosen that the concentration of 


? Reference [4] lists E° and 2.5026RT/F in absolute volts. 





succinic acid was about twice that of the acid 


cinate ion. In this way overlapping, as evidk 
by the magnitude of Msue, Was reduced. It sh | 
be emphasized, however, that this device is ot 
always advantageous, for my, which usually can vot 
be measured accurately, is inevitably increased 
Values of my, a quantity that approaches the 1) ve 
hydrogen-ion molality, my, as the ionic strenyth 
de ‘creases, were compute <d from pw H by 
log my =pwH—2Ayu TB yy). () 
The coefficient 7 in the denominator of the last term 
is the value found most suitable to describe the 
activity-coefficient term in mixtures of primary and 
secondary sodium succinates [1], and A and # are 
constants of the Debye-Hiickel theory. 
An exact expression for mg," is the following 


Msn Kymysuc~ (antilog pwH) (fasucfewfsuct). (11 


If the activity coefficients are again estimated by the 
the Debye-Hiickel equation, an expression for 
Mesuc™ 18 obtained: 

log Miysue 


log Mane = log k, t 


pwH+2A yu/(1+7By yp): l2 
Inasmuch as the dissociation constant is obtained by 
extrapolation to infinite dilution where uncertainties 
in the estimation of activity coefficients disappear 
My and Mgye" were used for my and mgy." in eq 8 
and 9. The ionic strength is given by 


- , , ) 
u L.5jm- My + Mguc™- 13 


The ionic concentrations required in eq 6 were 
accordingly calculated by eq 8, 9, 10, and 12 with the 
aid of the values of A, obtained in the earlier study 
[1]. The ionic strength needed to compute my and 
Mu" Was Obtained by arithmetic approximations. 
Both acid groups of succinic acid are fairly weak 
and my was only 1.6 percent of m for the most dilute 
solution studied. The extent of overlapping was 
greatest for the most concentrated solution, in which 
Moye" amounted to about 1.5 percent of m. 

Figure 1 illustrates the effect of the A, correction 
and of the method of estimating the activity coefli- 
cients on the determination of AK,. The circles and 
upper solid line (drawn by inspection) represent the 
final corrected values of pA, (—-log Aj). If the 
correction for bivalent succinate ion had been omit- 
ted, the lower solid line would have resulted, and the 
value of pK, found would have been about 0.02 unit 
too low. If 4 had been chosen instead of 7 as the 
coefficient of the Byy term in the denominator of 
eq 10 and 12, the extrapolation would have followed 
the broken line. However, the two upper lines con- 
verge satisfactorily, and the value of pA, at zero 
ionic strength would not have been altered. 

The slopes, —8, of the plots of pA{ versus wu wer 
measuredjfor each temperature and Bu added to A 
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cure 1. Extrapolation of pK; to zero ionic strength at 25 


broken line represents the best fit of the experimental values of pK if 4 

nstead of 7 in the denominator of eq 10 and 12 Similarly, the position 
ft ywer solid line illustrates the effect of omitting the correction for over- 
{the first and second acid groups 


computed from each experimental point. The re- 
sulting values of pA, were averaged. The mean pA, 
is listed in the fourth column of table 2, together 
with the average deviation of the 10 individual 
values from the mean at each temperature. The 
uncertainty in log A, is estimated to be less than 
+().003 unit. The dissociation constant is given in 
the last column of the table. The values of log A, 
used to compute the correction for overlapping are 
given in the second column, and log A, derived from 
a measurement of A,A, in the earlier study [1] 
appears in the third. The distribution of the in- 
dividual results for pAY on either side of the mean 
at 0°, 25°, and 50° is shown in figure 2. 

Earlier determinations of the first dissociation con- 
stant are summarized in table 3. No recalculation 
of any of the data has been attempted. Unfortu- 
nately, Jones and Soper [23], who also used the hy- 
drogen-silver-chloride cell but a different method of 
calculation, do not list their emf values. Hence, 
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PaBLe 2 Values of log Ky, log Ky, and K, from O° to 50 
Tem log K ~log A —- . 
perature (ref. 1 ref. 1 log K Ki X10 
( 

0 674 4. 2S 4. 2845-0. OO18 5. 20 

’ 5. 661 4. 2s 4. 263140. 0011 5 46 

10 5.049 4.24 4. 2449.0. OOO8 5. 70 

15 5. 643 4. 2%. 4. 2316-0. 0008 5.88 

20 5. 639 4. 22; 4. 2176.0. OOO8 6. 07 

25 5.636 4. Wh 1. 2006.+.0. 0009 6. 21 

tt) >. 642 4. We 4. 1980-0. OOO8 6.34 

5 5. 148 4.1% 4. 1914+0. 0013 6. 45 

0) 5. 4 4.1% 4 ISTR+O. OO12 6 50 

45 5. 670 4.19 4. 186040. OO15 6. 51 

WO 5. 680 4. 18s 4. 1863+0. OO10 6. 52 


TABLE 3. Summary of determinations of Log Ky 


* Method t log A 


Reference 


Conductance 


( 
1SS89, Ostwald [6] 25 4.177 
1891, Crum Brown and Walker [7] 25 4.17 
1891, Walden [8] 25 4. 168 
1904, Voerman [9] 2 4, 186 
| 0 $20 
1910, White and Jones [10] 2 4.184 
| 85 4.181 
1916, Boéeseken and Verkade [11} 25 4. 186 
1935, Jeffery and Vogel [12] 25 4.179 
Reaction rate 
1924, Duboux and Tsamados [24] Sucrose inver 73 4. 25 
sion 
1928, Dawson, Hall, and Key [25 I,-acetone 25 4. 201 
Emf, cells with liquid junctiot 
1924, Auerbach and Smolezyk [13] Q-cal 20 4.1/4 
1925, Britton [14 H-cal 1s 4.086 
1926, Mizutani [15] H-cal 18, 20 4.09 
1928, Simms [16 H-cal 25 4.195 
1928, Kolthoff and Bosch [17 H-cal . 1s 4.18 
1928, Gane and Ingold [18] H-cal 25 4.133 
1929, Olander [19] H-Q; Qat I8 - ; 20 
1931, Gane and Ingold [20 H-cal 25 1.19 
1934, Ashton and Partington [21 H-eal 25 +. 060 
1936, German and Vogel [22] Q-cal 25 1.10 
Emf, cells without liquid junction 
| 2 1. 162 
1936, Jones and Soper [23 H-AgCl 4) 4.135 
| 7 4, 132 
2 4. 218 
13] t 
1950, This investigation H-AgC] ° on 
; $19 
| uw +. Ist 
* H=hydrogen electrode, Q=quinhydrone, cal= calomel 


the difference of 0.045 at 25° and 0.054 at 50° be- 
tween the two determinations cannot be explained. 


IV. Thermodynamic Functions 
By the method of least squares, the following two 


equations were found to represent the values of log 
». 


Kk, given in table 2: 
log K,=A/T+ B+ CT=1206.25/T—3.3266 


0.011697 T, 14) 








and 
log K,=A*/T+ B*+C* log T 
2462.9/ T—47.9882+-17.755 log T. (15) 


The first has the form of the equation proposed by 
Harned and Robinson [26] to express dissociation 
constants as a function of temperature. It is based 
upon the observation that AF° is a quadratic function 
of 7, whereas the second equation presupposes that 
A/T° is a linear function of 7. By differentiation of 
eq 14 and 15 with respect to temperature, two ex- 
pressions for the temperature coefficient of log A, are 
obtained: 


d log Ky/d7T=1206.25/ 7*—0.011697, (14a) 
and 
d log K,/dT=2462.9/ T?—7.7108/T (15a) 
The corresponding second derivatives are 
& log K,/dT 2412.5/T (14b) 
and 
@ log K,/dT?*=7.7108/ 7? —4925.8/T (15b) 


The changes of free energy, (AF°), heat content 
(A/7°), entropy (AS°), and heat capacity (ACS) for 
the ionization process in the standard state were cal- 
culated by the following formulas 


AF 2.3026RT log Ay, (16) 
AH® =2.3026RT?(d log K,/dT), (17 
AS® =2.3026Rilog K, 
T(d log K,/dT)}, (18) 
AC, = 2.3026R(27T(d log K,/dT) 
T?(@ log K,/dT*)| (19) 


The values of these functions at 0°, 25°, and 50° ¢ 
computed from the two empirical representations of 
the relationship between log A, and T are compared 
in table 4. Cottrell and Wolfenden [27] have com- 
bined theircalorimetric data for AH° and its tempera- 
ture coefficient with the dissociation constant found 
by Jones and Soper [23] to calculate the thermody- 
namic constants at 25° C. Their results, labeled 
(C and W) are included in table 4. 

The two empirical equations, 14 and 15, evidently 
lead to appreciably different results for the thermo- 
dynamic functions at the extremes of the tempera- 
ture range. Unfortunately, the only basis for a 
choice between them is the fit of the “observed”’ 
data. The mean departure of the calculated log 
K, from the observed is 0.0005 unit for eq 14 and 
0.0011 unit foreq 15. Indeed, eq 14 usually provides 
the better fit of the observed dissociation constants 





TABLE 4. Comparison of thermodynamic functions 
from two different empirical equations (eq 14 and 15 
variation of log K, with T 





F quat ion AF SH AS Ac, 
ec 
mole mole deq™' mole deg-' mole 
14 22, 407 6, 384 5B OG 22 
15 22, 417 6, 828 57.1 148 
25° ¢ 
14 24, 013 186 69.8 ’ 
l 24, O07 138 70.0 148 
C and W « 23. 76 2 519 71.1 132 
we ¢ = 
14 25, 900 203 81.0 4 
15 25, GOS 553 81.9 148 


* From the value of AK, found by Jones and Soper [23 


for weak electrolytes. Hence, eq 14, 14a, and 14) 
were used to compute the thermodynamic quantities 
listed in table 5. The estimated uncertainty i 





AF® is 18 j; in AH®, 300 j; in AS®, 1.0 j; and in 
4U,°, 15}. 
Tapie 5. Thermodynamic functions for the first disso 
of succinic acid 
Femper AF sH AS AC 
atur , 
c mole mole deg mole deg mole 
0 22, 407 6, 38 5.6 122 
22. 706 5, 767 “0.9 125 
10 2, O16 138 63.1 127 
15 23, 337 4, 499 65.4 120 
1) 23, 669 4, S48 67. ¢ 131 
25 + 186 690.8 134 
0) 2, 513 72.1 136 
AS 1, 828 74.3 138 
0 1, 132 76.6 140 
45 425 7AS 142 
sw 203 S1.0 14 m 
es 
pa 
Cottrell and Wolfenden [27] estimated that A// of 
becomes zero at about 45° C. The dissociatio ere 
constants found by Jones and Soper [23], on th my 
other hand, appear to reach a maximum (and —log Th 
AK, a minimum) between 50° and 74°. The tempera- to 
ture, Tim, at which log K, passes through th: da 
minimum is readily calculated in terms of the con me 
stants of eq 14 and 15 by setting the first derivativ: me 
equal to zero. Hence, from eq 14, th 
es 
) an 4 1 ( d ) 
and from eq 15, aul 
tu 
7" eo Game " - vO 
Tot 2 3026 A* ( * - +} 


The data of this investigation place the minimu! 
value of —log K, at 47.9° C (eq 20) and 46.2° ‘ 
(eq 21). 
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Effects of Prior Static and Dynamic Stresses on the 
Fatigue Strength of Aluminum Alloys 


By John A. Bennett and James L. Baker 


Tests made on specimens of Alclad 248 


T sheet showed that prior static load had a 


marked effect on the fatigue strength in unidirectional bending when the stress amplitude 


was relatively small. 


From tests on bare 24S 


T sheet, it was found that a few cycles at a 


stress amplitude of 17,000 |b/in.? resulted in a large increase in the fatigue life at 20,000 
lb/in2 Damage tests for other combinations of stress amplitudes indicated that the damage 
was nearly a linear function of ratio of the number of cycles at a given stress to the number 


that would cause failure at that stress. 


\ new design of specimen and a new form of stress 


versus number of cycles to fracture (S-.V) diagram are described 


I. Introduction 


The phenomenon of fatigue failure in metals re- 
mains incompletely understood, despite the extensive 
research that has been done on the subject. In 
particular, the mechanism by which a small volume 
of metal is damaged by repeated stressing until a 
crack forms may be considered one of the outstanding 
mysteries in the mechanical behavior of materials. 
The purpose of the work reported in this paper was 
to evaluate some of the factors involved in this 
damaging process by using aluminum alloy speci- 
mens As there appears to be no reliable means of 
measuring fatigue damage except by a fatigue test, 
this evaluation was carried out largely by means of 
tests at two stress levels. 

Many investigators have reported experimental 
data showing the effect of dynamic stressing at one 
auplitude on the fatigue strength at a second ampli- 
tude, but the data for aluminum alloys are not in 
good agreement. Miner [1] ' found in testing 24S-T 
that the damage was proportional to the cycle ratio? 
within the experimental error. On the other hand, 


res in brackets indicate the literature references at the end of this paper 
ratio is defined as the ratio of the number of cycles at a given stress 
¢ to the number necessary to cause failure at that amplitude. Failure 
defined as either the beginning of a crack or complete fracture 


Stickley [2, 3] reported that the fatigue strength of 
17S-T rod could be increased by prestressing at a 
stress amplitude slightly below the test stress. Like- 
wise, Work & Dolan [4] found an increase in fatigue 
life due to stressing at small amplitudes. However, 
Dolan, Richart, and Work [5] found for 17S-T that 
the fatigue life under gradually varying load cycles 
was greater when the minor stress amplitude was 
18,000 Ib/in.? than when it was decreased to 16,000 
lb/in.? 

These results do not give any consistent picture 
of the progress of fatigue damage in aluminum 
alloys, and it was the purpose of this work to evalu- 
ate this progress in much the same way that it was 
done for SAE X4130 steel in the investigation re- 
ported in [6]. 

One of the most interesting phenomena associated 
with fatigue is the improvement due to understress- 
ing. In ferrous materials this improvement increases 
as the amplitude of the understress approaches the 
fatigue limit. Aluminum alloys have no fatigue 
limit, and the effect of understressing is not so 
marked, but it is known to cause an increase in fatigue 
strength under certain conditions. It seemed rea- 
sonable that during the course of a fatigue test at 
relatively low stress there must be a reversal in the 
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type of process taking place, with improvement 
occurring during the early part of the test, and 
damage later. In order to investigate this, one of 
the prestress amplitudes chosen was low enough so 
that the number of cycles to fracture would be very 
large. 

When a metal is subjected to fluctuating stress, 
there are presumably two opposing processes taking 
place, namely, work-hardening and fatigue damage. 
While the mechanism of fatigue damage may involve 
repeated plastic deformation of a microscopic volume 
of the metal, this is apparently distinct from the 
process that is usually referred to as work-hardening. 
It is to be expected that the greater part of this over- 
all hardening would be completed in a relatively few 
cycles at the start of the dynamic stressing. For 
this reason, it was decided to make tests after a very 
small number of cvcles in order to detect changes due 
to work-hardening 

A few tests were made to determine the changes in 
fatigue properties caused by static stressing, prior to 
testing, for comparison with the changes caused by 
dynamic stressing prior to testing. In order to make 
the effect as pronounced as possible, the tests were 
made on Alclad 248-T; the soft cladding would be 
expected to deform plastically at a lower stress than 
the core. producing a residual that would 
influence the fatigue properties 

Fatigue tests usually are performed by applying 
the fluctuating load until the specimen fractures. 
The life of the specimen, as determined by such a 
test, includes both the period of increasing damage, 
which culminates in the formation of a crack, and 
the period of crack growth, which ends with the 
fracture of the specimen. In order to study the 
first period without having the results confused by 
the second, it is necessary to determine as nearly as 
possible the point at which the crack forms. 


stress 


II. Materials and Test Methods 


All of the tests reported here were made on Krouse 
sheet bending fatigue testing machines. These 
machines are designed to load the specimen as a 
cantilever, the range of deflection being determined 
by the setting of an adjustable eccentric and the 
mean load by the adjustment of the movable vise. 
The amplitude of deflection remains constant during 
the course of the test, so that any change in the 
resilience of the specimen results in a change in the 
applied stress. 

The materials used in this investigation were bare 
and Alclad 245-T sheet. For the tests reported in 
sections II] and IV, the materials were obtained 
commercially. For the study of damage reported 
in section V, the 248-T sheet was obtained from the 
NACA Langley Aeronautical Laboratory from a 
stock maintained especially for fatigue research. 

The two types of specimens used are shown in 
figure 1, a and b. The first is the standard type of 
specimen recommended by the testing machine 
manufacturer, while the second is a modification 
adopted for the reasons givenin section IV. The 
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over-all dimensions of the two specimens we 
same, and the methods of test were identical 
specimens were cut with the long dimension p 
to the direction of rolling. 

There is a possibility of confusion of ter 
referring to the various surfaces and lines that | 
a specimen cut from sheet material; in this 
the term face will be used to refer to the surfa: 
the original sheet. The surfaces at right ang 
the faces and roughly parallel to the long dim 
of the specimen will be called the sides, and the 
section of a face and a side will be termed an 
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Fiat RE | Types of Specimen used for he nding falig rT, 
sheet materials. 
i, Standard specimen; b, modified specimen 

The edges of all specimens were manually polished 
with No. 400 Aloxite paper, the direction of polishing 
being parallel to the edge. 

In both types of specimens the maximum stress 
the reduced section was only about 25  peret 
greater than that at the point where the specime! 
was clamped in the vise. For this reason, ther 
was some difficulty at first with the specimens 


breaking at the vise. This trouble was 
entirely eliminated by placing several thickne 
greased paper between the specimen and the « 




















ng plates to distribute the load. However, it is felt 
hat in any future work with these machines it 
ould be advisable to increase the width of the 
pecimen at the vise to eliminate the need for this 
precaution. 

The load necessary to produce a given stress in 
he specimen was calculated from the usual canti- 
ever formula 

P Sd% 
6L 


where P is the load necessary to produce a given 
stress in the specimen; S is the nominal stress in 
the extreme fibers; d is the thickness of the speci- 
men; 6 is the width of the specimen at a given point; 
and ZL is the distance from the point where the width 
is measured to the point at which the load is applied. 
The standard specimens had a reduced section that 
tapered toward the crank pin where the load was 
applied. This resulted in a constant nominal stress 
in the reduced section, as both the width of the speci- 
men and the bending moment are proportional to 
the distance from the point of application of the load. 

The specimens of 4-in. radius had no straight side, 
so it Was necessary to use a special jig to determine 
the minimum value of the ratio 6/Z in the above 
formula. The jig is shown in figure 2. The pivot 
on which the vanes turn was at the same distance 
from the end of the specimen as the crank pin of 
the machine, and the extended edges of the vane: 
passed through this pivot. Therefore, if the vanes 
were rotated until they were tangent to the edges 
of the specimen, the point of tangency was the point 
of minimum 6/Z ratio. In use the jig was clamped 
to the table of a toolmaker’s microscope and the 
length of the jig alined with the longitudinal measur- 
ing screws. Then the vanes were rotated until only 
a narrow gap remained between each vane and the 
corresponding edge of the specimen. The width of 
the specimen was measured near the center of this 
gap by using the transverse measuring screw, and 
the distance to the center of the pivot was measured 
with the longitudinal screw. In the figure a broken 
specimen is shown in place on the jig to indicate 
how the typical fracture occurred near the point of 
minimum 6/L ratio. 


III. Effect of Prior Static Load on the Fatigue 
Strength of Alclad 24S-T 


The material used in this investigation was 
obtained commercially, and the surface was badly 
scratched and dented. These defects, however, did 
not appear to influence the fatigue test results. All 
of the tests were made in unidirectional bending with 
the mean load equal to half the maximum. The 
stress values given in this section refer to the maxi- 
mum nominal stress in the cycle. The standard 
'vpe of specimen (fig. 1, a) was used. 

In order to obtain the basic S-N curve, 24 speci- 

ens were tested at the following nominal stresses: 

000, 25,000, 30,000, and 35,000 Ib/in These 

sults are shown in the first line of table 1 and are 





Figure 2 Four-inch radius specimen in place on measuring 


jig 


\ broken specimen is shown in order to illustrate how the fractures occurred 
near the point of minimum b/L ratio 


represented by the curve (not the points) en figures 
3 and 4. 

The effect of prior static load was then determined 
for four conditions of loading, as follows: (1) static 
load corresponding to a nominal stress of 30,000 
Ib /in.? applied in the same direction as the subsequent 
fatigue test load; (2) the same load applied in the 
direction opposite to that of the subsequent fatigue 
test load; (3) static load corresponding to a nominal 
stress of 40,000 Ib/in.? applied in the same direction 
as the subsequent fatigue test load; and (4) the same 
load applied in the opposite direction. In each case 
the prior static load was allowed to remain on the 
specimen for approximately 1 minute, then the 
fatigue test was started as soon afterward as pos- 
sible. In most cases two specimens were tested for 
each load condition, the values given in table 1 being 
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CYCLES OF STRESS 
Figure 3 Effect of prior static stress (nomina maximum 


value 30.000 Ib/in2 on the fatique properties of Alclad 
24S-—T sheet in unidirectional bending 


Solid curve represents data for original material; squares are averages for 
specimens prestres in direction opposite to test stress, triangle ire average 





for specimens prestressed in same direction as test stress 
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Fatigue test in unidirectional bending 


Minimum stress 


Effect of prior static stress on the fatigue stre ngth of Alclad 24S-T sheet 


0, maximum nominal stress is listed 


Average number of cycles to fracture at test stresses listed 
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1), OOO Surne 4 10. 7F 6, 628 1 
10, 000 do wn 2, SAS 2 
0,000 | Opposite 2, 788 1,024 2 
0. 000 do 2 2, 166 45 2 
* Direction of prior static load relative to fatigue test load 
35 
4 SAME 
n O OPPOSITE 
2 so} 
= 30}+— | 
n 
2 
°o 
= 
- 
of 20)... 
= 
Q 
4 
no 
no 
ae 
x 
= 20}. 
n 
2 
3 '° ‘ + 
= 
4 
2 
‘ a 
0 10 10 
CYCLES OF STRESS 
Ficure 4 Effect of prior static stress (nominal maximum 
value 40,000 Ib/in2) on the fatigue properties of lelad 


294S—T sheet in unidirectional be nding 


Solid curve represents data for original material; squares are averages for 
specimens prestressed in direction opposite to test stress, triangles are averages 
for specimens prestressed in same direction as test stress 


the averages of the number of cycles to failure for 
the specimens tested under identical conditions. In 
figure 3 the results for the specimens tested after a 
static stress of 30,000 Ib/in.? are shown as individual 
points for comparison with the curve representing 
the fatigue properties of the original material, while 
the points in figure 4 represent the data for a prior 
static stress of 40,000 Ib/in It should be empha- 
sized that the stress values given here are calculated 
on the basis of elastic bending of the specimen. As 
there is considerable plastic deformation, particu- 
larly at the higher loads, these nominal values of 
stress are undoubtedly considerably different from 
the true stress. However, these values do indicate 
the magnitude of the bending moments applied to 
the specimens, which are significant from a practical 
stand point. 

The results indicate that there was little or no 
effect except at the lowest test stress used (20,000 
lb/in.2). At this stress (20,000 lb/in.*) the effect was 
quite marked; for example, with a prior static 
stress of 40,000 lb/in.* the specimens loaded in the 
same direction as the fatigue test load endured 10 
times as many cycles before fracture as those loaded 
in the opposite direction 








i) Range of No. of 30,000 Range of | No. of 35,000 Ra 
1 values tests Ib/in values tests Ib/in 

10 1.002 « 10) . 401 10 2?) 10 f xX 10 a) 
27 

20 2s4 2 wT i 222 

Avs 75 2 286) 1h 1 207 

“77 ‘ 2 OS lay 2 25 


IV. Effect of Specimen Geometry on Results 


of Bending Fatigue Tests 


As a result of the tests reported in the preceding 
section and other bending fatigue tests made pre- 
viously, it was thought that the form of specimen 
commonly used in this type of machine had certain 
disadvantages. Although it is designed to have a 
section in which the stress is nominally constant, this 
section is bounded by fillets that raise the actual 
stress sufficiently to cause most of the failures to 
originate at the bottom of one of the fillets. Accord- 
ingly, it was thought better to make the specimens 
with each side formed by a single circular are of larg 
radius, and the design finally adopted is shown 
figure 1, b. The reasons for selecting the value of 4 
in. for the radius were as follows: (1) This is the 
largest radius that could be used with essentially th 
same over-all dimensions as the standard specimen 
(2) it is sufficiently large so that the stress concentra 
tion due to the radius is negligible; (3) a sheet speci 
men having the reduced section formed by an are 0 
4-in. radius on each side is being used for axial load 
fatigue tests in the Engineering Mechanics section o! 
this Bureau, and it was thought that future compa: 
sons of data would be simplified by having a similar 
reduced section on the sheet bending specimens 
(4) it was desirable to have a radius that could lx 
used also for R. R. Moore and axial load fatigu 
specimens, and new type specimens for these ma 
chines have been designed with a 4-in 
(5) 8-in. diameter milling cutters are readily availabl 
so that the specimens could be milled with a sing! 
traverse of the cutter on each side. This factor ts 
important in reducing the cost of the specimens 

In order to compare the results obtained with th: 
two designs, 20 specimens of each type were prepared 
from a sheet of 248-T nominally 0.!02-in. thick. The 
results of fatigue tests made on these specimens ar 
shown in table 2 and figure 5. All tests were mad 
completely reversed bending. Three specimens 0! 
each shape were tested at each of the four tes! 
stresses to determine the approximate relationship 0! 
the fatigue strengths, then seven additional specimens 
of each design were tested at a stress amplitude 0 


f 


radius 


{ 


’ To avoid confusion it should be pointed out that the term ampilit 
here to mean the value of A in the equation for the stress at any time 


S=A (sin 2eff+B, 


where f is the frequency of the testing machine, and B is the mean st 
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000 Ib/in2 in order to obtain a more accurate 
mparison of the dispersion of the data. 

The fatigue strength of the 4-in. radius specimens 
is consistently higher than that of the standard 
ecimens, although the difference was not large. 
1e stress to cause fracture in a given number of 
cles was, on the average, 
ew type specimens than for the straight tapered 
the difference being greater the larger the 


‘ cs, 


stress amplitude. 


This diffe 


nA 


5 pereent higher for the 


‘rence might be ascribed 
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CYCLES OF STRESS 
Ficure 5. Comparison of fatigue properties of 24S-T sheet 


obtained with standard and 4-in. radius specimens. 


TABLE 


9 
“-- 


1. Effect of specimen shape on fatigue strength of 24S—T sheet 


Average number of cycles to fracture 


Stress 
ampli 
tude 


thiin2 
20, 000 
25, 000 
$0, 000 
$5, 000 


Standard specimen 


4-inch radius 
specimen 


A verage Average 
Number numberof Number number of 
of tests cycles to of tests cycles to 
failure failure 
; 5 STIX 1 ; 6, 44x10 
10 145 10 743 
3 146 ; 536 | 
; 13 $ 132 


B. Dispersion of data for a 


Standard deviation, «* 
Coefficient of variation, 100e0/ Nr 


stress of 25,000 Ib/in2 


4-inch 
radius 
specimen 


Standard 
specimen 


Ne 198x105 743 X10 
19 289 
(%) Os 30 


values of standard deviation given here were computed by using the 


Ja given on p. 14 of [9] 


ore common form of the equation is 


xX 


. ¥)2 
l 


s difference should be kept in mind in comparing these data with those 
ted by other investigators 


either to the lower stress concentration in the 
4-in. radius specimen or to an effect analogous 
to size effect, considering that the new specimen has a 
smaller volume of metal to which maximum stress is 
applied. If the latter factor were the chief considera- 
tion, it would be expected that the dispersion of the 
results for the tapered specimens would be less than 
that of the 4-in. radius type. However the coeffi- 
cients of variation in table 2 indicate that the oppo- 
site is true, and that there is quite a difference in 
favor of the latter type. On the basis of these re- 
sults it was concluded that the 4-in. radius specimen 
was a more satisfactory design than the tapered one, 
both as regards elimination of stress concentration 
and reduction of scatter of the results. Accordingly, 
specimens of this type were used in all of the subse- 
quent work reported in the following sections. 


V. Evaluation of Fatigue Damage in 


24S-T Sheet 


It was intended that this study of fatigue damage 
would follow approximately the same course as that 
described in the latter part of [6]. However, the 
data proved to be much less reproducible than the 
R. R. Moore tests on steel, so the results could not 
be so readily analyzed and only limited general con- 
clusions could be drawn from them. 

The determination of the point at which a crack 
starts in a fatigue specimen requires two techniques 
not ordinarily employed in fatigue testing, namely, 
stopping the test when a small crack has formed, and 
measuring the growth of the crack with increasing 
numbers of cycles. This provides information for 
estimating the number of cycles between the start 
of cracking and the first observation of the crack. 
The method used for stopping the test when a small 
crack had formed was that described by Foster [8]. 
This involved cementing #® small diameter wire on to 
the surface adjacent to the critical sections of the 
specimen. The wire carried the operating current 
to a relay controlling the testing machine motor, so 
that when the wire was broken the machine would be 
stopped. It was reported in the reference cited that 
the wire did not break until a crack had formed under 
it. This was not found to be the case with the sheet 
specimens in bending, although the conditions re- 
ported by Foster were duplicated as nearly as possi- 
ble. Wires cemented onto the face of the specimens 
broke in as little as one-fourth the number of cycles 
required to start a crack in the specimen. This 
caused a considerable delay for inspection of the 
specimen and cementing another wire in_ place. 
Accordingly, the wires were usually cemented onto 
the side of the specimen, a short distance away from 
the edge where the strain amplitude was not so great. 
This had the disadvantage that a somewhat larger 
crack was required to break the wire than when it 
was on the face, but the sensitivity was considered 
adequate. 

After a small crack had been detected with the 
fracture wire, measurements of the growth of the 
crack were made at intervals until final fracture. 
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Figure 6 shows typical curves of crack length (ex- 
pressed in terms of the ratio of the crack length to 
the width of the specimen at that point) versus 
number of cycles for three specimens tested with a 
nominal stress amplitude of 30,000 Ib/in The crack 
length was measured on the face of the specimen 
that had the longest crack, and the differences in 
the shapes of the curves in figure 6 are presumably 
due to differences in the development of cracks on 
the other face. The crack growth curves were 
extrapolated back to 0.01 crack-length ratio, and 
this was arbitrarily taken as the start of the crack. 
While this is obviously a rather approximate process 
for determining the start of cracking, the accuracy 
was considered adequate in view of the fact that the 
number of cycles from cracking to fracture was small 
in comparison to the number of cycles required to 
start the crack 

Several specimens were tested in the manner 
described above to determine the average number of 
evcles from cracking to fracture (Ne—Ne) (where 
Ne is the number of cycles at which the specimen 
fractured completely, and N¢ is the number of cycles 
at which a crack first formed) at each of four stress 
amplitudes, 20,000, 25,000, 30,000, 35,000 Ib/in.? 
These values were then subtracted from the average 
number of cycles to fracture for all specimens tested 
at the corresponding stress amplitude in order to 
arrive at the average number of cycles to form a 
crack, Ne. Both of these sets of data are shown in 
figure 7 and also are given in table 3. Three tests 
were made at a stress amplitude of 18,000 lb/in.* in 
order to obtain a better idea of the general shape of 
the S-N curve, but no determination of the start of 
cracking was made at this stress, as it was not to be 
used as a test stress in the evaluation of damage. 

The evaluation of fatigue damage was carried out 
by stressing a specimen for a predetermined number 
of cycles at one stress amplitude, then testing it to 
fracture at a second amplitude. The damage was 
defined as the percentage decrease in the number of 
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Crack growth curves for 24S-T specimens tested at 
a nominal stress amplitude of 30,000 lb/in2 


Figure 6 


Crack length ratio is equal to the crack length divided by the specimen 
width 




















TABLE 3. Data for S—N ¢ urve, 24 S—T' sheet 
A verage Calculated 
numberof Num Ay . number of 
Stress Number cyclesfrom ber of Bao ee ~y , cycles to pr 
umpli- (of speci-| beginning | speci om Ag duce crack 
tude mens of crack to mens digg « 
fr samen tested fracture, Ne Nt Nr 
Ne— Ne Ne—N 
lbjim2 
$5, 000 ry 27 x10 iF) 141 «10 1l4xK ht 
$0, 000 ) 70 7 $36 ve 
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CYCLES oF STRESS 
Figure 7. S-N curves for cracking and fracture in 248 


sheet specimens. 


cycles required to form a crack in the prestresse: 
specimen as compared with the original material, 


Ne—N 


CP 
<100, 
Ne 


D 


where )) is the percentage damage, and Nep is thy 
number of cycles at which a crack first formed in a 
prestressed specimen. None of the tests for the 
determination of N@ (calculated number of cycles 
to crack) had been made at the stress amplitudes 
chosen for the prestress, so the number of cycles to 
start of cracking had to be obtained by interpolation 
To make this interpolation as accurate as possible, 1! 
was desirable to find a method of plotting the data 
that would give a straight line. As shown in figure 8 
this condition was satisfied very closely by plotting 
the logarithm of N? against the reciprocal of th 
stress, so the interpolation was performed on this 
basis. 

Three prestress amplitudes were used, 17,()00, 
22,500, and 32,500 lb/in®. For the lowest value, th 
number of cycles to start a crack was not known, se 
the numbers of cycles of prestress were arbitrarily 
selected as 2 10°, 2*10° and 2X10". For the two 
higher stresses it was intended to apply the prestress 
for cycle ratios of 0.1, 50, and 90 percent. How 
the number of cycles corresponding to the inte: a 
values of cycle ratio were determined on the basis 0 
a group of tests made prior to starting the two- 
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, _ - T > TABLE 4 Summary of damage data, 24S—T sheet 
0 a SE ST SS EN EN See ee ee 
| 2s Se Ges a 
O ALL SPECIMENS 2X 1 cycles 2 10 cycles 2x10" eveles 
| prestress prestress prestress 
@ UPPER 60 PER CENT iaiaie 
OF SPECIMENS Number Number Number 
: : ; of spec Damage of speci- Damage’! of speci- | Damage 
r mens mens Inens 
Prestress amplitude 17,000 Ib/in.? 
6 hii Percent Percent ’ercen 
F190 20,000 2 448 427 ’ ¥ 
- 25,000) 65 i) 2 
a 10.000 2 ; 43 
‘ 65000 14 a 11 
oO 
> 
oO 
Prestress amplitude 22,500 Ib/in.2 
u 
© 
0.1-percent cycle ‘4-percent cyck S0-percent evclk 
a ratio ratio ratio 
W 
: 5 20.000 i ; f 77 2 42 
> 10 25.000 ; 17 4 73 2 +A 
z 30,000 ) ; 62 F 141 
Qa02 0.03 0.04 0.05 , 35.000 4] 7 67 
. . 
RECIPROCAL STRESS, (LB/IN, THOUSANDS) 
Prestress amplitude 32,500 Ib/in.? 
Fy rE 8. Graph illustrating the linear relationship obtained 
hen log N; is Pp otted against the rec tproc al of the stre 5s 0.1 y™ reent cv¢ le 42 pt reent cvy¢ le 75 perce nt cy | 
nplitude ratio ratio ratio — 
stress-level tests. Single-stress-amplitude tests were saan - 2 - ; p. 
made of additional specimens after the damage 30,000 14 9 2 k.. 
— bof } 25 ’ Mm 
? tests had been started, and these changed the 
average values of N*, so that the final values of 
; evele ratio were somewhat different from the in- 
' tended ratios. The value of N? originally obtained 
at 32,500 lb/in.? was too low, so that the final cycle 
se ratios are smaller than the intended values, whereas TT T 
the reverse was true at 22,500 lb/in*®. In addition, Q 
the highest ratio at 22,500 lb/in*. was altered as a w» 100 + y 
result of specimens that failed before the prestressing o |e 
was completed, as explained below. The final inter- 4 aE ae a ee oe wi 
wlated values of N22 which were used in deter- 
lated val f N* whicl let So | 7 
~~ $¢- —4——+—_ + —+ 
the mining cyele ratios were | | Pp 37 | 
na wn ’ T t -_ T T 
the V* 108 WwW | | | | Ss 
les Stress amplitude ae ( = 60 oO T =. me “— Game geen T a 
. cvcles Oo | | 
aes — | “i | 
to < ae -_ —I 
ol hiinr | 7 
on ibjin- | | 
if 22,500 1, 650 = 40 | mr | T T 
32.500 175 < / | 
ata = + adh ae ee | | aa a 
es < 4 Ry | 2 
ing ; = sel—) ate __| TEST STRESS, LOIN 
the lhe results of the damage tests are shown in figures - y @ +20,000 
his 9, 10, and 11 and are listed in table 4. A negative z -——1—_+—— § £26000 —— 
value of damage indicates that the average life of S / o 30,000 
wy) the specimens that had been prestressed was greater rw) 0 a +35,000 — 
. than that of the original material at the same stress a 
ry : , =e +—_—__—__+—____+— 
s amplitude There are two values of damage that rot | | 7 T 
ily exceed 100 percent; these are possible because the -eoL | | | 
Me calculation of damage is based on NZ. If the average 0 20 40 60 80 100 
ss value of Np for the prestressed specimens is less than PERCENT CYCLE RATIO aT ¢ 2 
the number of cycles from start of cracking to frac- = 52,500 LB/IN. 
d l then the average Nep will be negative, and the Figure 9. Evaluation of fatigue damage due to prestressing at 
of ( lated damage will exceed 100 percent, an amplitude of 32,500 Ib/in?, 
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100 RE ee eyeles at +22,500 lb/in.? introduced an apprec ab}; 
7) TEST STRESS,LA/IN | [7 error, as these results are quite out of line wit!. th, 
= we ted ae ee : others for this prestress amplitude. 
= 80 pp tne | |} / The effect of prestress at an amplitude of 17 (09 
wo a +35,000 E r. en lb/in.? was more pronounced, and the scatter t)ore- 

ratenny ae i = fore less important. As shown in figure 11, r 
r | Ps 4 was an increase in fatigue life as a result of this 
w 60 +—7—_+— +2 +——1—J prestressing, which was particularly marked at hy 
- Fa lowest test stress. The effect at all of the other (es; 
a » oe or ee me stresses was of approximately the same magnit uc 
a on oe Fi | =e eS so that only a single curve has been drawn on thy 
je Kr Y oO figure to represent these data. 
© + +— At + +—_—___4___}- — 
3 / VI. Discussion of Results 
at 20 pn + + + }—_—__+ - ——4 
ee ee ™ 
“si SS Te . 1 J The results of the tests described above indicat 
- / that both static and dynamic stress applied prior to 
w O 4 + the start of a fatigue test can have a marked effec 
° : | on the fatigue life at relatively low stresses. In th 
“ ps t a es + Fava — case of the Alclad material tested in unidirectiona 
-20 | | } | | | be oC - a a = — —_—— wer 
probably largely due to residual stress in the oute: 
stnneel » TR, Oe = s nil ms - fibers of the sheet. In the case of the bare material 
, tested in reversed bending, the improvement caused 


Fiaure 10. Evaluation of fatigue damage due to prestressing by a small number of cycles of stress slightly below 
at an amplitude of 22,500 lb/in’. the test stress was probably caused by cold-working 
In both cases it appeared that the beneficial effect 
For the tests involving high cycle ratios, some of was eliminated by the fatigue stressing at all amp! 
the specimens broke before the prestressing was tudes above the lowest value used in the tests. 
completed. The corrections applied to the data to The dispersion of the data was so great that 
account for these failures are explained in the would have been impossible to obtain results as pre- 
appendix. cise as those reported prev iously for steel [6] without 
It can be seen from the figures that the data were using an unreasonable number of specimens. It is 
somewhat scattered. No attempt has been made to not known how much of the difference in the scatte: 
draw curves illustrating the trends of the data on of the data was due to inherent difference in the ma 
figures 9 and 10, but a dashed line at 45° (represent- terial and how much was caused by the specimer 
ing a constant rate of damage) has been put on the geometry, but certainly the latter factor would b 
figures for reference. While most of the points lie important. It is much more difficult to obtain 
above this line, there is no consistency as to the producible conditions on the edge of a sheet than o1 


















relative position of the points representing the differ- the surface of a cylindrical specimen. 
ent test stresses. It appears also that the method of 
computing the cycle ratio for the largest number of VII. Conclusions 
” As a result of flexural fatigue tests made on 245-T 
e aluminum alloy specimens the following conclusions 
= were reached. 
4 1. The application of a static bending load to th 
> specimens had a significant effect on the fatigue lif 
- of the specimens under unidirectional bending, if th 
_ | | | | | ; stress amplitude in the fatigue test was small. For 
= a large static load in the same direction as the fatigue 
w -200 t t t load, the fatigue life was increased, while the revers’ 
> 2 was true if the directions were opposite 
3-300 sence te 2. A new design of specimen, having each sid 
0 + 25,000 formed by a circular are of 4-in. radius, was foun 
9-400 ae to be more satisfactory than the straight tapere 
w _—? ; design. The new specimens were cheaper to | 
-~500 facture, had less stress concentration in the tes! 
» 2s 5 6 © @ 4 6 8, 20 section, and gave more-reproducible results than t! 
MILLIONS OF CYCLES OF STRESS AT + 17.000 LB/IN standard type. 
Ficure 11. Evaluation of fatigue damage due to prestressing 3. Fracture wires were found to be a satisf: mY 
at an amplitude of 17,000 Ib/in? means of stopping these tests when a small fatigw 
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ck had formed, but their use involved some 

iculties. 

i. The dispersion of the results of flexural fatigue 

ts on aluminum alloys was much greater than 

t of rotating beam test results on steel previously 

vorted. 

5. The relation between the reciprocal of the stress 

d the logarithm of the number of cycles to the 

irt of cracking was found to be linear. 

6. For prestress amplitudes of 32,500 and 22,500 

in.2, the precision of the data was not sufficient 
to justify an exact statement regarding the progress 

f fatigue damage. The assumption of a linear rela- 
tion between cycle ratio and damage was within the 
experimental error. 

7. Prestressing at a stress amplitude of 17,000 
lb/in resulted in an improvement in fatigue life, 
that is, the damage was negative. For specimens 
subsequently tested at a stress amplitude of 20,000 
lb/in.*, the improvement in fatigue life after either 
2,000 or 2,000,000 eycles of prestress was more than 
100 percent 
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Vill. Appendix 


1. Correction for Premature Failures 


The proportion of premature failures for the four conditions 
in which they occur are 


Cycle 
ratio or 
number of 


Prema- 
ture 
failures 


Stress amplitude 


cycles 
lb/in? Percent 
17,000 2x 10° 14 
22,500 54AQ 14 
22,500 ‘ 98% 10 
32,500 75% 7 


For the first, second, and fourth conditions listed, the data 
were corrected for the premature failures by considering that 
these specimens had a value of Ne=0 at the test stress As 
the values were actually less than zero, this assumption had 
the effect of making the average N- at the test stress too high 
However, for the conditions mentioned, where the proportion 
of premature fractures was small, this error was not considered 
significant 

For the third condition in this tabulation, the proportion of 
premature failures was so large that the method used in 
handling the data for the other three conditions was not 
considered satisfactory Accordingly, a new value of cycle 
ratio was determined from the following considerations; as 40 
percent of the specimens broke before the prestress was 
completed, those that were tested might be considered to be 
the strongest 60 percent of the population and the average 
V- for this portion is larger than that of the whole popula- 
tion. This new value, N;, was determined by interpolation 
between the values for the top 60 percent of the specimens 
tested at + 25,000 and + 20,000 |b/in upper line on fig. 8 
On this basis the eyele ratio of the specimens tested was 80 
percent, 


WASHINGTON. March 27, 1950. 
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Boiling Points of Aqueous Solutions of Dextrose 
Within the Pressure Range of 200 to 1,500 Millimeters 


By John L. Torgesen,' Vincent E. Bower, and Edgar R. Smith 


By using twin ebulliometers of the Swietoslawski type, with water as the reference 
standard, data were obtained on the vapor-pressure—temperature relationship for aqueous 


solutions of dextrose 


The pressures ranged from 200 to 1,500 millimeters and the concen- 


trations of the solutions from 10 to 60 percent of dextrose by weight. The Vapor-pressure 
temperature relationship is expressed by equations of the form 


logio Pp 


A 


B 
C+t 


in which p is the vapor pressure in standard millimeters of mercury exerted by a solution of a 
given concentration at the temperature, ¢, in degrees Celsius. A, FP, and C are constants 
The boiling-point elevations of aqueous dextrose solutions at a given pressure are expressed 


by equations of the form 


logy At=a-+t 


Bc + yc? +6c3, 


in which At represents the boiling-point elevation in degrees Celsius at a given pressure, c 
is the concentration in weight percentage of dextrose, and a, 8, y, and 6 are constants 


I. Introduction 


The boiling points and boiling-point elevations of 
aqueous solutions of dextrose within the pressure 
range of 200 to 1,500 mm and the concentration 
range of 10 to 60 percent of dextrose by weight are 
reported in this paper. This work is part of a pro- 
gram sponsored by the Corn Industries Research 
Foundation to obtain physical data for materials of 
importance in the manufacture of various corn 
products. 

The method used was a comparative dynamic one 
for measuring successively, with the same ther- 
mometer, the boiling point of a given solution and 
that of water in twin ebulliometers connected to a 
manostat. <A series of corresponding boiling points 
of the solutions and of water at various pressures was 
thus obtained. The values of ¢,, the boiling point of 
a solution of dextrose, and ¢,, the corresponding 
boiling point of water at the various pressures were 
related by equations of the form 


a+bt,+e @, (1) 


in which the constants a, 6, and ¢ for each concentra- 
tion of dextrose were evaluated by the method of 
least squares. The values of ¢, and the correspond- 
ing pressures taken as reference standards were 
selected from the compilation of Osborne, Stimson 
and Ginnings {1} * and have previously been tabulated 
in convenient form [2] 

No equation was found to relate the boiling points 
of the solutions to their concentrations, at a given 
pressure, with the precision of the experimental 
data. For this reason the boiling points of solutions 

' Research Associate at the National Bureau of Standards, representing the 


Corn Industries Research Foundation 
? Figures in brackets indicate the literature references at the end of this paper 
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of even values of concentration were obtained graphi 
cally by plotting At, the boiling point elevation given 
by the difference between ¢, and ¢,, with respect to 
the concentration (or At/concentration versus con- 
centration). From the resultant series of curves 
one for each of the standard reference pressures, the 
values of At for even concentrations of dextrose were 
obtained, and the boiling points of the solutions were 
evaluated by adding the corresponding boiling point 
of water. These boiling points of the solutions were 
tabulated with their corresponding pressures and 


the constants A, B, and C in the Antoine equation {3} 


B 
log p A= ? 2 


were evaluated to obtain the relationship between 
vapor pressure and temperature for each of the 
several even concentrations. (The symbol “log’’ is 
used in this paper to denote the logarithm to the 
base 10.) Equation 2 is explicit in temperature 
when written in the form t= B/(A—log p)—C. Also 


dp P\ A log py 
dt Bloge — 


The best equation that was found to represent th 
relationship between the measured boiling point 
elevations and the concentrations is of the form |4 

log At=a+Be+ye*+ de’, t 
in which the constants a, 8, y, and 6 were evaluated 
by the method of least squares. But a comparison 
of the computed with the observed values of tl 
boiling point elevation at a given pressure shows 
deviations that are several times larger than 
precision of the original data. 
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II. Apparatus and Materials 


xcept for minor changes, the apparatus was the 
as described in a previous paper [2]. Two 
e barometric ebulliometers, of the type devel- 
oped by Swietoslawski [5], were sealed to a manostat. 
I, the ease of the solution ebulliometer, a flask 

taining water was interposed to inhibit the loss 
of water vapor from the solution to the manostat. 
Pressures below atmospheric were obtained with a 
mechanical pump, and those above were obtained 
1y the addition of nitrogen from a cylinder. 

lhe solutions of dextrose were prepared from NBS 

Standard Sample dextrose previously dried in a 
vacuum oven at 60° to 70° C for 6 hours and weighed 
into known weights of water. All concentrations, 
expressed in weight percentage on a vacuum basis, 
are corrected for the vapor and liquid hold-up that 
prevailed in the ebulliometer during the boiling 
procedure. Crystalline a-dextrose unde ‘rgoes mutaro- 
tation in aqueous solution to form an equilibrium 
mixture of a- and 8-dextrose [6]. The rate of the 
conversion is dependent on the temperature and 
pH of the solution, and the equilibrium amounts of 
a- and 8-dextrose present depend on the temperature 
and concentration of the solution. Hence the 
solutions studied contained a mixture of a- and £- 
dextrose. However, since both forms possess identical 
molecular weights, no effect on the boiling point 
or the vapor pressure of a solution is to be expected 
as a result of different ratios of the two isomeric 
forms at different temperatures and concentrations. 
The solution ebulliometer was charged with a 
known weight and volume of solution, the water 
ebulliometer with the proper volume of distilled water 
and the filling tubes sealed off. Boiling tempera- 
tures were measured at as near the same boiling 
rate as could conveniently be achieved, the rate being 
controlled to return approximately 20 drops of 
condensate per minute to the boiler as observed in 
the drop-counter placed in the return tube. This 
boiling rate approximates the midpoint of the flat 
constant-temperature portion of the curve obtained 
by plotting boiling temperature with respect to 
boiling rate [5] and was determined by initial experi- 
ments. The comparative measurements of tem- 
perature were made with a platinum resistance 
thermometer (coiled filament type) and a Mueller 
thermometer bridge by the method described in 
another paper from this laboratory [7]. The tempera- 
tures were measured to 0.001 deg C, with an average 
reproducibility of 0.002 to 0.003 deg C and an 
estimated accuracy of 0.005 deg C on the Interna- 
tional Temperature Scale. It was not necessary to 
hold the bridge at constant temperature, as the 

inperature corrections to the resistances were 
practically identical for the temperatures of the ref- 
erence and measured substances in these comparative 
Measurements. 

‘he ratios of the boiling points of the solutions 
hose of water at normal atmospheric pressure 
nat the beginning and end of a series of measure- 
ts did not differ by more than 0.008 percent in 


the extreme case of the most concentrated solution. 
The average difference for the several concentrations 
studied was 0.003 percent. Refractive index measure- 
ments on the original and boiled solutions were 
practically identical. It is believed that no sig- 
nificant change in concentration occurred during the 
boiling procedure, although the solutions of higher 
concentration assumed a slight yellowish tint, par- 
ticularly after boiling at the higher pressures. Spe- 
cific rotations were not measured. 


Ill. Experimental Results 


The correction for the water vapor and liquid 
hold-up in the ebulliometer, to be applied to the 
concentrations of the solutions as originally pre- 
pared, was determined experimentally by using 
solutions of sodium chloride. For this purpose, a 
stopcock was sealed temporarily at the bottom of the 
solution ebulliometer for rapid withdrawal of samples 
while boiling. Solutions of sodium chloride of dif- 
ferent known concentrations were placed in the 
solution ebulliometer and boiled in the usual fashion 
at normal atmospheric pressure. When boiling 
equilibrium had become established, as evidenced by 
steady boiling temperatures, a 25-ml sample was 
quickly withdrawn into a flask that was surrounded 
by ice water to minimize the loss of water vapor. 


TABLE 1. Determination of ebulliometer hold-up 
Original Boiling 
_— concentra- Boiling concentra- id. 
harge tion of rate tion of Hold-up 
NaCl NaCl 
Percent Drops/min Percent g 
65. O86 11. 760 ? 12 11. 873 0. 62 
65.119 11. 760 20 11. 891 71 
65. 135 11.760 0 11. SYS 7 
70. 008 22. 195 12 22. 391 6l 
69. 963 22. 195 20 22. 413 (is 
70. O87 22. 195 30 22. 440 76 


The samples were analyzed for sodium chloride by 
evaporation to dryness in platinum dishes on a 
steain bath, followed by heating to incipient fusion 
and weighing the sodium-chloride residue. The dif- 
ference in concentration between the original and 
the boiling solution, together with the knowledge of 
the amount of original charge in the ebulliometer, 
provided data for the calculation of the hold-up of 
liquid and vapor while boiling was taking place at 
normal atmospheric pressure. It was assumed that 
the liquid hold-up remains essentially constant at 
all pressures involved. The vapor hold-up, which 
depends on the pressure and temperature, could be 
closely estimated from the knowledge of the free 
volume of the ebulliometer. Since this method of 
determining the hold-up may be of interest, the 
results are given in table 1. From these results, the 
amount of liquid and vapor hold-up at 20 drops a 
minute, which was the boiling rate employed in this 
work, was taken as 0.68 g, of which 0.61 g is at- 
tributed to liquid hold-up on the walls of the ap- 
paratus, and 0.07 g is a calculated value for the 
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Pasie 2. Corresponding boiling points of dextrose solutions and water 
il 
e 
10.44- 21 2f- 2.50 40.31- 50.4 59.6. 3.23- 
percent Water percent Water percent Water percent Water percent Water percent Water percent W in 
solution solution solution solution solution solution solution e 
re 
‘ ( ( r ( ( c ( ( ( ( “< c 
100. O44 0. 724 100. SH7 08. TRS 101. 286 9. 889 101. 686 99. 686 108. 237 100. 068 104. 806 100. 059 105, 386 a cn 
100. 083 98. 732 100. 568 99. 783 101. 263 09. 876 101. 685 09. 684 1. 246 100. O72 104.813 100. 066 101. 148 us ve 
4. 182 oS. 868 04. 207 @. 514 05. 4533 4. 163 95, 258 &. 330 99. 271 6. 170 GS. 302 3. 711 96. 579 ] 
oF. a a RAI 04. 273 a3. 520 87. 802 86. 490 87. 725 RS 879 108. 323 100. 140 oO. SSO 86. 462 91. 387 Sf) te 
a6. 414 86. 109 86. G27 86. 207 81. 624 80. 358 80). 308 78. 523 99. 144 96. 044 BO. 696 76. 489 MS44 79 ’ 
i ¢ 
aol 433 Nf 128 al 40 al 2) v2. 84 91. 513 93. 033 91.122 O4. 684 91. 673 69. O77 65. 092 76. 686 71 } ist 
75. O46 75. #52 6. 675 75. 2 101. 389 09. O65 101. 483 99.472 80. 448 86. 522 oO. 914 86. 488 105, 298 Wy 
75. O77 75. OM 76. 704 76. 021 101. 653 100. 229 101. 701 09. 6R6 82. 844 80. 011 104. 634 99. 882 110. 689 104 rive! 
100. 261 09. 929 4. 7F 4. 22 101. 653 100, 233 101. 608 U9. 68S 75. 432 72. 705 111. 205 106, 362 115. 148 109 +! 
100. i 09. G28 . 940 4. 202 106. 877 105. 417 106. 970 104. 808 1038. 161 Y9. ORS 117. 229 112. 147 119. 269 11 - 
105. 559 105, 218 78. 142 77. 455 111. 307 109. 79s 111. 431 109. 306 108. 162 99. 990 122. 282 117. 38 122. 843 116. TRS 
105, 562 105, 220 87. 3S aH, HOR 115. 20 113. 683 15. 416 113. 248 108, 270 105. O11 125. 670 12a). 369 126. 482 120. 304 LABI 
109. 860 100. 507 4. 550 %. 793 118, 782 117. 212 19. 115 116. 895 112.778 109. 442 116. 143 111.123 105, 324 ws 
109. 8A2 109. 511 1). 616 W. 832 122. 419 120. 805 122. S2e 12). 24 116. 880 113. 457 104. M2 8. 828 
113. 867 113. 508 100. 791 100. 007 101. 522 100. OGS8 116. 730 114. 42 120. 461 116. 972 
113. 873 113. 512 105. 994 105, 1Sé 110. 091 107. O87 123. 804 1D). 324 
117. 356 116. GSS 110. 298 109. 463 101. 731 99. 722 103. 191 100. 081 
117. 34 116. O04 114. 495 113. 649 ‘ 
12). S38 12). 461 117. 897 117. 030 
1a). S44 1D). 462 121. O07 12). 118 
118. 216 117. 84 
118. 224 117. 8%) 117. 891 117. 023 } 
114. 205 113. 843 114. 405 l i) 
114. 214 113. 851 110. 208 109. 469 
110. 4%) 110. 0@ 105. 852 105, 044 
110. 458 110. 104 100. 689 moO 
106. 190 105. S44 
106, 20 105, B55 
100. 2 w aM 
100. 230 WO. 806 
100, 2 99. 902 
4. 021 ut fay press 
4. 04 Ws. 712 ° 
86. 212 85. G05 devia 
SH. 22S a5. O22 
at th 
76 Os 75 ow \ alu: 
76. 124 75. 825 
62. 585 62. 207 press 
62. 641 62. 344 — 
77. 031 76, 732 press 
Ch 
Si), JSS 45. GRD 
4. 04 O83. 76% and ( 
100. O80 wy. 747 ; 
able 
and | 


vapor hold-up at 760-mm pressure and the boiling 
temperature. At 200- and 1,500-mm pressure the 
calculated vapor hold-up is 0.02 and 0.12 g, re- 
spectively, giving a total hold-up at these extreme 
pressures of 0.63 and 0.71 g, respectively, assuming 
that the liquid hold-up remains constant over this 
pressure range. 

The comparative boiling points of dextrose solu- 
tions of several concentrations and water, in the 
order as actually measured, are given in table 2. 
The concentrations listed are those prevailing during 
the boiling procedure at 760-mm_ pressure. The 
change in concentration accompanying the change 
in vapor hold-up with variation of the pressure is 
small but can be significant, particularly for the 
more concentrated solutions. For example, a 10- 
percent solution changes in concentration by —0.01 
percent at 200 mm and by +-0.01 percent of dextrose 
at 1,500 mm. A 60-percent solution changes by 

0.04 and +-0.04 percent at these extreme pressures. 

By means of eq | the boiling points of the dextrose 
solutions were expressed as functions of the boiling 
points of water over the pressure range involved. 
The constants a, 6, and ¢ for each of the solutions are 


given in table 3. The average and maximum devia- 
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sures 


tions of the experimental points from those caleu ate 
lated by eq 1 with the appropriate constants ar at 
given in the last two columns. The average devia- 

tions increase from +0.002 deg C in the case of the 

more dilute solutions to +0.013 deg C for the most 
concentrated solution. An explanation for the i 

creasing experimental deviations at higher conce! 

trations of dextrose is that the boiling equilibrium 

in the latter cases is established with less certaint) 

The pumping action in the ebulliometers is probabl) 

less efficient because of the increased Viscosity of th 

solutions and, also, the concentration of the solutior - 
in the boiler may vary due to incomplete mixing of . 
the condensate with the highly concentrated so son 
tion. as 

The boiling temperatures of the solutions at 

standard reference pressures were calculated by tay 
means of eq 1 with the constants from table 3 


0 


difference between the boiling temperatures 0! 100 
solution and water at a given pressure represents 10 
boiling-point elevation, Af, at that pressure. | - 
each reference pressure, At was plotted with resp wn 
to the concentration on a scale that permitted 504 

16 


plotting and reading of temperatures to + 0.005 di 
and of concentrations to +0.01 percent of dextr 














ilues of Af at even concentrations were read from TABLE 4 Rowling points of dextrose solutions at standard 
. ; ~ referer ress 8 
e curves and added to the corresponding boiling en 
ints of water to obtain the boiling temperatures of 
e solutions listed in table 4 at the reference pres- Boil Boiling point of solutions 
¥ ing 
res. A separate set of the constants A, B, and ¢ Pressure | point 
eq 2 was evaluated to fit the data for each solution a 10 20 10" 40 50" fn) 
ven in table 4. Values of the constant C were 
termined by the method of averages, followed by —_— ( ( ( ( ( ( y 
evaluation of constants A and B by the method of ao) 2) Se Sa Sr Se SS oS 
ist squares. The values of the constants are 230.13 | 75 | 75.286) 75.631) 76.101 | 76.718 | 77.607 79. 235 
~ rr 2 . 355. 22 sO 80). 280 SO. O48 81.130 81. 767 82.772 84. 338 
ven in table 5 lhe over-all average deviation of ats | es | enaes| eness | ante: | oneee| teee| oe tee 
e vVvapor-pressure-temperature relationships ex ee Pe ee ee ee ees eres oe 
633. 99 YS 95 10 4 Ff 6, 220 06. 926 VS. O11 09. 680 
‘ > j ” 7H0. 00 100 100. 319 100 OL. 265 101. G&S 103. 006 104. 804 
i BLE 3 Roi ing points ol dertrose solutions in terms of the On6, 06 105 10 29 105 7 106 303 107.042 108. 183 109. 93 
bowling points of water 1, 074. 58 110 110.340 | 110.77 111.343 112.102 113.273 | 115.070 
1, 268.0) 115 115. 352 115. 793 116. 385 117.163 118. 365 1D). 212 
1+bf «+ct 1, 489.14 120 120. 36. 120. 819 121. 428 122. 225 123. 460 12 ‘ 
Concet M 
i n a 1 f ‘ A verage man , . 7 
oak deviation ian were evaluated by the method of least squares 
lhese constants are given in table 7 for several refer- 
Percent ( ( ence pressures over the range 200 to 1,500 mm. As 
). 34 0. 337% 0. 997929 0. 2011K10 0. 002 +0. 008 2 . . . a . 
—¥- ++ othe 3377 —— po the deviations given in the last two columns are con- 
-+~ sr abe 000 = +. 00 siderably larger than the experimental precision, 
10.3 N27 1.00 sH2 OO Oo . a. poe 
0.54 | 1.9000 | 1.007553 511 007 o12 these equations were not used to evaluate At. They 
1 69 3. 2872 1. OOSH18 1.0054 009 +. 026 
63.2 4.7004 ISKGAS ) O42 013 23 
TABLe 5 | apor-pressure—tem pe rature relationsh ps of 
dextrose solutions 
pressed by eq 2 is +0.04 mm. The larger maximum 
deviations shown in the last column of table 5 were ___B 
at the highest pressure, 1489.14 mm, in every case. ( 
Values of the temperature and rates of change of 
: . . Cones Ave May 
pressure with temperature at even values of the Rae ay , B ; oH <r 
pressure are given in table 6. lextrose ati feviation 
The relationship between boiling point elevation 
: get rereen nm Hg mm Hg 
and concentration was obtained by using the data in — . Tr one on +- ~ 
table 3 to calculate the boiling points of the solutions m0) 7 1633.713 | 224. 618 0 4 
’ “ue . Ww) r it7. 2Y 22 ay o 2 
and boiling point elevations at the reference pres- 10 7 1636. 304 224. 764 ol 02 
\ : . : w ? 1622. 430 225. 20) Oo 11 
sures. From the tabulations of At and concentration ro canr agn | Sim 980 oe 17 
at a given pressure the constants a, 8, y, and 6 in eq 4 
TABLE 6 Values of pressure, temperature, and rates of change of pressure vith te mperature for dextrose solutions 
Concentration of dextro s it 
' Peres 1) Percent 30 Percent 40 Perce Percent 60 Perce 
rempera , rempera ipidt rempera ip/at Ie i ipidt Pempera ip/dt Pempera ip/idt 
ure 7 ture ah ture " ure - ture ture 
He ( ' ( c nm/°¢ ( nm/°¢ ‘ nm/°¢ ( nm/°¢ ( n/?¢ 
fl 1.87 4.02 217 4.91 52. 5S 100 3. 11 1. 8A Lo 4. Nt ». 42 is 
(6. 72 SN 67. O8 &S4 67.50 & 82 68. ON x79 69. O02 8.7 70. 51 s 
4 76.17 12. 47 7 2 12. 45 76. 09 12.41 77. 61 12. 3 7s ‘ 12. 31 sO. 14 2. 24 
83.2 87 i2 15. 84 84.11 15.79 84. 7E 15.7 88.77 15.¢ 87. 3¢ , 
RS. ON ) SU. Me 19. O% RU. AT 9.03 Cia | 18. 9 91. SS 15.8 3. 21 1S 
Ls a ed 22. 04. 21 22. 20 04.73 22. 14 a5 42 22.0 uw wh 271.9% VA. If 21 
” ~~ OO 2 “0 YS. 438 25. 2 Gs. 06 5 16 8. 67 25. 06 100. 77 24. ¥2 102. 46 24. 72 
‘ 100. 32 27.08 100. 73 27.01 101. 27 26. 93 101. OS 2%. 82 03. 10 MH. OF 104.81 2H. 4 
800 101. 77 28.2 102.18 on. 18 102. 72 28.10 3.44 7. OS 104. 5 27. 82 106. 29 27. SY 
105. 14 1.14 105. 55 1. 06 106. 11 0. 97 106. S85 10. 84 107. 99 6 108. 74 wm. 40 
TLD 108. 21 $3. OF 108. 64 Sy 109. 20 3.77 109. 5 3. 6 111.11 3. 43 112. 89 bb. 14 
10 111.04 6. 73 112. 47 6. 63 112.0 iti 112. 81 6. 37 113.990 wi. 14 15.79 52 
a0) 113. 67 4 114.11 0 34 114.69 19, 2. 115. 46 19. 06 116. 65 tS] 118. 48 ms. af 
100 116. 12 42.12 116. 57 42. 01 117.15 41. 89 17 41.71 119.1 41.44 121. 00 41.0 
100 118. 42 4. 7t 118. 87 44.04 119. 47 4.50 120. 26 44. 41 121. 49 4.02 iz it Sat 
0) 120. 59 ‘7.3 121. 05 $7. 22 121. 65 47.08 122. 4 46. 88 123. 70 46. 57 125. 59 46.11 
600 122.¢ 49.91 123.11 49. 78 123. 72 49.6 124. 4 49. 41 125. 79 49.08 127.70 is. 50 


Extrapolated value 
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TABLE 7 Boiling-point elevation of dextrose solutions 
log At=a+A¢ Cc*+8c% 
Pr u . g 4 Average maximum 
deviation deviation 
mm lig c ( 
233. 72 ), GRUR4 5. 113375X<10-2 7. 377463 10 5. 526503 X 10 0.017 0. 032 
‘55. 22 e274 }. 191070 7. 435453 5, 463234 +. 023 + O45 
25, st TSTI 5. 239106 7. 477459 5. 427242 +. O28 + 0 
60. OO 44020 5. 235447 7. 42722 5. 345087 +. O31 +. O61 
074, 5S 1824 251193 7. 400808 SOHO +. 033 +, OO7 
450.14 RASA 45 225459 7. 470027 5. 404016 + O33 + (7 


are useful, however, in calculating the boiling point 
elevations of dextrose solutions to within 0.05° to 
01°C 

Molal boiling point elevations of aqueous solutions 
of dextrose have been measured by Juettner [8], but 


only at normal atmospheric pressure. A com- 
PARLE Vo ho -point elevation of dextros ) or 
60 mm H 
Mi J 5 | 


parison of his values with those of this investig 
is given in table 8. Apparently, there are no 


published measurements of the vapor pressures 


aqueous solutions of dextrose in the range re pe 
in this paper 
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Wavelengths for Calibration of Prism Spectrometers 
By Earle K. Plyer and C. Wilbur Peters 


and other comy ounds 


Several absorption bands of poly stvrene, 1,2.4-trichiorobenzene 
ave been measured in the infrared region from 1.5 to 24 uw on grating spectrometers. These 
bands have been determined with sufficient accuracy for use in calibration of prism instru 
ments \ table is included that gives the cell thicknesses used in the measurements. The 
emission lines of mercury for the region from 0.5 to 2.4 yw are included In order that the 
table may be of most value, a number of bands that have been determined by previous ob- 
servers have also been included 


Many infrared absorption bands have been care- 
fully measured, and their reported wavelengths [1] * 
are useful for calibrating spectrometers. Techniques 
of calibration utilizing vibrational bands with re- 
solved rotational structures, including those of am- 
monia, carbon dioxide, and water vapor, have been 
deseribed by Oetjen, Kao, and Randall [2]. When 
one undertakes to calibrate an infrared spectrometer 
he soon finds that additional reference lines or bands 
would be extremely helpful and that the information 


available has some serious disadvantages. In cer- 


tain regions there are an insufficient number of 
University of Michieas 
* Figures in brackets indicate the literature references at the end of tt paper 


standards. In some cases where a band has a 
tional fine structure that is not resolved by a p 
it is not possible to locate any individual band 
sufficient precision to justify its use for calibra 
Furthermore, some absorbing materials requir 
practically long absorption cells. To avoid thes 
advantages and increase the number of calib: 
points, additional lines and bands in the region 
the visible to 24 «» have been measured by 

sources or absorbers convenient for the calibt 
of prism instruments. These include: (1) pol 
rene films and trichlorobenzene, which provide st 
ards between 15 and 24 uw; (2) AH—4 mercury 

in the region visible to 24 uw; (3) toluene at 2 
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20, 


at 


y- 











1) methyleyclohexane [3], 7 to 12 wu. Wavelength 
lata on all except the mercury source were obtained 
n order to meet the objectives of this investigation. 

This information, combined with previously re- 
ported values for carbon dioxide, methane, water 
vapor [1], and methanol [4], form a reasonably com- 
ilete set of secondary calibration points, supplement- 
ng the bands of ammonia and water vapor now 
n use, 

All the measurements recorded in table 1 and 
marked on the characteristic spectra that follow 
were made with grating instruments. Those in the 
region 1 to 3 w were made at the National Bureau of 
Standards with a 15,000-line-per-inch grating and 
ead sulfide detector. The polystyrene bands in the 
}4-u region were measured for the authors on the 
Ohio State University grating spectrometer by John 
5. Burgess. The bands from 3.6 to 24 uw were 
measured by the authors at the University of Mich- 
igan. These grating measurements are the basis for 
he values marked on the characteristic spectra 
shown in figures 2 to 8. The characteristic spectra 
as shown in figures 1 to 6 were taken with an appro- 
priate prism instrument. 

In making the grating measurements, the resolu- 
tion was adjusted to be of the same order as that 
obtained from good prism spectrometers for the 

gion appropriate to a given optical material 
Fairly wide slits were therefore required. The band 
positions were determined with a constant-energy 
background continuum to avoid the confusing shift 
that occurs when the background is not constant, 
that is if the spectral energy record shows a steep 
slope, the apparent positions of absorption minima 
may be shifted from the correct wavelength. 

The bands shown in figures 2 and 4 have been 
measured with high resolution, so that an accuracy 
to +2A could be obtamed. For the calibration of 
nstruments using lithium fluoride prisms, this 
legree of accuracy is desirable. 

In the spectral region from 4 to 14 yg, the grating 
measurements have been made to 0.001 yu of the 
wavelength value, and from 14 to 24 uv to 0.01 4 
The following paragraphs give a more detailed 
account of the procedures in handling the experi- 
mental materials and an explanation of specific 
features of the graphic and tabular data. 


Figure 1 shows the mercury lines, and the wave- 
lengths recorded in table 1 have been supplied by 
Curtis J. Humphreys of this Bureau from his meas- 
irements and those of other investigators |5]. These 
spectra were recorded with a Perkin-Elmer 12—A 
spectrometer with a lithium fluoride prism. The 
nes were measured with different slit widths, but 
he drum readings of the maximum intensity of a 
ne did not vary appreciably with the slit width 

hen the lines were not complex. In the table, 

ies that were not resolved are joined by brackets, 
ul the wavelengths given on the curve are the 
verage values for these groups of lines. 

The polystyrene films are made by first preparing 

10-percent solution of polystyrene in toluene and 
ourmg some of this solution on a smooth, clean 








TABLE 1 Calibration wavelengths for prism spectrometers 
Wave 
Wave- | number State Description Substance 
length vac 
cm 
0. 54607 | 18307.8 | Emission AH lamp Mercury 
576061) 17327 do do Do 
7007 || 17264. 4 do do lo 
1, 01398 USSU. 4 do do Do 
1. 12866 SS857.7 do do Do 
1, 35703) 7367.0 lo do Do 
1. 36,28) 7311.0 do do Do 
1. 30506 7166.3 lo do ) 
1. 52042 6536, 2 do do 0 
1. 6606 6020. 3 Liquid 0.5 mm cell 1.2.4-Trichlorober 
7m 
1. 69202) G08. 5 Emission AH—4 lamp Mercury 
HO4101 5000. 4 lo do Do 
1.70727\ 5855.7 | do 1k 
1. 7100N) 843.0 1 lo Ik 
1. 81307 14.0 i do D> 
97009 74 j 1 > 
2. 1526 1644.2 Liquid 0.5-mm cell 1,2.4-1 ! t 
2. 24929 4444. ¢ I ! AH lamp Mercury 
12t $322. 9 Liquid ( mm cell 1,2,4-'1 t 
254 1200. 1 I AH-4 | Mercur 
4030 4160 Liquid cell 1,2,4-1 t 
2.43574 f101.¢ i 1 1) 
2. 4044 107.9 1 i D 
2. 5434 U3 d " D 
2.7144 683. 0 Vapn { cell Methanol 
224 sim 4 S 2h-p f Polys 
a4 10ND. 6 i , is 
Qt 00) j j I 
i wrt 1 hk 
WY wes i i Lo 
it 20 G m cel Met! 
4188 2u24 Solid 25-y fil Pol y 
7S Pt ] j I) 
4. 255 4 (ra \ t ( th 
1. Sit 2054 Val Mett 
138 1V45.8 Sol ips f I ysty 
4 Sol 1) 
6, 238 142 i I> 
6. 69 1404 I) 
iN I i Oo. Mett 
1 672 1034.2 Vay cn M i 
724 102s Solid p fil I tyr 
1 Ty) i > 
1.4 | Liq Th Nett hex 
s M428 j j I) 
aA ?. | ‘ Carb liox 
4. 20 Aug. ¢ Solid Poly 
4. OS tu G Cart x 
is O45. 8 Liquid { ics Unknow t 
eal grad + 
hlorober 
‘ ! 2 1,.2.4-1 } } 
in] “) i 1h 
2. st i of m cell 2 I } } 
Lins it 
21. 52 Hi Th ell Toluer 
21. si SS. 12 ll 2,4-1 ! 
22. 7t 439 i > 
23.8 419.2 Vapor A pl Wa 


“Broad bands 


piece of fairly fine plate glass, which is resting on a 
level surface. A doctor blade (number 15 for a film 
of 25 w in thickness) is drawn over the solution to 
make the film more nearly uniform in thickness. 
The film is then allowed to dry at room temperature; 
the surface should be protected from dust and other 
particles while drying. To remove the dried film 
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from the glass plate water is applied, and the film 
readily floats off the glass. After the film is redried, 
it is ready for use. The thickness may be checked 
with either a Swedish gauge or a micrometer. In 
figure 2 are shown the bands of polystyrene in the 
region of 3.4 » as recorded with the lithium fluoride 
prism. Seven bands are observed in this region, and 
six have been selected for use in calibration. In 
figure 3 is shown the spectrum of a film of polystyrene 
from 2 to 16 4. The grating instrument wavelength 
values of the bands selected for calibration have been 
marked on the curve, which was a record obtained 
from a Baird infrared spectrophotometer. 

The trichlorobenzene samples, both highly purified 
and technical grades, were obtained from Hooker 
Electrochemical Co., Niagara Falls, N. Y. The 
authors are indebted to R. H. Kimball of that com- 
panv for furnishing highly purified samples of both 
1,2,.3- and 1,2,4-trichlorobenzene. In figure 4 the 
spectrum from 1.5 to 2.5 u of highly purified 1,2,4- 
trichlorobenzene as obtained from the Perkin-Elmer 
12—-A spectrometer with a lithium fluoride prism is 
shown. The seven bands that have been selected 
for the calibration have the wavelengths marked on 
the curve. These bands were not very broad even 
when measured on the grating instrument, and the 
region of maximum absorption could be determined 
to the fourth decimal place of wavelength in microns. 
Original grating instrument recorder traces of two 
bands of 1,2,4-trichlorobenzene in the 2-y region are 
given in figures 7 and 8. The wavelength of the 
band at 1.6606 uw was determined from the mercury 
emission lines superimposed on the spectrum. To 
achieve this a Cetron mercury lamp, similar to the 
S-—1 sun lamp, was used as the source. By focussing 
one of the hot electrodes on the entrance slit, the 
absorption and emission spectra could be measured 
simultaneously. In the determination of a second 
band in this region, figure 8, water vapor lines that 
appear on both sides of the band were used as stand- 
ards. On these two figures the wavelengths of the 
reference lines have been indicated as well as the, 
wavelengths of the two bands of 1,2,4-trichloroben- 
zene 

The spectrum of technical grade 1,2,4-trichloro- 
benzene from 14 to 24 uw as recorded by the Perkin- 
Elmer 12-A spectrometer is shown with the wave- 
lengths noted in figure 5. The technical grade con- 
tains some 1,2,3-trichlorobenzene, and some of the 
bands in the spectrum are due to the 1,2,3-isomer 
In order to distinguish the bands that arise from the 
1,2,4-trichlorobenzene and those from the 1,2,3- 
trichlorobenzene, the spectra of highly purified 
samples of each of these isomers are shown in figure 
6. The band at 20.56 u is seen to be due to 1,2,3- 
trichlorobenzene. However, the band at 15.48 xu 
does not appear in the spectra of either of the purified 
samples and thus probably arises from another 
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chlorobenzene contained in the technical 
1,2,4-trichlorobenzene. Although the results 
not shown here, a mixture of about 85 percen 
1,2,4-trichlorobenzene and about 15 percent 
1,2,3-trichlorobenzene was also measured from 


‘tai 
=f 


to 24 uw. No changes in the positions of the ban |s 
that were selected for calibration points were o}- 
served in the prism measurements when the spectr 
of the mixture was compared with that of the t 


nical grade 1,2,4-trichlorobenzene. The value of 


22.76 uw is a prism instrument determination and js 
probably accurate to +0.02 uw. As the path from 
the source to the detector for the grating instrument 
is greater than that for the prism instrument, thy 
atmospheric water vapor absorption is also greate 
This absorption in the region of the 22.76-u band 
makes it difficult to determine the position of th: 
band, and thus the prism spectrometer value is 
given. The bands that are broad and not as ac- 
curately determined are listed in the table and ar 
denoted as broad bands. 

In table 1 is given a list of bands and emission lines 
that extend from 0.5 to 24 uw. All data not obtained 
in the current set of experiments are quoted from th 
sources listed. The absorption band of water vapo: 
at 23.85 uw is made up of several lines, and the listed 
value was obtained with low resolution of the grating 
spectrometer. The wavelengths are given in microns 
and the wave numbers in vacuum in em™!. The 
reduction of the values of wave number to vacuum 
was accomplished by using 1.0002706 as the index of 
refraction of air. The index of refraction of air is not 
a constant, but the variation in the value is small and 
would introduce inappreciable changes in the wav 
number. The list includes bands measured in this 
work and other bands, which should give sufficient 
points for calibration of prism instruments. Addi 
tional wavelengths beyond those shown in the tabl 
are needed in some regions, and it is hoped that 
further measurements will be made in the futur: 
This work has not included any measurements by 
yond 25 yw, but there are a large number of wate: 
vapor lines that are very well suited for calibration 
in this region [6, 7]. 
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rhe inserts were obtained by diluting the trichlorobenzene 1:4 in carbon disulfide, 
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Infrared spectra from 14 to 24 pw of highly purified samples of 1,2,4-trichlorobenzene and 1,2,3-trichloroben zene. 
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High-Pressure Apparatus for Compressibility Studies and pro} 
. ° ° evlu 
its Application to Measurements on Leather and Collagen hol 
of t 
By Charles E. Weir te: 
con 
The design and construction of apparatus to be used to measure volume changes of vess 
solids (or liquids) between 1,000- and 10,000-atmosphere pressure is described in detail plug 
Calibration of the equipment and its use in determining the compression (— AV/ Vo) of leather the 
are discussed. The compression of leather between 1,000 and 10,000 atmospheres is reported. rare 
The compression between 2,000- and 10,000-atmosphere pressure is approximately 7 percent a ; 
and appears to be relatively unaffected by moisture content, type of tannage, or sample — 
variation. The compression of all leathers tested is described by the equation A/V thre 
1.23 « 1075 (P— 2000) — 5.60 10- (P— 2000)? + 1.35 10-" (P— 2000 in 
nigh 
I. Introduction performed in few laboratories. To obtain informa- dian 
tion on the apparatus and techniques involved i with 
A Bureau program of measurement of physical such measurements, preliminary experiments wer mae 
constants of leather and collagen has resulted in a conducted on equipment that has been used for hard 
determination of the expansivity of leather and many years in the Geophysical Laboratory of th the | 
collagen [1].'_ These experiments also demonstrated Carnegie Institution of Washington, for measur bore 
that the shrinkage of leather, heretofore considered ments at pressures as high as 12,000 bars. Sub- prod 
as analogous to a melting, was in reality a phenome- sequently similar equipment has been constructed , 
non occurring over a range of temperatures and at this Bureau and used in conjunction with a hydrau- ieee 
resulted in an increase in real volume of the leather- lic press in the Geophysical Laboratory. ‘This shies 
water system studied. Subsequent studies on the report deals with a description of the apparatus and the | 
rate of shrinkage as a function of temperature, a part of the data obtained. ore 
tannage, liquid medium, ete., have been reported ppl 
9° ; : : -e : ped 
ad Vie " » S » r . r 
[2.3]. In view of the results an investigation of the II. Apparatus "he 
effect of pressure on leather and the shrinkage th 
process was indicated, including a determination of In theory the apparatus required is simpl mae 
the compressibility of leather for which only esti- practice, however, limitations imposed by the strength . 


mates have been available [4]. 
Experiments involving appreciable pressures are 
extremely exacting mechanically and have been 


Figures in brackets indicate the literature references at the end of this paper 


of materials require exacting machine work and t! 
best of materials. An experiment consists essent\:!l) 
of forcing a piston into a vessel containing the samp! 
under study immersed in a suitable liquid, and re 
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» the depth of penetration of the piston at known 
fernal pressures. 

[he apparatus required consists of 

|. A thick-walled pressure vessel having a smooth 
ve and a second opening provided to permit 
asurement of the internal hydrostatic pressure. 

». A leak-proof piston, or plug, which is forced into 

e bore. 

; Aram designed to drive the plug into the vessel. 

t+, A leak-proof plug carrying an_ electrically 
nsulated lead from the internal pressure-gage to the 

ternal pressure-measuring equipment. 

5. The pressure-measuring equipment. 

6. A hydraulic press to actuate the ram. 

7. Assorted extractors to remove washers, etc. 

The pressure vessel and accessories to be described 
are made in accordance with designs perfected by 
workers at the Geophysical Laboratory and are the 
result of wholehearted cooperation by these workers. 
The equipment is very similar to that described in 
considerable detail in 1919 by Adams, Williamson, 
and Johnston [5], but it is believed that a more 
recent highly detailed description of the apparatus, 
its fabrication, and use may be of interest and value 
to potential workers in this field. 

The pressure vessel consists of a cylinder of SAE 
1340 steel—a carbon, nickel, chrome, manganese, 
molybdenum, alloved steel having ‘‘deep-hardening”’ 
properties—5 in. in diameter and 12 in. long. The 
evlinder, partly shown in figure 1, contains an axial 
hole % in. in diameter and 7 in. long. At the bottom 
of this bore a 2-in.-long hole is drilled, which de- 
creases from ‘, to is in. in diameter and serves to 
connect the pressure gage with the contents of the 
vessel. An opening for the pressure gage and its 
plug is made into the side of the vessel to connect with 
the bottom of the \¢-in. diameter axial hole. The 
gage opening consists of an outer hole 1% in. deep 
designed to fit a bolt 1% in. in diameter having 12 
threads per inch; a central portion % in. in diameter, 
lin deep having a flat bottom that meets the wall at 
right angles; and an innermost portion %. in. in 
diameter extending slightly beyond and connecting 
with the axial '\s-in. diameter hole. After initial 
machining, this vessel was heat-treated to procure a 
hardness of Rockwell C-40. The scale produced in 
the heat treatment was removed, and the %-in. axial 
bore that was distorted was ground and lapped to 
produce a fine finish 

Before the pressure vessel could be used for high- 
pressure work it had to be subjected to a pressure 
seasoning followed by a refinishing operation. After 
the initial machining, the bore of this vessel was ob- 
served to increase in diameter by several mils on 
application of pressure of approximately 10,000 atm. 
The vessel, however, can be rendered stable, at least 
n this pressure range, by a seasoning process consist- 
ing of subjecting the vessel to an internal pressure far 
in excess of any subsequently contemplated pressure 
For this seasoning, special plugs, which will be 

cribed later, were used to fill the openings in the 
vessel. In seasoning, pressure was applied in incre- 
ments of 2,000 atm, pressure being maintained for 30 
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FIGURE 1 Section of pressure vessel 




























min after each increment, until a maximum calculated 
internal pressure of 23,000 atm was attained. At the 
conclusion of this operation all fittings were found to 
be damaged, while the vessel had increased noticeably 
in diameter, externally as well as internally. In the 
refinishing operation it was necessary to enlarge the 
bore to 0.696 in. compared to the initial value of 
0.630 in. The effectiveness of the seasoning treatment 
is shown by the fact that the vessel has been in use 
for several months, withstanding approximately 100 
applications of pressures as high as 10,200 atm, with 
no measurable increase in diameter of the bore 

The movable leak-proof plugs [6] and associated 
washers are shown in figure 2. Plug A is the type 
used during seasoning and is designed so that the 
washers may undergo large lateral expansions to 
follow the considerable distortion occurring in season- 


ing. Plug B is used in compressibility measure- 
ments. Both plugs are made of Stentor steel, an 


oil-quenching manganese steel, and consist essentially 
of a body and stem each *% in. in length when finished. 
The stem is %. in. in diameter, and is threaded for the 
upper * in. of its length with 24 threads per inch to 
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FIGURE 2 Sections of leak-proof plugs and washers 
permit extraction. It is very important that the 
junction of stem and body be filletted to prevent 
“pinching off’ of the stem by the washers. In 
fabrication, the plug is machined to within 25 mils 
of the desired size, and heat-treated to procure 
extreme hardness (at least Rockwell C-60). The 
hardened piece is ground to size and polished to 
produce a fine finish, which is designed to minimize 
high stresses set up in grooves existing in the ground 
surface. The body of the plug is finished to be 
to 1 mil smaller than the bore of the vessel. Best 
results have been obtained with the smaller clearance 
of 5 mil 

After compression measurements, the frictional 
force between the washers and the bore is of the 
order of 800 atms, and the plug must be extracted 
forcibly. In extraction, large stresses are set up in 
the stem of the plug, and difficulties have been en- 
countered through breaking of the hardened stem. 
The failures appeared principally when the stem had 
been threaded prior to hardening, and were probably 
the result of small cracks in the threads formed in 
quenching. The best procedure has been found to INCH 
grind the threads on the hardened stem. No 
failures to date have occurred in plugs made in this 
manner. It might be advantageous to draw the 
stem somewhat after finishing, but the necessity of 
resorting to this treatment has not been demon- 








Fiat RE; Section of ram 





strated. 

The ram, used to drive the movable plug is shown 
in figure 3, and consists of a cold-rolled steel head 
shrunk onto a glass-hard Stentor steel shaft. The 
head serves to center the ram in the press and is 
relatively unimportant, since all stresses are borne 
by the shaft. The rod contains a recess at its lower 
end, which is finished to provide ample clearance of 
the stem of the plug that fits into it. Considerable 
clearance is advisable to prevent snapping the stem if 
moderate tilting occurs. The diameter of the shaft 
must be such that lateral swelling under compression 
will not cause binding in the bore of the vessel. In this 
instance the shaft is 5 mils smaller than the bore. 
Barring explosions, the ram travel is seldom greater 
than 2 in., and the length of the shaft should be kept 
to a minimum to prevent failure from buckling. 

The pressure-gage plug, the insulated electrical 
lead, and its associated washers are shown in figure 4. 
This is the most complex fitting required and has 
caused more failures (including explosions) than all Ficure 4. Section of gage plug and insulated lead 
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hers combined. The plug and lead are both made 

SAE 3440 steel hardened to Rockwell C—40. 

The electrical lead consists of a rod 4% in. long, 

ing % in. in diameter for 3% in. of its length, 

ving a shoulder 0.262 in. in diameter and *% in. 
ng, and ending in another section %, in. in diameter 

d % in. long. The tip of this latter section is 

proximately > in. in diameter to facilitate silver- 

dering to the gage coil. 

The gage plug is 1% in. in diameter and is threaded 
or 14% im. with 12 threads per inch. The end of the 
wlt is machined to two concentric cylindrical sur- 
faces, respectively, % in. in diameter and 's in. long, 
and 0.615 in. in diameter and \ in. long. These 
surfaces are connected with a conical frustum \ in. 
long, which produces a seal with the packing washers. 
A recess 1 in. deep and %¢ in. in diameter, having a 
flat bottom with no fillet at its wall, is machined in 
the end of the bolt, and a (y-in. diameter hole 
penetrates the remaining portion of the bolt. A 
-in.-deep hole ‘2 in. in diameter is drilled in each 
face of the bolt, and an oversized pin is driven into 
each hole to serve as a convenient point for securing 
electric contact. 

The insulating washers shown around the elec- 
trical lead are in order from the bottom of the 
cavity: \-in. limestone, %.-in. outside diameter, 

-in. inside diameter; \s-in. soapstone, %¢-in outside 
diameter, 29-1 inside diameter; ‘s-in. Neoprene 
rubber, 0.317-in. outside diameter, %%2-in. inside 
diameter; Me-in. steel, 4-in. outside diameter, %%-in. 
inside diameter; ‘-in. soapstone and \-in. limestone 
as before. This assembly is seated by using a small 
jig and an arbor press. 

Two Bakelite washers not associated with the 
pressure seal are useful in preventing short circuits 
of the electrical lead, which may arise frum tilting 
of this lead under pressure. One washer, %2-in. 
outside diameter and %%»-in. inside diameter, serves 
to center the shaft on the low-pressure side, while 
the other, a thin shell 0.312-in. outside diameter 
and 0.262-in. inside diameter, protects the shoulder 
from the plug on the high-pressure side. 

In seasoning, a two-piece fitting was used, one 
being a solid hardened Stentor steel plug having 
the shape of the end of the plug, the other a 1-in. 
bolt to hold the plug in place. 

In operation the worker is protected by placing 
shields around the plug. Although no tendency has 
appeared toward failure of the threads, even at 
seasoning pressures, several failures have occurred 
in which the electrical lead was ejected with high 
velocity. These failures were all traced to errors 
in construction or to failure of steel parts. No 
difficulties, aside from electric short circuits, have 
been experienced with the arrangement described. 
The short circuits have resulted from slow extrusion 
of the limestone insulating washers into the pressure 
vessel, and use of a fine-grained lithographic lime- 

ne apparently has eliminated this difficulty. 


The pressure gage consists of a coil of No. 38DSC 
nganin wire having a resistance of approximately 


}ohms. It is wound in two layers on thin paper, 


shellacked, and baked for 15 hrsat 125°C. Half-inch 
lengths of platinum wire are silver-soldered to the 
ends of the coil and in turn to the insulated electrical 
lead and the pin in the face of the plug. The gage 
coil is mounted on the end of the electrical lead, 
which fits through the coil and provides support. 

The properties and seasoning of manganin pressure 
gages are discussed in detail by Adams, Goranson, 
and Gibson [7]. Gages used in these experiments 
were calibrated by using the water-ice VI transi- 
tion, which occurs at 8,710 atm at 20° C [8]. As 
the pressure coefficient of manganin is essentially 
linear [9], this one point plus the known resistance 
at 1 atm provides a calibration sufficiently accurate 
for the present purpose. 

The freezing-point apparatus consists of a tube 
containing distilled water and a few particles of 
ground glass ‘to facilitate crystallization. The neck 
of the tube is drawn to a fine capillary, and the tube 
is inverted and slipped into a steel tube containing 
several milliliters of mercury. The mercury serves 
to seal the tube and permit transfer of pressure to 
the enclosed water. Calibration is performed by 
subjecting the freezing-point apparatus to pressure 
in excess of the expected freezing pressure, freezing 
being noted as a rapid and sudden decrease in 
resistance of the gage coil. The pressure is then 
lowered until a moderate rate of rise of resistance 
corresponding to melting of ice VI is noted. After 
a period of approximately 1 hr, an equilibrium 
pressure is attained, but extreme care must be 
exercised to insure coexistence of water and ice VI 
at this point. Generally, several equilibrium meas- 
urements are required on freezing as well as melting, 
and such measurements should be in good agreement. 

The electric equipment used for measuring the 
resistance of the gage coil was a G-1 Mueller Bridge. 
This is an unthermostatted Wheatstone bridge of 
the shunted-decade type, having a range 0 to 
51.111 ohms available in steps of 0.0001 ohm. The 
method of measurement differs from that used by 
other principal workers in this field [5, 7] and was 
suggested by E. F. Mueller of this Bureau. 

As the pressure gage is approximately a 100-ohm 
coil, the bridge cannot be used directly to measure 
the resistance of the coil. Two similar gage coils 
are used, with one coil—preferably the smaller of 
the two—being connected directly into one of the 
arms of the bridge where it serves essentially as a 
comparison resistor. This coil is immersed in the 
same liquid as is used for pressure transmission, and 
is mounted in a well in a large metal block attached 
to the side of the pressure vessel. The second coil 
is the pressure coil and is mounted inside the vessel, 
as described previously. One battery lead to the 
bridge is made at the junction of the two coils, as 
shown in the circuit diagram of the electric set-up 
in figure 5. In operation the bridge measures the 
difference in resistance of the two gage coils, one of 
which is subjected to changes in pressure and tem- 
perature—ambient temperature—whereas the other 
is subjected only to similar temperature changes. 
As the resistance of manganin increases with pres- 
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Fieure 5 Schematic diagram o, pressure-measuring equip- 


ment 


sure, the pressure coil should have the slightly larger 
resistance to eliminate the necessity for passing 
through a zero resistance reading. A commutator 
and mereury contact switch were included in the 
circuit but were not used in these experiments. 

When two 1%-v dry cells were used, the total 
bridge current was 12 ma, and the gage coil current 
was 6 ma. With this current a deflection of 's mm 
for each 0.0001-ohm step was obtained on the gal- 
vanometer scale set up at a distance of 2 m. An 
atmosphere was found to be equivalent to approxi- 
mately 0.00026 ohm, and the sensitivity was there- 
fore such that a pressure of less than 1 atm could be 
read directly 

The maximum variation in resistance of the gage 
coils noted at 1 atm was equivalent to +3 atm, 
and might be caused by the effect of temperature 
a maximum temperature variation of + 2° C occurred 
during all measurements—on the bridge or on the 
coils or by the effect of pressure on the gage coil 
resistance. In the absence of precise temperature 
control this variation was ignored, as it represented 
a negligible source of error No difficulties due to 
contact resistances in the switches were ever ex- 
perienced 

No description of the accessories required will be 
riven A complete set of extractors for every part 
is required, but the design of such instruments is 
arbitrary. The extractor generally takes the form 
of a shaft threaded externally at one end and either 
externally or internally at the other. A large nut 
fitting the former end is screwed down against a 
convenient bearing surface and serves to extract the 
desired part, which is screwed into the latter end of 
the shaft. It is of the utmost importance that suit- 
able provision be made for extracting every piece 
before pressure is imposed, for even though parts 


may have large clearances they are often foun 
be solidly joined after a single application of press 

The dial gage used to measure piston travel 
attached to the top platen of the hydraulic p 
This gage had a useful range of 0.400 in. and 
graduated in 0.0001 in. Suitable gage blocks y 
used to extend the range of the gage. Calibrat 
showed a maximum error of +0.0001 in., and 
corrections were applied to the readings 

A press suitable for the measurements was 
available at this Bureau, and through the coop 
tion of the Geophysical Laboratory of the Carn 
Institution of Washington the equipment was 
up for use with the 500-ton, hand-operated hydra 
press previously described [5]. 


Ill. Experimental Method and Calculations 


Measurements consist of readings of the depth of 
the plug in the vessel—as denoted by the dial gag 
and the corresponding pressure [5]. Observations 
are made in order from the maximum pressure at 
every 1,000 atm. ‘To insure uniform behavior of 
the packing washers, all measurements are mad: 
with the plug moving into the vessel. This situation 
is obtained by decreasing the pressure to 50 atm 
below the desired value, and, after a pause to permit 
attainment of thermal equilibrium, pumping back 
as closely as possible to the desired pressure. The 
dial gage readings obtained are interpolated to even 
thousands of atmospheres to simplify the calculation 
This procedure is repeated until the “referen: 
pressure,”’ here 2,000 atm, is reached. Occasionally 
observations are made at 1,000 atm also, but fric- 
tional forces often render these values unreliabl 
A similar measurement is performed with the speci- 
men replaced by a steel bar having similar volum: 
The two sets of data are used to calculate the com- 
pression of the specimen in terms of that of the steel 
bar, which is known [10]. 

The equation required for the calculation was given 
by Adams et al. [5], but the terminology used here 
is that given by Adams [13]. The equation is as 
follows 


k—k’=S Vy) AL l+aP)+Al’aP—P’)4 


In this equation P refers to pressure in atmosph 
and V volume; subscripts zero denote 1 atm, sub- 
scripts r and / refer to steel reference bar and 
fining liquid, respectively, and primes denote 
experimental reference pressure of 2,000 atm; 
the compression (V—V9)/V> AV/V, at pres 
P; k’ is the compression at 2,000 atm (V’—V. 
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term 
term 
arisil 
and 
iq Ul 
madi 


Correction terms used in calculations of compressi- 


bility 


Pressure ! varsol 


atm 

10, 000 00472 0.0130 0.0104 
9, 000 00414 Olly 009! 
&, 000 00352 0104 OO78 
7.000 00296 37 OO91 OO65S 
6, 000 00237 0078 0052 
OOS 
0052 
0039 
0026 
0013 


0039 
0026 
0018 
0000 
—. 0013 


5, 000 OOLTS 
4, 000 OO119 
}, OOO 00060 
» 000 00000 


1, 000 00060 


§ is the cross-sectional area of the bore at 1 atm 
245.88 mm’); @ is the pressure coefficient of this 
area (1/Ay) (€A/dP)=1.3* 10~-*/atm; AL is the differ- 
ence in travel of the piston at pressure P reckoned 
from the reference pressure of 2,000 atm in the meas- 
irements on steel bar and specimen; and Al’ is the 
absolute difference in the position of the piston at the 
reference pressure in the runs on reference steel and 
specimen. 

It is clear that if V’=V! and if Al’=0, three 
terms on the right side of eq 1 disappear. These 
terms then may be considered as correction terms 
arising from inequality of volumes of reference steel 
and specimen and differences in volume of confining 
iquid. The contribution of these terms may be 
made negligible if suitable care is exercised in select- 
ng the volume of the reference steel and in filling 
the pressure vessel. 

The values of the correction terms of the right side 
of eq 1 are given in table 1. The value of @ was 
obtained from application of the equations given by 
the theory of elasticity for deformation of a cylinder 
inder internal pressure [11]. The compression of the 
steel reference bar follows from the compressibility 
of steel given by Bridgman [10] converted into at- 
mospheres. Compression values for Varsol were 
determined experimentally on the Varsol actually 
ised, and were calculated by setting V,=V’=0 in 
q 1 

lhe volumes of the specimens studied are measured 
vy the method of hydrostatic weighing at atmos- 
heric pressure. The materials of interest have 
arge compressibilities, which are functions of 
pressure, and the values of the volumes at 2,000 atm 
equired in eq 1 are not easily obtained by successive 
ipproximations involving successive calculations 
ind extrapolation to the known values at 1 atm. 
lo place the determination of values of V’ on a 
athematical basis, eq 1 was rearranged in the form 


( ! =) i |+s ait +aP)AL+ 


a(P—P’)al’ +( al’ mse. 


Ss l 


1+(k, 


All quantities on the right sid~ of eq 2 are known 
explicitly except V’. 
The compression, —(AV/V__., L us been found to be 


well represented by a cubic i: p.essure, as given by 
AV 7 
—>>- —aP+ bP?+cP"*. (3) 
\ 0 
The right side of eq 2 and 3 are equated, and the 
value of V’/V,) was obtained from the data by the 
method of least squares. The value of V’ so ob- 
tained was used in equation 1 to calculate experi- 
mental values of k—k’. 


IV. Preparation and Treatment of Speci- 
mens 


Measurements were made on three types of 
commercial leather: chrome, retan, and vegetable 
tanned, as well as two types of collagen: hide 
collagen and kangaroo tail tendon. 

Test samples of leather and hide collagen consisted 
of a stack of disks % in. in diameter and 4 in. long, 
held together by means of fine wire. The sample of 
tendon was a bundle of tendons of similar dimen- 
sions. The disks comprising the sample were sub- 
jected to the following treatments prior to test. 

Two samples of chrome tanned hydraulic packing 
leather denoted as A and B, cut from different areas 
of the same hide were tested. The specimens were 
degreased for 24 hr with chloroform, washed in 
running water for 24 hr, and conditioned at 72° F 
and 50-percent relative humidity for 2 mo prior to 
use. 

After measurements on these conditioned samples, 
the Varsol was removed by evaporation at room 
temperature for 1 week, and exposure to a tempera- 
ture of 60° C for 24 hr. The samples were then 
dried at 100° C for 18 hr. Sample A was then tested 
and the Varsol removed again, as described pre- 
viously. Both samples were then redried at 100° C 
for 24 hr in vacuum (<_2mmHg) and retested. 

One sample of retan upper leather was degreased, 
washed, and conditioned as described previously. 
Tests were made only on the conditioned sample. 

Two dried samples of vegetable tanned leather 
obtained from neighboring areas of the same hide 
were tested. One sample, denoted as A, was pre- 
pared by degreasing, washing, and drying at 100° C, 
as described. The second sample, B, was washed, 
dried in vacuum for 16 hr, degreased, and then 
dried at 100° C. 

Partially purified hide collagen was obtained by 
dehydration of a freshly flayed cowhide with acetone. 
The dehydrated material was split to remove the 
grain and flesh layers, and the resulting corium was 
conditioned and used for test. 

Tendons were conditioned 
ceived no chemical treatment. 

All samples were weighed prior to test and were 
then subjected to impregnation with Varsol under a 
vacuum of approximately 20 mm Hg. Following 
compression determinations the samples were re- 
weighed in Varsol. From the loss in weight in 


prior to test and re- 
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ABLE 2 Compressibility of conditioned chrome-tanned leather 
I y 0) 


Sample A Sample B 
1 9.595 mil; 1 1.333 ml 1 9.520 ml; ¥ 9.281ml 
22% 1.367» gen 922% 1.362; g/cm? 
Py 
Ru Ru 

2 ; ‘ Average l 2 A verage 

10.00% 06U3 0. 0701 0. 0674 0. 0899 0. 0692 0. 0697 0. 0702 0. 0700 
000 0630 0641 O618 0643 0633 O44 0642 0643 
SK) O558 OS7T4 O553 O577 OM OS7S 0573 O576 
7,000 O4s81 0495 (487 0493 492 O510 O50 O505 
6,000 O308 0419 0399 0413 0407 0424 O418 0421 
5000 O319 0329 0310 0827 0821 O33 0845 0340 
ALL 221 (222 rat2 0226 (22 PAL Oz 233 
Onn O12 O120 Ola O121 0120 O121 O125 O12 
Jinn (euw) (au) (eu) (ann (nun wan (naw tt 

Alt it sue A nea ilu 
ail. 0. 0888 0. 06% 0. 0708 0. 0682 0. 0602 0.071 0. 0885 0. 0698 
000 632 0630 O82 O82S8 On) O65S OH2S8 4 
5.0K OM O54 OM. O559 O53 0590 OMS O577 
un 402 O482 0492 O48S O480 O51 O4o8 O504 
6,000 O40 040 0406 040 O405 (432 0410 0421 
Om 031 0321 OBS OS18 347 0322 0334 
4.000 (22 218 22 0223 0223 0239 O22 (233 
7 121 O12 121 0120 0120 0122 124 2 
Jian (nan Pe te) (HuM) (aa) (aa (hen (Mum) (aun 
] At ck “ 


Varsol and the density of the Varsol, the volume of mental error the compressions of samples A 
the samples was obtained. The density of Varsol are not significantly different. 
was obtained by successive weighings of a steel bar 


in air, water, and Varsol, respectively, and was 2. Effect of Moisture on Compression 
found to be 0.7811 g/ml 


Measurements were made on sample A al 


V. Results and Discussion drying it at 100° C, and on both samples A an 


after drying them at 100° C in vacuum, as describ 


1. Reproducibility The results are given in table 3. 


In order to evaluate the reproducibility of the TABLE 3 Effect of moisture content on compres 
: chrome-tanned leathe 
measurements, duplicate experiments were made in yme-tanned leath 
many instances. A particular set of measurements 


was obtained in the following way. Sample A 

An initial measurement was made on the steel Sample A dried a — 
reference bar, and was followed by a measurement at 100° ¢ aaa ‘anand 
on sample A and then on sample B. The second Vo=7.882 ml eer CL | ie 
determinations were made by measuring sample A, ae a fs: 4 vacuum whe vat 
then sample B, and finally the reference bar. As as Vo=7.59 “ee 
only two measurements were made in 1 day these Ru Aer 5. a 
determinations required 3 days for completion. The ? 
values of compression could be calculated on the ; 
basis of either of the two measurements on the steel 
bar, and both calculations were performed. It is 
considered that the best calculation is that based on Om os = hae - 
the run nearest in time to the reference run. Values 9,000 0827 of2 Ori42 

" &.000 0570 0568 O5T2 

based on the farthest removed reference run are 7000 0502 O477 0494 
given as alternate values. The results of this experi- — a4 poe amas = 
ment are given in table 2 4,000 0240 023! 02z9 

Runs 1 and 2 of the table refer to the experiment 4 po oan pon wn 
described, and two other runs performed at a different 1,000 o119 o11s 0134 
time on sample A are included for comparison. The 
agreement is good between runs, and within experi-  2¢on arm “1 © The Felative volume change, —A V/V, reek 
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It is surprising that the removal of moisture causes 
ttle or no change in compression of the leather. <A 
ital loss in weight of 18.6 percent is involved between 
onditioned samples of table 2 and the dry samples of 
able 3. If it is assumed that this total loss in weight 

presents moisture [12], and the compression of a 
omposite sample containing 18.6 percent of water 
nd 81.4 percent of leather is calculated from the 
ompressions of table 3 and known values for water 

3], values are obtained that are of the order of 10 
ercent higher than the measured values of table 2. 

kxven if the moisture present in the leather is con- 
sidered to be “‘bound water’ under high pressure, it 
would appear that the compression of the leather 
should change upon removal of the water, except in 
the one case in which the compressibility of the 
combined water and the leather are nearly identical 
over the whole pressure range. The probability of 
such coincidence seems small. This finding deserves 
irther experimental study on this and other hygro- 
scopic materials containing varying amounts of 
morsture, 


3. Compressions of Other Types of Leather and 
Collagen 


The results of measurements of compressions of 
two types of collagen and retanned and vegetable- 
tanned leather are given in table 4. 

It is to be noted that the compression values of the 
vegetable-tanned samples agree but appear to be 
slightly lower than values obtained for chrome- 
tanned samples. Values for the retan leather, on 
the other hand, appear somewhat higher than 
values for chrome leather, while the values obtained 
for the collagens agree well both with each other 
and with the chrome leather. The divergences 
observed with vegetable and retanned leather may 
be real, but conclusions in this respect will not be 
made until further more extensive measurements 
are conducted 


4. Compressibility of Leather 


To obtain an explicit expression for compression, a 
power series in pressure may be fitted to the experi- 


TABLE 4 


Dry vegetabie-tanned crust leather 


Retan 
conditioned Hide collagen Kangaroo tendon 
Sample 1 Sample 2 
1 9.380 ml J 8.796 ml \ 6.311 ml 
Pressure Vo=9.536 ml Vo=9.071 ml 9.033 ml \ 8.556 ml V’ =6.139 ml 
\ 9.297 ml J 8.840 ml pn°c = 1.3105 g/cm! 922°C = 1.32ts g/cm p22°C = 1.3286 g/cm 
22°C = 1.4065 g/cm p22 1.39%. g/cm 
ak—k k—k . 
um 
10,000 0. O67 0. 0660 0. 0748 0. 0701 0. 0706 
000 On08 OS7e OnS4 of29 Oo45 
&,000 0538 0530 0610 0566 0570 
7,000 0465 0467 0524 O4AN4 0491 
6,000 0387 0890 0449 O414 0407 
5.000 a2u9 0300 0354 O86 (B17 
4,000 0219 0216 0257 0236 0234 
8,000 ola O117 0142 0122 0126 
2,000 0000 0000 0000 0000 0000 
k’ is equivalent to the relative volume change, —A V/ Vo, reckoned from 2,000 atm 
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Compressibilities of collagen and leather 





























mental data. Terms of higher degree than the 
third in pressure have been found to contribute 
little to the efficacy of the fit, and a cubic equation 
in pressure appears to describe the data satisfactorily. 
Data obtained for different samples, tannages, and 
moisture contents might all be fitted with such 
equations, but since there are only slight differences 
in compressions of these materials the equations 
evolved would be expected to agree fairly well and 
to exhibit small variations in coefficients. Therefore, 
an equation of third degree in pressure was fitted 
only to the best data—the average values for sample 
A in table 2—-which are considered as representative. 
Values obtained by use of this equation will closely 
approximate the experimental values obtained on all 
leathers tested. The equation obtained by use of 
the method of least squares is as follows 


—AV/V, = 1.231075 (P-2000) — 5.60 « 10-" (P- 
2000) ? + 1.35107" (P-2000) *. (4) 


The compressibility obtained by differentiating the 
above equation is found to be 


_1 dV _) 93% 10 


dP 1.12K10-°C 


2000)+ 


4.05 107" (P— 2000)’. (5) 
A previous estimate of the compressibility [4], 33x 
10~°/bar, is considerably larger than that given by 
the equation above. However, equation 2 shows the 
compressibility to be a function of pressure and the 
value of 3310~°/bar was obtained at an undeter 
mined pressure much lower than the minimum pres- 
sure of 2,000 atm at which eq 2 applies. Extrapola- 
tion of eq 4 to atmospheric pressure yields a com- 
pressibility considerably smaller than 33> 10~°/bar. 


VI. Conclusions 


The design and construction of appar< tus for meas- 
uring compression of solids (or liquids) in the pres- 
sure range 2,000 to 12,000 atm has been described. 
The results of measurements made on a natural high 
polymer, leather, show that the apparatus is capable 
of vielding data of good reproducibility. 





The data obtained on leather are indicative of a 
compressibility that is intermediate between the 
compressibilities exhibited by most solids and those 
of liquids. The compression is not materially af- 
fected by moisture content. Values obtained for 
chrome-tanned leather and two types of collagen 
are in good agreement, but the compressions of vege- 
table and retanned leathers appear to be slightly 
different. Compressions of samples taken from dif- 
ferent areas of the hide are in substantial agreement. 


a 


for advice and assistance in connection with th : 
sistance measurements. 
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Arc and Spark Spectra of Technetium co 
nt 
“115 : att 
By William F. Meggers and Bourdon F. Scribner “: 
re 
Four milligrams of highly purified technetium, loaned by the United States Atomic da 
Energy Commission, were used to obtain a description of the are and spark spectra charac- ide 
teristic of this fission product. Solutions containing 50 to 200 micrograms of technetium ing 
were dried on copper electrodes and excited by electric ares or sparks. A stigmatic concave | 
grating of 22-feet radius was employed to photograph the spectra from 2200 to 9000 
angstroms, within which limits more than 2,300 lines characteristic of Te atoms or ions art 
were recorded. Wavelengths were measured relative to iron standards, relative intensities Co 
were estimated on a scale of 1 to 1000, and almost every line was definitely assigned either to 
to neutral Te atoms or to singly charged (Te*) ions. The measured wavelengths range a0 
from 2261.30 to 8829.80 angstroms with average probable errors rarely exceeding +0.01 — 
angstrom. The average results of wavelength measurements and of intensity estimates wa 
for 2,121 lines in Te 1 and Te 1 spectra are presented. The strongest Tc 1 lines have wave- pre 
lengths 4297.06, 4262.26, 4238.19, 4031.63, and 3636.10 angstroms. The strongest Tem an 
lines have wavelengths 2543.24, 2610.00, and 2647.02 angstroms. This description of Ter Th. 
and Tem spectra will serve for spectro-chemical identification and for structural analyses 
of these spectra Spe 
are 
I. Introduction finding, in material extracted from columbite, thre: Am 
K spectrum lines for atomic number 43 and five / thr 
In 1869, when the periodic chart of the atoms was spectrum lines for atomic number 75. These two for 
first proposed, Mendeléeff [1]' saw that the chemical atomic numbers are identical with the homologs o! and 
elements homologous to manganese were missing in manganese, and the names masurium and rheniun fine 
the two succeed*ng periods; he called these eka- were proposed for them. The discovery of rheniun of 
manganese and dwi-manganese. This was the first was abundantly confirmed, and within a few years ela 
hint of the possible existence of the element under this new element became a commercial metal. |! istic 
discussion. In the following half century many 1931 the are spectrum of rhenium was described and for 
attempts, both chemical and spectroscopic, were interpreted by Meggers [4], who then inquired of th ana 
made to detect these homologs of manganese, but discoverers if a sample of masurium was availabi the 
without success. Moseley’s [2] discovery in 1913 ofa for spectroscopic investigation. The Noddachs nd 
quantitative relationship between Roentgen fre- replied that they had succeeded in concentrating spec 
quencies and atomic numbers proved that atomic some masurium, and that they woulda .nd the firs! 
numbers 43, 61, 72, 75, 85, and 87 were still unknown, milligram that could be spared. Unfortunately that 
and provided a new method of detection by means of milligram never materialized, and it later appear ' 
Roentgen spectra. This method was applied in 1925 that the extraction of element 43 from natura 
, ah . . . . - I e 
by Noddack, Tacke, and Berg [3], who reported sources would be unlikely or even impossible (5) all 


Figures in brackets indicate the literature references at the end of this paper 
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? In a letter dated Dec. 22, 1931, from W. and I. Noddak to W. F. Me 
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\fter the discovery of artificial radioactivity in 
134, it became theoretically possible to manufacture 
iy desired element by artificial nuclear transmuta- 
on. Thus element 43 was first prepared in 1937 

the Berkeley cyclotron by bombarding molyb- 
num with neutrons or deuterons. It was iden- 
fied, and its radiochemical properties were in- 
estigated by Perrier and Segré [6], who named the 
ement technetium [7] because it was the first new 

element to be produced artificially. The amount 
of technetium produced in this manner was so small 
that it could only be detected by means of its 
characteristic radioactivity, which is enormously 
more sensitive than either roentgen or optical 
spectra. However, soon after the discovery of 
uranium fission in 1939, it was found that element 
{3 occurs among fission products io the extent. of 6.2 
percent [8]. In 1948 the isolation of milligram 
amounts of element 43 from uranium fission was 
described by Parker, Reed, and Ruch [9]. With 
these samples the whole-number atomic mass was 
determined by Inghram [10] to be 99, the first optical 
spectrogram was made by Timma [11], and the 
Roentgen emission spectrum was photographed by 
Burkhart, Peed, and Saunders [12]. Thus far the 
published information on the optical spectra of Te 
consists of estimated relative intensities and approxi- 
mate wavelengths of 102 lines (4297.2 to 2461.7 A) 
attributed [11] to element 43. In the light of 
present information critical examination of those 
data discloses that possibly 46 of the 102 lines are 
identifiable as Ter, and 32 as Teu, but the remain- 
ing 24 do not belong to Te. 

In September 1948, the National Bureau of Stand- 
ards applied to the United States Atomic Energy 
Commission for the loan of 5 mg of {Te with which 
to make descriptions and analyses of the atomic and 
ionic emission spectra of technetium. A 3-mg sample 
was delivered in January 1949. With this sample 
preliminary results were obtained for wavelengths 
and relative intensities of about 2,000 Te lines. 
These results were reported at the Oak Ridge 
Spectroscopy Symposium on March 25, 1949 [13] 
and at the Buffalo meeting of the Optical Society of 
\merica on October 29, 1949 [14]. In July 1950 
three additional milligrams of {Te were borrowed 
for the purpose of completing this description of are 
and spark spectra and for the investigation of hyper- 
fine structure of certain spectral lines. The purpose 
of this paper is to present the wavelengths and 
relative intensities of 2,121 lines that are character- 
istic of Te atoms and ions. These data are suitable 
for spectrochemical identification and for structural 
analyses of the first two spectra of Te. Details of 
the structural analyses of Ter and Tem spectra, 
and of the investigation of hyperfine structure in Te 
spectral lines will be given in subsequent papers. 


II. Technetium Spectrograms 


lhe technetium samples came in the form of oxide. 
[hese samples were dissolved in dilute ammonium 
livdroxide to make solutions of known concentration, 
ruions of which could then be added to spectro- 


scopic electrodes so that a definite amount of Te was 
deposited on each. Fifty or one hundred micro- 
grams of Te were dried on the flat-machined end of a 
copper rod \ in. in diameter and 1 in. long. Elec- 
trodes with approximately rectangular ends, 6 by 
2mm, were later employed to concentrate more of the 
light along the optical axis. Two such loaded elec- 
trodes were used when oscillatory discharges were 
employed for exciting the spectra, but when the 
direct-current arc was employed the best results 
were obtained by placing 200 yg on the anode. With 
a desire to economize material, the first spectrograms 
were made in the blue and ultraviolet with 50 yg on 
each electrode, but to obtain stronger spectrograms 
it was necessary to use 100 on each. Also it was 
desirable to choose light sources that would burn 
the entire added sample off the ends of the copper 
electrodes without predominant excitation of the 
copper spectra. Furthermore, it was necessary to 
choose two sources with sufficiently different excita- 
tion characteristics so that comparisons of spectral 
line intensities would permit positive assignment of 
the lines either to neutral atoms or to ionized atoms. 
These considerations led to the choice of two sets of 
excitation conditions available in the Applied Re- 
search Laboratory Multisource unit. 

The spark spectrum was obtained from an oscil- 
latory condenser discharge with the following circuit 
constants: capacitance 5 uf, charging voltage 940 v, 
inductance 100 wh, and resistance 0.4 ohm. The are 
spectrum was obtained from an overdamped con- 
denser discharge with circuit constants: capacitance 
50 uf, charging voltage 940 v, inductance 400 yuh 
and resistance 25 ohms. The gap between the 
electrodes was set at 4 mm. 

Before any Te spectrograms were attempted, 
experiments were made with Mo spectrograms, and 
source conditions were determined that gave a five- 
fold enhancement of Mou lines over Mor lines 
when the spark spectrum was compared with the are. 
These source conditions described above were found 
to be suitable for an unambiguous separation of 
Tet and Tew lines. 

It was immediately seen that the distribution of 
Terand Ter lines resembled the distribution of first 
and second spectral lines of other metals; the second 
spectrum predominates in the ultraviolet and weakens 
in the visible, whereas the first spectrum is weak in 
the ultraviolet but predominates in longer waves. 
Because few Ter lines are seen in the blue, the Te 
spark spectrograms were not extended beyond 
4700 A, 

Although the are-like excitation provided by the 
ARL Multisource unit was satisfactory for the short 
wave region, it was disappointing for the visible and 
infrared, partly because it gave weak spectra and 
partly because it excited the spectra of atmospheric 
nitrogen and oxygen with such intensity that many 
portions of the Te spectra were masked. Conse- 
quently, the spectral range 4500 to 9000 A was 
reobserved with a direct-current are between ‘s-in. 
diameter Cu electrodes, the anode being charged 
with 200 we of Te. These direct-current are spectro- 
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grams were made with an applied potential of 220v 
and a current of 7 amps 

The spectrograms were made with a concave 
grating of 22-ft radius, ruled with 15,000 lines per 
inch on an aluminized pyrex disk of 6-in. diameter. 
This grating was placed in a Jarrell-Ash instrument 
provided with a concave mirror of similar form, so 
that the spectrograph gave stigmatic slit images. 
A fixed slit of 20-u width was used, and the maxi- 
mum resolving power for sharp lines in the first 
order spectrum was about 50,000, the reciprocal 
linear dispersion being about 5 A/mm. 

Three positions of the grating and 20-in. plate 
holder recorded the spectrum between 2100 and 
9000 A with successive exposures on three types of 
photographic plates. The region 2200 to 4700 A 
was photographed on Eastman 103-uv_ sensitized 
plates, the region 4400 to 6900 A on II-F sensitized 
plates, and 6500 to 9000 A on I-N hvpersensitized 
plates. 

By sliding a slotted diaphragm in front of the slit, 
a series of adjacent spectrograms were made as fol- 
lows: Iron are, copper arc, technetium are on copper, 
technetium spark on copper, copper spark, iron are. 
Spark spectra were omitted in the visible and near 
infrared regions. The exposure durations were 
usually 1 sec for the iron arc, 45 see for the copper 
blanks, and 30 sec for copper with technetium added. 
Preliminary tests with Mo added to copper showed 
that the sample was almost completely burned off in 
All exposures with Te deposited on Cu were 
made with the electrodes mounted in a cylindrical 
brass cell provided with a quartz window and with 
an air exhaust through a glass-wool filter designed to 


30 sec 


At least three satisfactory spectrograms \y 
made for each of the three spectral ranges mentio: 
This gave at least six independent determinati 
for lines appearing in the overlapping portions 
these ranges. The positions of Te spectral | 
relative to iron standards were measured to 0.001 mo 
with an excellent 50-cm comparator constructed 
the Gaertner Scientific Corp. In all, approximate. |, 
20,000 spectral line images were bisected to obt n 
the wavelengths of about 2,500 Te and impurity 
lines. All lines appearing in the spectra of Cu plus 
Te that were not duplicated or masked in the Cy 
comparison spectra were measured, even though 
some were easily recognized as impurities. Th 
final values of all the impurity lines served as checks 
on the seale of wavelengths and on the amazingly 
high purity of the Te samples. The principal impur- 
ities detected spectroscopically were calcium, mag- 
nesium, and silicon, which are common contaminants 
in chemical operations. Only the stronger lines of 
lithium, sodium, and potassium were present, and 
the same is true of iron, nickel, chromium, molybde- 
num, aluminum, bismuth, and platinum. No trac 
of ruthenium, rhenium, or any other metals w 
found. In the identification and elimination of im- 
purity lines, consideration was given to relative 
intensities as well as to wavelengths. Thus, the line 
2496.77 A, with intensity 200 in Te, cannot be con- 
fused with boron, 2496.778 A, because the stronger 
line of boron, 2497.733 A, had only intensity 10. A 
few doubtful coincidences of Te and impurity lines 
are retained in table 1, and it is not likely that many 
Te lines have been masked by impurities. The chief 
omissions may be real Te lines hidden behind the 





as 


trap, and thus conserve, most of the vaporized Te principal lines of copper, 3248, 3274, 5106, 5153 
for future recovery 5218, 5700, 5782, 7933, and 8093 A. 
rasLe | Are and spark spectra of technetium 
= wavelengt! mplex: d uble: h=hazy; //=very hazy: w=wide: W=very wide. s=shaded to shorter waves; /=shaded to longer wave 
» Intensity \ Intensity r Intensity N Intensity | 
1 Are Spark A Arc Spark A Are Spark A Are Spark 
S820. 80 10% S308. 14 10 7TO81L. 86 l 7779. 53 l 
S751. 80 l 8254. 48 2 7968. 27 2 7777. 17 3 
8737. 95 Seclt S248. 05 3 7965. 41 5 7746. 60 2 
8722. 70 2 8237. O8 de 7963. 86 l 7737. 36 3 
8719. 97 2 8228. 65 l 7949. 64 2 7705. 70 l 
8716. 77 l 8225. 09 36 7885. 96 2 7697. 36 RO 
8707. 21 st) 8211. 29 10 7874. 72 5 7684. 43 s 
8652. 76 3c 8206. 49 } 7871. 16 50 7681. 73 3 
8543. 63 ~ 8205. 24 10 7863. 11 l 7675. 32 l 
8537. 66 2 8170. 51 10 7861. 36 30 7658. 31 2 
8531. 07 7 8160. 89 2 7858. 24 2 7624. 49 7 
8522. 56 2 8153. 52 l 7856. 32 7 7618. 39 3 
8514. 70 2 8146. 52 l 7854. 35 l 7579. 20 10 
8481. 37 l 8142. 60 2 7848. 35 le 7573. 99 15 
8404. 00 l 8127. 78 l 7817. 67 150¢ 7550. 98 l 
8383. 68 3 S119. 87 3 7816. 72 5 7548. 40 i 
8346. 59 hi 8068. 80 2 7807. 11 l 7545. 45 3 
8314. 09 2 8045. 30 2 7798. 18 4 7543. 32 7 
8315. 49 H 8042. 45 2 7792. 96 100cW 7540. 19 | 100¢ 
8300. 17 10cu 7099. 68 10 7782. 92 2 7534. 87 6 
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TABLE 1 Arc and spark spectra of technetium—Continued 


r Intensity Intensity d Intensity r Intensity 

1 Ar Spark A Arc Spark A (re Spark A Are Spark 
1790. 47 3 1616. 86 8) 5 4328. 77 3 2 1146. 21 l 
1789. 48 2 4609. 17 te oh $324. 85 l $145. 02 SOx 50 
1787. 69 2c 1608. 00 2 l $323. 94 l 1141. 28 3 2 
1785. 59 1 1602. 72 2 l 1323. 35 3 2 $139. 85 5 i 
1784. 69 l 1594. 06 2 l 1319. 94 lh 1139. 10 3 2 
1783. 92 2 1593. 36 20 10 1318. 65 2 2 $134. 8O 3 2 
1780. 56 l 1584. 86 2 l 1314. 03 l $130. 44 2 I 
1779. 84 l 1579. 09 l 1312. 48 l $128. 27 10 6 
1776. 33 2 1578. 46 10 6 $310. 21 3 2 $124. 22 SO 50 
1773. 86 2 1570. 85 2 l $308. 59 2 l $119. 28 } 3 
1771. 53 150 1564. 57 30 20 1307. 78 l l $1115.08 100 60 
1762. 35 2 1563. 80 I l 1305. 82 } 3 $113. 26 l 
1758. 06 l $562. 94 l $305. 52 l $111. 37 2 l 
1756. 09 l 1558. 73 l $297. 06 500 300 4110. 21 10 6 
1752. 71 1 1557. 06 6 } 1294. 35 5 3 1106, 69 2 l 
1749. 60 3 $552. 85 5 } $289. 54 2 l 1106. 21 l 
1746. 02 | $552. 41 l 1286. 31 l $105. OS 2 l 
1741. 24 I $552. 21 3 2 $278. 91 10) 20 4101. 24 2 l 
1740.59 200 $545. 56 l 1277. 53 l 4100. OS 2 | 
1739. 53 l 1542. 10 ! 3 $274. 95 3? 2? 1098. 77 l 
1736. 48 % 1539. 54 100 50 1266. 16 l l 1095. 68 200 100 
1719. 28 100 $1537. 36 2 l 1262.67 100 50 1094. 17 2 | 
1717. 76 10 $522.85 200 100 1262. 26 400 200 1093. 68 3 2 
1716. 76 l $515. 98 10 6 1258. 61 | 1091. 54 2 l 
1714. 20 2 $501. 74 l $257. 33 l 1090. 13 2 l 
1706. 90 6 $498. 54 lh $243. 00 2 2 1089, 23 l 
1702. 83 l $495. 03 10 5 $239. 55 2 l 1088. 70 150 80 
1694, 28 3 $487. 05 10 20 1238.19 300 150 1083. 54 3 2 
1689. 77 l $481, 52 10 5 $232. 09 2 l 1083. 20 l 
1689. 34 j $473. 65 l $230. 84 6 } 1079. 81 lh lh 
1687. 97 2 $465. 61 2 l $230. 22 } 1077. 71 2 l 
1686. 37 l 1463. 82 2 l 1229. 96 ! 1069. O1 2 l 
1683. 58 | $454. 81 5Ca? j $220. 32 2 l 1060, 84 l 
1678. 88 2 $452. 24 2 l $218. 61 30 15 1057. 81 | 
1673. 05 2 41434. 98 2 ] $215. 51 | | 1057. 24 2h 1} 
1672. 82 2 1432. 54 l $214. 22 l 1056. O8 th Sh 
1672. 17 9 $429. 60 30 15 $212. 42 l 2 1054. 42 2h lh 
1669.30 100 10 $425. 47 lh 2h $211. 43 2 l 1053. 18 5 3 
1664. 32 2 l $420. 32 1 $206. 58 9 6 1051. 94 5 3 
1660. 21 30 10 $412. 02 2 l $201. 86 2 l 1049. 10) 150 SO 
1656. 99 | 1383. 10 2 l 4199. 46 2 l 1041. 76 2 l 
1654. 64 2 l 1373. 97 2 l 4198. 31 l th 1041. 52 2 l 
1653. 59 l 1370. 73 2 l 189. 61 3 2 1039. 23 10 6 
1648. 34 70 20 4370. 10 2 l 1186. 52 50 30 1037. 09 3/ 
1647. 57 2 I 1369. 26 l 1176.28 100 60 1037. O1 th-+@ 
1646. 61 l $368. 92 ] $174. 46 l 1032. 91 l 
1645. 47 2 l 1363. 03 2 l $172. 53 80 50 1031. 63 300 200 
1643. 29 j 1361. 95 l 4170. 28 50 30 1029. OS | l 
1637. 51 90 10 $359. 25 j 3 1169. 67 10 6 1026. 19 3h+a 2h 
1635. 74 2 l 1358. 49 ] 2 1167. 43 6 } 1025. 42 l 
1633. 13 2 l $344. 58 2 l $165. 62 80 50 1020. 77 10 7 
1631. 33 l $342. 22 1 1161. 53 l 1017. 23 / 3 
1630. 57 20 10 1336. 87 10 5 1160, 02 ] 1016. 68 3 | 
1624. 97 2 l 1332. 46 3 2 $153. 64 | 1012. O1 6 | 
1622. 71 3 2 $331. 81 l $147. 61 2 l 1010. OS l 
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2394. 49 l 2368. O1 l 
2393. 42 l 2367. 67 l 
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Ill. Results 


The results of this investigation are presented in 
table 1, where wavelengths and relative intensities 
of 2,121 spectral lines characteristic of technetium 
atoms and ions are given. In this table the wave- 
lengths represent the averages of two to six inde- 
pendent determinations Several hundred faint 
lines that were measured only once on the strongest 
spectrograms have been omitted from table 1, but 
some of these may be reported later if any confirm 
predictions from atomic energy levels. 

The average probable error in these wavelengths is 
believed to be less than 0.01 A There are three 
reasons for this belief, namely, computed probable 
errors of lines measured on six spectrograms, agree- 
ment of the wavelengths of impurity lines with their 
accepted values, and the average difference between 
the measured values of classified lines and their 
values computed from atomic energy levels derived 
from these observations [14] 

In addition to six-figure wavelengths, estimated 
relative intensities of are and spark lines of Te are 
given in table 1. These intensities range from 1,000 
for the strongest line to 1, which represents the weak 
lines whose reality is unequivocal because they were 
observed two or more times. In the spectral range 
in which are and spark intensities were compared 
(4700 to 2200 A) Tet lines have nearly double 
intensity in the are, whereas Te lines have about 
treble intensity in the spark. Most of the Tem and 
Tet line-intensities thus differ by a factor of 5 or 6, 
and their assignment is unambiguous. <A few lines, 
mostly faint ones, appear to have about the same 
intensity in are and spark, but naturally a few approx- 
mate coincidences of Ter and Tem lines may be 
expected 

Occasionally the intensity numbers are 
panied by literal symbols having the meanings given 
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\ Intensity Xv Intensity 
1 Are Spark A Arc Sp 
l 2316. 72 l 
l 2315. 32 2 
l 2313. 91 l 
l 2300. 97 l 
2 2300. 44 l 
l 3 2299. 96 l 
l 3 2298. 14 I 10 
l 2292. 14 l 
2h 2291. 64 | 

2 2285. 47 l 
l 2282. 41 | 
l , 2280. 33 2 
l 2272. 45 l 
l 2271. 14 ld ! 
l 2266. 24 2 

2261. 30 


in the table. Some of these symbols are especial! 
significant for Te spectra because they suggest |) 
perfine structure due to interaction of valence elec- 
trops with atomic nuclei. This hyperfine structur 
was observed in manganese spectra and in rhenium 
spectra, unresolved by gratings in the former cas 
but easily resolved in the latter. The hyperfin 
structure of Te lines appears to be intermediate bi 
tween Mnand Re. Many Tc lines, especially among 
longer waves, appear as flat-topped images, four or 
five times as wide as the slit. These lines are unde: 
investigation with interferometers. 

Besides a general resemblance of Te spectra wit! 
those of Mn and Re, we can now define in detail th: 
analogous spectral lines of the three elements o 
theirions. For neutral atoms of Mn, Tc, or Re the reso- 





nance lines are known as d° s* °S,,,—d°s' p' ®P 5,., 

and the respective wavelengths are 4030.76 
1033.07, 4034.49 A in Mn; 4297.06, 4262.26, 4238.19 
A in Te: and 3460.47, 3464.72, 3451.88 Ain Re. For 


singly ionized atoms of Mn, Tc, or Re the resonance: 
lines are known as d° s'*S;—d° p' *P4 3,0, and the r 
spective wavelengths are 2576.10, 2593.73, 2605.69 A 
for Mnt; 2543.24, 2610.00, 2647.02 A for Ti 
but not vet revealed for Ret Full details con- 
cerning the interpretations of lines in the first 
two spectra of Te will be divulged in another 
paper. The above-mentioned resonance lines of ‘T’ 
and Tc* are probably useful for detecting this el 
ment in concentrations of 1077; the fact that they 
cannot be found in any spectral tables excluding 
fission products may be regarded as proof that tech- 
netium is not spectroscopically detectable in na 
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Dissociation Constants of 4-Aminobenzophenone Calcu- 
lated From Ultraviolet Spectra at Several Temperatures 
By Elizabeth E. Sager and Iris J. Siewers 


Many of the physical-chemical properties of compounds of low solubility are difficult 
to determine because of the necessarily small concentrations of material. In many cases 
advantage may be taken of measurements of the absorption spectra if the bands occur in 
favorable ranges of the spectrum. The diphenyl! ketones fall in this category and have been 


the subject of recent studies 


The ultraviolet absorption spectra of several concentrations of benzophenone, 4-amino- 
benzophenone, and 4,4’-diaminobenzophenone have been measured and their molar absorb- 
encies determined. The dissociation of 4-aminobenzophenone, a weak base, has been care- 


fully studied at temperatures from 10° to 40 


C in 5-degree steps. Dissociation constants 


for each temperature, as well as activity coefficient terms, have been calculated. The heat 


of dissociation has also been determined 


I. Introduction 


The acidic or basic properties of many slightly 
soluble compounds have not been determined because 
of experimental difficulties in applying the usual 
titration or electromotive-force methods to very low 
concentrations. Diphenylketones, such as benzo- 
phenone, 4-aminobenzophenone and 4,4’-diamino- 
benzophenone, are practically insoluble in water and 
fall in this category. They are colorless in aqueous 
solution, but their aromatic structure suggests that 
their ultraviolet spectra may be used to study some 
of their chemical properties. Accurate spectro- 
photometric measurements of concentrations as low 
as 10-* or 10°° molar can now be made at many 
temperatures. The observed changes in spectral 
absorbancy with changes in hydrogen-ion concen- 
trations can be used to calculate dissociation con- 
stants and other related thermodynamic quantities. 

The molar absorbancies of the molecular and ionic 
states of benzophenone, 4-aminobenzophenone, and 
+,4’-diaminobenzophenone were first determined and 
compared. The parent compound, benzophenone, 
has no group to ionize, and its spectral curves are 
identical whether the material is in aqueous, acid, or 
basic medium. 4-Aminobenzophenone is a weak 
base and does not dissociate appreciably in aqueous 
solution, as demonstrated by identical curves of the 


free base in water and in various concentrations of 
alkali. Moderately strong acid is required to effect 
complete dissociation, as evidenced by the entirely 
different absorbancy curve in the near ultraviolet. 
4,4’-Diaminobenzophenone has two groups whose 
dissociation ranges overlap, and extensive experi- 
mental data will be required to determine the 
constants 

A comprehensive study of the dissociation of 4- 
aminobenzophenone has been made at several tem- 
peratures, and the constants at 10°, 15°, 20°, 25 
30°, 35°, and 40° C have been calculated and are 
reported in this paper. Hydrochloric acid was used 
to control the hydrogen-ion concentrations at various 
stages during the dissociation. The contribution of 
the 107°-M base to ionic strengths of 107° and 10-?-M 
hydrochloric acid is so small that it can be neglected 
thus the hyvdregen-ion concentration may be taken 
as that of the known amount of hydrochloric acid 

The limiting spectral absorbancy curves of the free 
base and of the ionized compound do not change with 
change in temperature in the range studied, which is 
very fortunate. At a given temperature, then, all 
changes in absorbancy for one dissociation series may 
be attributed solely to the dissociation reaction 
Another favorable factor is that 4-aminobenzo- 
phenone dissociates over a pH range of 1 to 3, where 
the temperature has little influence upon pH values 
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For instance, the pH of 0.005-M hydrochloric acid 
on a comparative scale at 10° C is 2.332, and at 
10° C it is 2.334, a negligible difference. 

The concentration constants for the dissociation 
are determined in the following manner. The degree 
of dissociation and resulting ratios of free base to 
ions are calculated from the absorbancy data. The 
hydrogen-ion concentrations are known from the 
amounts of hydrochloric acid in the solutions. The 
products of the hydrogen-ion concentration and the 
give the concentration 
constants, A By plotting -log,A,. as a function 
of ionic strength, practically straight lines with 
slight slopes are obtained that may be extrapolated 
to infinite dilution for estimated -log,A,, or pK, 
values for each temperature. The activity coeffi- 
cient terms for the weak base can also be expressed 
as a function of ionic strength. Employing the 
usual thermodynamic equations, the heat of dissocia- 
tion, A/7, may also be calculated. 


ratios of free base to ions 


II. Experimental Procedure 
1. Materials 


Benzophenone was obtained from Eastman Kodak 
Co., and 4-aminobenzophenone and 4,4’-diamino- 
benzophenone were obtained from Dow Chemical 
Co. The benzophenone, a white crystalline com- 
pound, melted at 48.6° to 49.2° C. The 4-amino- 
benzophenone, a pale yellow crystalline compound, 
melted at 123.0° to 123.4° C. The 4,4’-diamino- 
benzophenone, a yellow crystalline material, melted 
at 244.1° to 244.4° C 

Hydrochloric acid of reagent grade was used to 
make stock solutions, 10~'-./, for the majority of 
the dilutions. Stronger acid was required in some 
ceases for determining limiting absorbancy curves. 
Sodium hydroxide, approximately 10~'-./, was pre- 
pared free from carbonate, from which weaker con- 
centrations of alkali were made. Conductivity 
water was used to dissolve the compounds under 
investigation and in making all stock solutions and 
subsequent dilutions. 


2. Equipment 


A model DU Beckman spectrophotometer was 
modified with a constant-temperature cell compart- 
ment of our own design.' The cell compartment is 
essentially a metal-box air bath, jacketed on all 
sides, top, and bottom, through which water cir- 
culates from a constant-temperature bath” It was 
designed especially to accommodate large eylindri- 
cal absorption cells of 38-mm inside diameter, 
described in previous publications [1, 2, 3]. The 
cells are Pyrex cylinders of 1-, 2-, and 5-em length 
with openings for filling and emptying and are fitted 


with removable end-plates of crystalline quartz. 
The constant-temperature cell compartment was designed with the aid of 
W. 8. Hertstein, draftsman, of the Instrument Shop of this Bureau, who made 
many helpful suggestions. I, W. Lund, also of the Instrument Shop, constructed 
the metal box, metal holders, and other parts necessary to modify the instrument 
A I-gallon constant-temperature bath of American Instrument Co. design 
was used 
Figures in brackets indicate the literature references at the end of this paper 


The cell assemblies are held together with meta 
containers and screw caps, using Bakelite an 
rubber gaskets. All cells, end-plates, and meta 
parts were made at this Bureau. One complet 
cell assembly measures about 44 mm in over-a 
diameter and 60 mm in over-all length. To alloy 
passage of the radiant energy from the exit sli 
through the cells and solutions to the photoce 
compartment, two round openings about 22 mm j 
diameter were made in opposite sides of the wate: 
jacketed box. Only two cells can be accommodated 
Each cell is held in a metal cradle attached to a 
sliding red, by which mechanism they can be placed 
alternately in the path of radiant energy. Th 
metal housing with the source of energy was moved 
back about 1 in. so that it is not in contact with 
the metal walls of the cell compartment. 


3. Temperature Control 


maintained at 


The water bath was the desired 
temperature, between 10° and 40° C within +0.| 
deg C. Water, initially at the bath temperatur 
circulates through four lengths of hose attached to 
the two inlets and outlets of the cell compartment 
and cover in the room at 25° C. Consequently, the 
temperature within the cell compartment is not 
exactly that of the water bath when it is above or 
below 25° C. The temperatures reported are read 
ings taken from a thermometer with only the 
mercury-bulb tip inserted in a rubber stopper in the 
removable top. Care was taken not to allow more 
than the tip of the thermometer inside the box during 
spectrophotometric observations, so that extraneous 
light could not be introduced. In order to shorten 
the time for the solutions in the cells to attain the 
desired equilibrium temperature, several volumetri 
flasks containing the solutions for a series were 
placed in wire baskets that were suspended in thy 
water bath. The cells were emptied, washed, and 
filled in the room at 25° C, and 30 min or more was 
allowed for the solutions in the cells to attain constant 
temperature. 


4. Spectrophotometric Observations 


The transmittancy readings were made at wav 
length settings 2 my apart over the ultraviolet rang 
of 210 to 400 mau, inclusive. Usually the maxima 
and minima of the bands were first read, then thy 
series of observations recorded at even wa velengths 
from 210 through 400 mu. A check of the values 
at selected wavelengths was then made in reversi 
order. In this manner any change or trend in tly 
transmittancy readings could usually be detected 
A period of 30 to 40 min was required to make 160 
or more observations. The reproducibility of trans 
mittancy observations was usually within 0 
percent. 


5. pH Measurements 


In this study pH measurements were used only a 
a check on the prepared solutions, as hydrogen-io! 
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oncentrations were the known hydrochloric acid 
oncentrations. All solutions were measured with 
ommercial glass-electrode assemblies at 25° C. A 
‘ambridge pH meter was used for most of the meas- 
irements, but Beckman and Leeds and Northrup 
nstruments were also employed. 


{. Molar Absorbancies of Benzophenone, 
4-Aminobenzophenone, 4,4’-Diaminoben- 
zophenone 


Spectrophotometric measurements of various con- 
centrations of benzophenone, 4-aminobenzophenone, 
and 4,4’-diaminobenzophenone were made using 1-, 
2- and 5-em absorption cells. This was done first to 
check the validity of the Beer-Lambert law of ab- 
sorption, which states that at any given wavelength 
or wave number 


ay—A,/bM (1) 


in which dy is the molar absorbancy index, A, is the 
specific absorbancy of the sample (—logy trans- 
mittancy), 6 is the depth of solution in centimeters 
through which the radiant energy passes, and 
is the molar concentration of the absorbing com- 
pound. One may vary the concentration of the 
compound or change the depth of solution by using 
different cell lengths, and if the same curves of molar 
absorbancy index are obtained over the spectral 
range of observed wavelengths, then A,/b.M is a 
constant and there is conformation to the absorption 
law 

The three curves in figure 1 represent the molecu- 
lar, or undissociated, states of benzophenone, R; 
t-aminobenzophenone, RNH,; and 4,4’-diamino- 
benzophenone, R(NH,), in aqueous or alkaline solu- 
tion at 25° C. Molar absorbancy index is plotted as 
a function of wave number with the corresponding 
wavelengths shown at the top of the figure. Benzo- 
phenone has a maximum molar absorbancy index of 
17,300 at 38,760 em™! (258 my). 4-Aminobenzophe- 
none shows a maximum absorbancy index of 15,900 
at 30,100 em~' (332 my) with a secondary band of 
11,900 index at 40,800 em= (245 my). 4,4’-Diamino- 
benzophenone shows greatly increased absorbancy 
with a main band of 23,500 index at 29,600 em™! 
338 mu) and a secondary band of 13,000 index at 
12,400 em~! (236 my). 

The same curve is obtained for benzophenone 
whether in aqueous solution, or acid or alkaline 
medium, as it has no group to dissociate. Also the 
irves are the same for the two amino compounds 
vhether in aqueous solution or in alkali, demon- 
‘trating the fact that in water they do not dissociate 
ippreciably but remain as free bases in the molecu- 
u state. To avoid confusion, the curves for the 
onic states or the dissociated states of the amino- 
enzophenones are not given in this graph. How- 
ver, they are similar to the curve for benzophenone 
ith slightly increased absorbancy and slight shift 
f their bands of maximum absorbancy to higher 
ave numbers 








Fiaure 1, Comparison of the ultraviolet absorption spectra of 
he nzophenone, R, j-aminobe nzophenone, RNHb, and 1 4’- 
diaminobe nzophenone, Ri NH,)» 


IV. Calculation of Dissociaticu Constants 
of 4-Aminobenzophenone from Spectral 
Data 


1. Spectral Absorbancy Data 


A typical set of transmittancy measurements over 
the ultraviolet range of the spectrum for a dissocia- 
tion series is shown in figure 2. They are of 4- 
aminobenzophenone, 5<107~°-M in l-cm cells at 
various hydrogen-ion concentrations at 25° C given 
in table 1. The observed transmittancies in per- 
centage are plotted as a function of wavelength. 
The free base, or undissociated 4-aminobenzo- 
phenone, is represented by curve 17. The com- 
pletely dissociated compound is plotted as curve 1, 
and the 15 intermediate curves represent the partially 
dissociated base at various stages at known hydrogen- 
ion concentrations. It is readily seen that the two 
limiting curves, 1 and 17, intersect each other at 
three points, giving at least two well-defined isos- 
bestic points at 242 my and at 283.5 my, which hold 


S _ 








nA 


FiGurReE 2. Transmittancy curves of 4-aminohe nzophenone as 
the free base (curve 1? 
and 15 intermediate stages during the dissociation 


, the dissociated compound curve I 


’ 
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Volar absorbancy curves of 4-aminobenzophenone 


S804 ated compound (curve 1), 


Figure 3 
as the free base (curve 17), the d 


and 15 intermediate stages during the dissociation 

for the entire series. The constancy of these points 
of equal transmittancy or absorbancy is very im- 
portant in that it indicates purity and stability of the 
compound, and also that only one reaction is taking 
place, in this case, that of dissociation. 

The same data are presented in figure 3, with molar 
absorbancy index 10~* plotted as a function of 
wave number with the corresponding wavelengths 
given at the top of the figure. Curves 1 and 17 
represent the dissociated and free base, respectively. 
As the concentration of hydrogen-ion is increased 
the broad band of the free base at 30,100 em™! 


decreases, until upon complete dissociation there is 


almost negligible absorbancy. One might assume 
that sufficient acid would suppress this band com- 
pletely, but high concentrations of acid were found 
to cause other chemical changes, indicated by some 
changes in absorbancy at higher wave numbers and 
a slight shift in the isosbestic points. The limiting 
curve for the dissociated state was very carefully de- 
termined using 0.3-, 0.4-, 0.5-, 0.6-, 0.81-, 1.0-, and 
2.0-M hydrochloric acid. At 0.6-.M acid, the absorb- 
ancy at wave number 30,100 was lowest and there was 
no evidence of displacement of the isosbestic points. 

The secondary band of the free base, curve 17, 
with maximum absorbancy at 40,800 cm™~', increases 
with increase in hydrogen-ion concentration and 
shifts toward the lower wave numbers as it ap- 
proaches the absorbancy of the dissociated com- 
pound, curve 1. Practically all measurements in 
this portion of the ultraviolet range are on steep 
sides of the bands, as evidenced in figure 2, and do 
not have the accuracy of the measurements at lower 
wave numbers. 


2. Method of Calculation of Concentration Constants 


The acidic dissociation of the ammonium ion of 
t-aminobenzophenone may be expressed as follows 


RNH; = RNH,+H’*. 


The concentration constant, A,, may first be calcu- 


lated, neglecting for the present the activity « 
efficient terms for the free base and ions, according 
the equation 

K [RNH,}[H*] 
IRNH?] 
in which [RNH,], [H*], and [|RNH+#] are the conce 
trations of free base, hydrogen-ion, and dissociated 
compound, respectively. A familiar way of « 
pressing this relation is 


K (1 ~ a) (H*) 
a 

in which @ represents the degree of dissociation into 
ions and 1—qa the remaining free base 

The amounts of either base or ions are propor- 
tional to their molar absorbancy indices as in eq | 
For brevity let a(R), ay(R*), and ay’, represent 
respectively, the molar absorbancy index at any 
W avelength of the free base, the ions, and mixtures of 
the two during the dissociation. 

\s differences in absorbancies are used to caleulat 
a or the amount of RNH} formed at any stage in 
the dissociation, a graphical illustration of th 
differences is given in figure 4. The total change in 
absorbancy from free base to ions is given as th 
outer curve labeled ay(RNH,)- ay(RNH, L Th 
differences for intermediate stages are represented by 
the 12 remaining curves, one of which is labeled 
ay(RNH,) ay’. One can readily see that the 
differences are positive for the large band nearest the 
visible, zero for the isosbestic point at 283 mau, all 
negative between 283 and 242 my, and _ positive 
again for the data at the lower wavelengths. Th 
degree of dissociation, or a, is calculated according to 
eq 4, and the amount of free base remaining, 1 —a 
according to eq 5 as follows 


yg! R) 
(1 yg' R) 


,; 
a vu 


ay(R*) 


a 


ay(R*) 
dy(R*) 


Ayr" 
ay (R) 
Theoretically, the values of @ should be the sam 
across the spectral range regardless of wavelength 
It is obvious that the best estimates will 
where the differences in absorbancies 
are greatest. Data too ne«r the isosbestic points ar 
worthless in such calculations. It is wise to examin 
the a values across reasonable wavelength ranges in 
more than one band, to detect a possible trend. Such 
trends have often been noted [1, 2], particular!) 
where impurities may be involved, or where anothe: 
reaction can conceivably take place simultaneous!) 
Dissociation constants calculated from absorbanci 
at one or two wavelengths are not necessari 
reliable. In this case the authors selected the way 
length range of 314 to 350 mu, inclusive, to calculat 
an average a from 19 wavelengths. A comparis 
was made with values calculated from 256 to 2' 
my in which greater variations occurred becau 


selected. 
be obtained 
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bservations were made on steep slopes of bands, 
nd there was satisfactory agreement. 

The ratios of free base to ions may be calculated 
lirectly by combining eq 4 and 5, thus 


a dy ay(R 
ay(R) dy 


his ratio multiplied by the hydrogen-ion concentra- 
on gives the concentration constant A,, as in eq 
or 3. p,, the negative of the common logarithm 
I ean then be calculated. 


3. Presentation of Data 


'ypical absorbancy data and calculated values for 
the dissociation of 4-aminobenzophenone at 25° C 
ure given in table * Two series show reproduc- 
bility of results from two different runs, using 
freshly prepared solutions throughout. The molar 
concentration of hydrochloric acid is given in column 
|. Molar absorbancy index at 332 muy is listed in Figure 4. Differences in absorbancy of undissociated 4-amino- 


column 2 Column 3 shows a calculated from he nzophenone and the partially and completely dissociated 








compo ind 


absorbancies at one wavelength (332 mu), and 
column 4 gives an average @ calculated from ab- 
sorbancies at 19 wavelengths from 314 to 350 mu, average a, however, should minimize the error 
inclusive. Inasmuch as the transmittancy observa- Values of a below 0.1 and above 0.9 yield unusually 
tions are subject to error in the third figure, this large or small ratios of free base to ions and were 
error is carried in the subsequent calculations. An not included in the calculations of the constants 
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TABLE 1. Dissociation of 4-aminoben zo phenone at 25° ¢ 


Amount dis Amount dis 

sociated into | sociated into Ratio of free 

ions, A332 ions, A314 to base to ions 

1 wavelengtt 350, incl 

Average 
19A 


Molar con Molar ab 
centration sorbancy 


of HCI index, A332 


H 


0. 0000 
0001 
O065 
aol 
Oo1ls 
Oo26 
Oo4l 
0052 
Ove 
o102 


O61 
0256 
0405 
1020 


6440 





The concentration constant, A,X 10°, is given in 
column 6, and —logy» A,., or p&K,, in column 7 
The increase in pA, with increase in ionic strength 
is slight, but nevertheless there is a trend. The 
arithmetic mean of the values of pA, for each series 
is 2.18, corresponding to 0.0066 moles/liter for K.. 

A summary of calculations from spectrophoto- 
metric data that was taken of the dissociation of 4- 
aminobenzophenone at 10°, 15°, 20°, 30 and 
40° © is given in table 2. Each series of solutions 
for the dissociation at one temperature was com- 
pleted before changing the temperature of the water 
bath for the next Practically the same 
amounts of b vdrochloric acid were used for each 
series to facilitate direct comparison of the amount 


or 
oo 


series 


on of 4-an 


Papi ssocial 


It is readily seen from the table that for any given 


hvdrogen-ion concentration the degree of dissocia- 


temperature. The 
and pk, 


with 
constant 


tion decreases increase In 
concentration 
correspondingly decreases 

Of interest is the magnitude of error that any 
Variation in temperature may have upon the absorb- 
ancy of a given solution. As the limiting curves, all 
free base or all ions, were found to be substantially 
unchanged over the temperature range studied, any 
change in absorbancy can be attributed solely to the 


thereby increases, 


change in constant of the dissociation reaction with 
change in temperature. A solution of 4-aminobenzo- 
phenone in 0.0064-\7 hydrochloric acid, which at 
25° C is practically half base and half ions, was placed 
in the cells, and the temperature of the bath was 
changed in the following manner. An equilibrium 


inobenzophenone at 


of free base dissociated into ions at a particular h 
drogen-ion concentration at different temperatur: 
The spectrophotometric observations were made 
usual over the ultraviolet spectral range of 210 
400 my, to assure accuracy of data and adherence 
isosbestic points, but the degree of dissociation 
the average calculated from the 19 wavelengt 
between 314 and 350 mu. It was feared that erra 


results might be obtained from the measurement 


at 10° C, as a slight film formed on the outer fa 
of the quartz end-plates as they were removed fri 
the cell compartment at 10° to the room at 25 
for filling and emptying. However, the results 
not seem to be out of line. At 15° C the film y 
not visible to the naked eye 


10°, 15 


temperature of 25° C was first attained, and obse 


( 


tions were made at 2-mu wavelength settings from 


314 to 350 mg, inclusive. The bath was then chang 
to give equilibrium temperatures in 5-deg st 
first up to 40° C, then in reverse order down to 10° ¢ 
and again up to 25° C The cells w 
not removed from the cell compartment, and a! 
3 days was required for the experiment. The chat 
in absorbancy were determined and the stabilit 
the solution also proved The results are show 
figure 5, in which specific absorbancy of the solu 
is plotted as a function of wavelength. In this 
ticular case there is roughly an 8-percent increas 
absorbancy with a 5-deg increase in temperat 
It is evident that for accurate results close temp 
ture control is necessary for all intermediate st 


as a check. 


in the dissociation. 
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V. Other Related Thermodynamic Quan- 
tities 
1. Estimation of pA, and K, from 10° to 40° C 


True thermodynamic constants, A,, may be cal- 
culated at each temperature studied by eq 7 


RNA 


me /RNH 

which the activity coefficients, /, are expressed 
with the appropriate subscripts 

When pA, is plotted as a function of ionic strength 

each temperature, it is found that straight lines 
may be drawn through the values, giving slight 
slopes. Extrapolation to infinite dilution should 
vive the pA, values, and the activ ity coefficient terms 


TaBLe 3. Activity coefficient terms, fr 


IRNHy JRNH 


may then be expressed as a linear function of ionic 
strength. The results are shown in figure 6. The 
experimental error is at least +0.01 pA unit, and 
the lines were drawn as visual estimates. pA, was 
found to be 2.36, 2.28, 2.24, 2.17, 2.12, 2.06, and 2.02 
for 10°, 15°, 20°, 25°, 30°, 35°, and 40° C. K, is 
then 0.0044, 0.0052, 0.0058, 0.0068, 0.0076, 0.0087, 
and 0.0096 in the same order. 


2. Activity Coefficient Terms 
The straight lines in figure 6 may be used to select 
pk, values at any ionic strength between 0 and 0.03 
Corresponding A, values may then be derived. As 
the ratio of K, to AK, is equal to the activity coefficient 
term fut+/rxu,//exn?, this term, called fe, may now 


be calculated. Some terms are shown in table 3 for a 
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j-aminobenzophenone 


Figure 6 pA 


stages in dis ral ’ hudrochlo 


or fa, cal ilated for several te mperatures 








few ionic strengths at the seven temperatures at 
which measurements were made. 

There are little data in the literature on activity 
coefficient terms of weak bases at low ionic strengths 
However, if the terms at high ionic strengths are 
extrapolated to lower strengths, the results of 
Randall and Failey |4], Harned and Robinson [5], 
and Harned and Mannweiler |6] obtained by electro- 
motive-force methods and Weil and Morris |7], who 
used spectrophotometric methods, are in agreement 
with our results. 


3. Heat of Dissociation 


lf one now plots the log AK, values for each 
temperature as a function of 1/7, a practically 
straight line is obtained as shown in figure 7. Thus 
the heat of dissociation can be assumed to be prac- 
tically constant, and AH may be calculated according 
to.eq Ss 


d log K/dT=AH/2.3 RT?, (8) 


in which AH represents the molar changes of heat 
content for the dissociation, R is the gas constant 
8.3144 joules deg™' mole', and 7 is the temper- 
ature in degrees Kelvin (degrees Celsius + 273.16). 
AH is found to be 19,000 joules deg! mole. 


4. Basic Ionization Constant, K, 


One may derive the basic ionization constant, K5, 
from the usual relationship 


K,—A,/K,, 


in which K, is the activity product for water. For 
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Figure 7 pA, alues at temperatures from 10° to 40 
plotted asa Junction of l 7 


example, at 25° C Ky, is equal to 107-'*/(6.6 10 


which yields a value of 1.610°". p&p, is then 


approximately 11.8. 
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A New Method of Radioactive Standard Calibration 


By Howard H. Seliger 


By proper arrangement of experimental observations and statistical analysis it has 
been possible to make Geiger-counter measurements of different activity samples inde- 


pendent of intermittent disturbing effects 
involving more than 1.4 


The results of 5,328 individual measurements, 
10° counts made by using a Latin square arrangement, are shown 


to have an error distribution identical with that expected from the statistical nature of the 
disintegration process A completely worked out example of a calibration procedure is 
given, in which it has been possible to achieve a higher degree of accuracy in only 30 percent 


of the original counting time 


I. Introduction 


Anyone who has had occasion to make Geiger- 
counter Measurements on radioactive isotopes has 
found at times that his results were not self-con- 
sistent; or more precisely, that readings taken at 
time A could not be duplicated within the expected 
Poisson deviation at time B. Sometimes upon fur- 
ther investigation it was found that these variations 
were due to fluctuations in the electronic voltage 
supply, slight changes in sensitivity of the Geiger 


counter itself, or other transient effects. M 
often, if it were not too time-consuming, the orig 
data were discarded and the measurements repeat 
this time with the fingers crossed. In this pap 
procedure for arranging and analyzing a series 
measurements is deseribed that eliminates u 
mittent disturbing effects and that permits a de! 
mination to be made as to whether the variance 
the set of readings is a reasonable one to exp 
on the basis of a Poisson distribution of cou 
The end result is that all readings of a set are trea 
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s though they were taken at the same time with 
he same equipment. 


the Chi Square 


(Il. Part 1 of the Solution 
Test 


if we know that radioactive disintegrations follow 

Poisson law, we can predict the variation within 
vhich individual values should lie. Our problem 
s to determine the probability that the variation of 
in observed set of readings is a reasonable one to 
xpect if the variations are due only to the random- 
ess of the disintegration process. This probability 
distribution of the ratio of the observed variation to 
the expected variation has been calculated and is 
known as the Chi square distribution.' The Chi 
square statistic is defined as the sum of squares of 
the deviations divided by the mean, or 


(1) 


where bi Is the observed counts per reading, A is 
the expected, or average, counts per reading and 


denotes the summation over all the readings. 


Then = (X,—A)? is the sum of the squares of the 
' 

deviations from the mean count, and therefore an 
estimate of no* calculated from the individual counts. 
When the measurements do follow a Poisson law, the 
average count is known to be an estimate of o* for 
the measurements. Thus the expected value of this 
ratio (eq 1) when the data follow Poisson’s law is n. 
This is, however, subject to statistical variations. 
Tables are available showing for various n the 
probability limits for this ratio (see footnote 1). 

Let us now apply this simple test to a set of con- 
secutive readings taken with a Geiger counter and 
scaler, with a constant activity source. All readings 
are taken for a constant time interval, which does 
not enter into the calculations. 


Total 
counts 


Reading 
number 


Total 


counts 


Reading 
number 


14, 805 
15, 378 
15, 021 
15, 040 
15, 335 


$35,622 


‘ und VanVoorhies, Statistical procedures and their mathematical 

p. 404 (McGraw Hill Book Co., Inc., New York, N. ¥ 1940): Fisher, 
stistical methods for research workers, p. 58 (Oliver and Boyd, London, 1944 
A. Brownlee, Industrial experimentation, 3d °d., p. 39 (British Ministry of 
ply, Directorate of Royal Ordnance Factories 's 


As we have 10 readings, there are 9 degreesof freedom. 
If we look up a table of x’ for n=9, we find that the 
probability of getting a ratio of 28.7 is less than one 
chance in a thousand.? Therefore, we conclude that 
the above set of readings do not follow a Poisson 
distribution, and that data taken with the counter 
and scaler are not reliable. Conversely, if this ratio 
came out within the values 14.68 and 4.17, which are 
the 10-percent and 90-percent probability limits, 
respect ively, we could conclude sat isfactory operat ion 
on the part of our equipment and continue with our 
measurements. Thus even this simple application of 
the Chi square test is an extremely useful one for 
checking equipment, and in the following sections 
we shall show how it and the method of Latin squares 
are used together for more complicated data analysis 


III. Part 2 of the Solution 


The Latin square was originally developed for agri- 
cultural field trials. If A different strains of corn 
are to be compared, the experimental area is sub- 
divided into a checkerboard of A rows and K col- 
umns, making available A? plots. The plantings are 
applied so that each strain appears once in each row 
and in each column of the A? plots. This requires 
that there be A plots assigned to each strain of corn 
Every row and every column of the area contains a 
complete set of the A’ strains under test. Intuitively 
it may be seen that such an arrangement will be ef- 
fective in reducing the error of comparisons among 
the strains since if any one of the columns or rows is 
more fertile than the others, all A strains obtain this 
advantage and the relative performance of thé sev- 
eral strains is unaffected. If there exist fertility 
gradients in the area that would ordinarily interiere 
with the accuracy of the comparisons, the Latin 
square method of allocating the plantings to the 
plots has ensured a more nearly equal opportunity 
for all strains. 

A typical 4 by 4 Latin square arrangement is shown 
in figure 1, a. The four plants of any given strain 
appear in the four rows and four columns. The four 
plots of any one strain may in consequence show 
considerable variation among themselves, and should 
not be used directly as an estimate of the error in the 
comparison of the different strains. However, the 
comparison of the averages for the strains has been 
made more accurate by insuring that all strains 
sample the various strips of fertility. Thus differ- 
ences among rows and columns no longer contribute 
to the error of the comparison. It is necessary then 
to find another way of estimating the error of the 
average for each strain, as the customary method of 
calculating the dispersion from plots receiving the 
same strain is no longer applicable. However, be- 
fore we proceed to outline the new method for cal- 
culating the dispersion, let us carry on with the 
development of the Latin square arrangement. In 
figure 1, b, is a Latin square arrangement for the 
Greek letters a, 8, y, 6. Now let us superpose figure 
1, b, upon figure 1, a to obtain figure 1, ¢. If the 


Brownlee, Industrial experimentation, 3d ed., p. 161 (British Mir 
torate of Royal Ordnance Factories, 1948 


the Latin Square 


K. A 
f Supply, Direc 
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(a) 
Aa} BB 
C8 | DY 
DB | Ca 
BY | A& 
Ficure 1. General Lat 
s, Latin square for strains A, B, C, D; b, Latin square with Greek lett 


four alpha plots in figure 1, ¢, are examined, it will 
be seen that these four plots contain a representative 
plot from every row and every columnand every strain, 
the strains still being the roman letters A, B, C, D 
The same is true for the beta, gamma and delta 
plots. The averages for these Greek sets therefore 
should show agreement, and the variation 
among these four averages does in fact provide the 
estimate of error relevant to the comparison of the 
averages for the four strains. In the same manner, 
Arabic numerals can be assigned to the plots, so that 
the four plots marked 1 in figure 2 contain represent- 
atives from every row, column, strain, and Greek 
letter. This is also true for the plots marked 2, 3, 
and 4. These four sets also should agree and can be 
used as a further estimate of the experimental error 
The Latin square shown in figure 2 is now unique, or 
completely orthogonal, in that this is the only way 
in Which the individual simple Latin squares composed 
of Laiin letters, Greek letters, and Arabic numerals 
as shown in figure 1, can be superposed subject to 
the restricting conditions above We would expect, 
therefore, that there should be associated with these 
classifications as many degrees of freedom as there 
would be if only one number or letter appeared in 
each of the 16 squares. 

It is customary to summarize the statistical compu- 
tations for experimental arrangements of this type in 


cood 
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(c) 


‘ 


r 


(b) 
CY | 08 
AB | Ba 
BS | AY 
Da |} cB 
n square arrangement 
uperposition of b upon a, subject to the same restricting conditions 


a tabulation of the analysis of variance of the five 
classifications that have been enumerated. 


I] Jegrees ol 
freedon 


Variation among the four row averages 3 
Variation among the four column averages 3 
Variation among the four strain (Latin 

letter) averages 3 
Variation among the four Greek letter 

averages 5) 
Variation among the four Arabie number 

averages 3 


They are 


As seen with each of these is associated three degrees 


the usual (n 


(corresponding to 
deviations). The 


computing standard 


of freedom 
divisor in 


total of 15 degrees of freedom is that properly as- 


signed to 16 observations if a standard deviati 
were computed for the 16 results ignoring all class 
fications among the plots. Now, not only do the 
degrees of freedom add up to the correct total, but th 
sum of the squares of the deviations within ea 
classification, when added up over the 5 classifi 
tions will be found equal to one-fourth of the sum 
the squares of the deviations of the 16 observatio 

















C28 
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D3B| Cea Azy 




















Ar rangement of Latin square to indicate sequence 
of readings. 


about the grand average. (The factor one-fourth 
arises because the five components have been com- 
puted from averages of four, in contrast to the use of 
individual observations). It is therefore possible to 
partition the over-all variance into these five por- 
tions. Two of these are the rows and the columns, 
which by experimental arrangement, have been 
removed from participation in the experimental com- 
parisons. One portion (Latin letters) directly re- 
fleets differences among the strains and the remain- 
ing two portions (Greek letters and Arabic numbers), 
vield independent estimates of the experimental 
error. Completely worked out examples of this 
form of Latin square will be found in Brownlee (see 
footnote 1). The net result is that the computation 
leads to determining the ratio of the mean sum of 
squares for each classification to the mean sum of 
squares for error, as defined by Brownlee (see foot- 
note 1). This ratio, called F, has the expected value 
of unity in all cases where the classification under 
examination contributes nothing to the variation of 
the results, as would have been the case for instance 
if only one strain of corn had been sown in all 16 
plots. The computed value of F will in such a case 
be expected to fall in the neighborhood of unity, the 
reasonable limits of variation being determined by 
the number of degrees of freedom attached to the 
numerator and denominator of the ratio. The 
theoretical distribution of the F ratio has been 
tabulated (see footnote 1) and is an indispensable 
guide in passing judgment on a particular value of F. 
These tabulations are listed as Tables of Variance 
Ratio for various significance levels, corresponding 
to the probability levels in the x? test in section I. 


IV. Chi Square Test Modified for the Latin 
Square 


In section II the principle of the Chi square dis- 
tribution was demonstrated, and in section IIT the 
e of the Latin square for the separation of dis- 
rbing factors from the error calculations and for 
e calculation of variance ratios was outlined. The 
xt steps in our sequence are: (1) To apply the 
itin square to our particular problem of calibrating 


unknown sources of different activities against a 
known standard source; (2) to modify the calcula- 
tion of the Chi square test to fit the data taken in 
our Latin square arrangement. 

Perhaps the application of the Latin square to our 
particular problem will be easier to follow if we start 
with the end result and then discuss the logic of the 
arrangement as given. Consider again the Lativ 
square in figure 2. The 16 subdivisions of the large 
square are numbered consecutively from 1 to 16 in 
columnar order in the upper left-hand corners. 
These numbers represent the sequence of the readings 
taken. The Latin letters A, B, C, and D represent 
a primary standard source and three unknown 
sources, respectively. We shall consider the stand- 
ard source and the three unknown sources as our 
“strains”. As the observed counting rates due to 
the sources A, B, C, and D are in general different, 
(due to the fact that the absolute disintegration 
rates of the sources are different) a variance among 
the four source averages (Latin letters) is to be ex- 
pecced and therefore should not be included in any 
error calculations. However, the sums of sets of 
four readings as for instance the sum of the a’s, 
include only one reading each of A, B, C, and D, 
and therefore the sum of the a’s theoretically should 
be equal to the sum of the 8’s, y’s, and 6’s. This 
being the case, we can perform our Chi square cal- 
culation on the sums of four readings instead of on 
the individual readings. Again, from section III, 
we have the following sets of sums of four readings: 
Rows; columns; like Latin letters (sources); like 
Greek letters; like Arabic numbers. 

The variation among the sources (Latin letters) is 
immediately ruled out for the reason already given 
From the sequence of the measurements as given by 
the numbers in the upper left-hand corners of the 
squares in figure 2, it is seen that the sums of the 
columns will be especially sensitive to short-term 
intermittent disturbing factors of all kinds. Thus 
the variance among the four column averages should 
not be included in any error calculation. More- 
over, these intermittent fluctuations, if they occur, 
are analogous to the fertility differences in our agri- 
cultural example, so that the relative performance 
of our four sources should be unaffected (if we con- 
sider sums and not individuals). 

Now we have remaining three sets of four sums 
each for the calculation of error. These are the rows, 
the Greek letters and the Arabic numbers. The 
next step is to modify our Chi square calculation to 
fit these data. 

We know that in our particular case of radioactive 
disintegrations o*, the standard deviation squared of 
our Poisson distribution, is equal to V, where N is the 
total counts involved. The estimate of no? can be cal- 
culated from the three sets of four sums. In this case 
there are 3 degrees of freedom from each set, making 
a total of 9 degrees of freedom. 

Table 1 shows a complete set of data taken in the 
calibration of three unknown sources B, C, and D 
against a laboratory standard source, A. — In order 
to indwate clearly the sequence of the readings and 
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TaB.Le 1. Data as taken in an actual calibration of sources. 

MNO. [SAMPLE|oUNTS| ‘SEC GROSS TET 
a |34480| 1277 | 2700 | 2646 
2 s | 27364] 9731 | 2812 | 2758 

3 c | 25650| 864 | 2969 | 29:15 

4 0 | 25774| 9568 | 3008 | 2954 
5 s | 27060] 947: | 2857 | 2803 
| 6 | A |30540| 1109 | 2754 | 2700 
7 o | 27478! 914 | 3006 | 2952 

8 c | 25628] 850 | 3015 | 2961 

9 c |256i4| ese | 2985 | 2931 
10 | D  |26404] 8781 | 3007 | 2953 
T a | 26082] 9654] 2702 | 2648 

i2 | B | 25612 | 9032] 2836 | 2782 

3 | D | 25650] 864 | 2969 | 2915 

4 | c | ese96| see2| 2967 | 2913 
is | B | 25654] 902.1 | 2844 | 2790 

i6 | A | 26510] 980 | 2705 | 2651 


























the*method of summation, these data have been 
arranged in the form of a Latin square in table 2. 

The analysis is performed in the following manner: 

1. Caleulate the average row sum of the four sets 
of row sums. Notice here that in order to do this, 
all 16 readings are included, and this grand sum is 
divided by 4. Therefore the average row sum 1s 
equal to the average Greek letter sum and also to 
the average Arabic number sum. Let us call this 
value 


=(4 113.18. 


2. Calculate the sum of squares of the deviations 
of the individual row sums from the average sum of 
four 
(0.23*+-0.067 +-0.13*+-0.30°=0.1634). 
3. Calculate the sum of squares of the deviations 
of the individual Greek letter sums from the average 
sum of four. 


(0.25? +0.042+0.572+-0.367=0.5186). 

4. Calculate the sum of squares of the deviations 
of the individual Arabic numbers sums from the 
average sum of four. 

0.32? +- 0? + 0.267 + 0.06" = 0.1736 


Oa 


TABLE 2 


Data from table 1 arranged in a Latin squa 
make the summation procedure easier to follow. 





ROW 
TOTALS 





! 5 9 13 
26.46 | 2803 | 293i 29.15 112.95 
Aia_| B26 | c3y | 048 
2 6 10 14 

2758 | 2700 | 2953 | 2913 | 113.24 


Bay A35 D2a ciB 











3 7 T 15 
29.15 | 2952 | 2648 | 2790] 11305 
C25 DIY AsB B3a 

4 8 12 16 | 
29.54 | 296! | 27.82 | 265) | 113.48 | 


038 | Caa Bis A2y 


COLUMN 
TOTALS W273 | 4.16 1314 | 2.69 








ARABIC |! + 7+12+14]3+5+10+16/4+6+9+415 |2+8411+13 





TOTALS] 11293 | 13.22 | 3.75 | 1282 





GREEK |!+8+10+15 |4+5+11+14/2+7+9+16/3+6+12Ha 





TOTALS} 113.50 | 113.18 | 11292 | 3.12 








LATIN 
OR A B C D 
SAMPLE 
TOTALS | 10645 | lii.33 | 117.20 | 117.74 























5. Sum these three sums of squares. 


(0.1634 +-0.5186 +-0.1736—0.8556). 

6. We must digress for a moment and consider 
now the basis for this calculation. The expected 
standard deviation squared of a number of events 
that follow a Poisson law is N, when N is the total 
number of events. However, as is seen in table 1, 
our sample deviations are calculated from observed 
counting rates. In order that our ratio of observed 
to expected deviations be a unitless quantity, we 
must express the sum of squares of the deviations as 
total counts. This can be done in the following 
manner: 

In each individual measurement we observe 
least 25,600 counts, which is 200 message registe! 
counts on a seale of 128. Therefore, we can, with 
negligible error, assume that for a set of four read- 
ings the total number of counts involved is 
4 25,600=102,400. If we divide this value 
the average counting rate for a sum of four read- 
ings, we get a value with the units of secon 
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rresponding to the average time per reading. 
his value squared will, when multiplied by the sum 
squares of the deviations that has the units of 

counts)/(seconds)]*, give units of (counts)*, the 

ume as our denominator. 

Our calculated ratio is now 


__.. 102,400? 
0.8900 X< 
: kav 


6.8 
102,400 84 


In our table of x’ we find, for 9 degrees of freedom 


X"70 6.39, 


Thus our ratio has greater than a 50-percent prob- 
ability of occurence on the basis of a Poisson law, 
and we conclude that the resulting calibrations are 
perfectly acceptable. 

As an illustration of the efficiency of this method 
in eliminating transient disturbances from the meas- 
irements, let us perform an F test on the column 
totals. The mean square column deviations are 
1.4046/3, as here we have 3 degrees of freedom. 
The mean-square-error deviations are 0.8556/9, 
since there we used 9 degrees of freedom. Thus the 
ratio 4.92, as we can see from tables of Variance 
Ratio has reached the 5-pereent significance level. 
Evidently in this case the Latin square arrangement 
was successful in removing this effect from the vari- 
ance of the sample averages. It is the elimination 
of these disturbances that accounts for the increased 
accuracy over the old direet-comparison method 
and the subsequent 70 percent shorter counting 
time. This reduction in counting time is due to 
two factors. First, the number of runs used was 
reduced. Where previously a standard source and 
three unknown sources were run alternately, so that 
for 4 runs on each of the three unknowns the stand- 
ard would be run 13 times, making a total of 25 runs 
for the entire set of readings, we now have only 16 
runs in which the standard is measured as many 
times as the unknowns. Second, as the extraneous 
nonPoisson) variation was present in the former 
method, a much larger counting time was necessary 
than that based on a Poisson distribution in order to 
obtain the precision called for. The Latin square 
arrangement permits the counting time to be cut in 
half, and the precision of the measurements is still 
greater than in the old method. 

The standard deviation of any of the sums of four 


0.8556 3 
8; 308. 
> 9 


The standard deviation of any individual count- 


¢ rate is 
0.8556 3 
[ 0. ee | 0.154. 
4x9 


For example, the average counting rates of sources 
A, B, C, and D, which are the objectives of this 
entire procedure, are, respectively, 26.61, 27.83, 
29.30, and 29.44+.15 c/s. These results can also 
be calculated as shown in Brownlee (see footnote 1). 
However, the method as demonstrated is more 
direct and easier to calculate. 


V. Summary 


In this paper two things have been accomplished 
First, by the arrangement of the experimental data 
in the form of the orthogonal Latin square, thereby 
eliminating intermittent effects, the accuracy has 
been increased to such an extent that only 30 percent 
of the original counting time is required. Second 
it is now possible to determine the actual exper!- 
mental standard deviation of any particular average 
reading, and in the process of obtaining this standard 
deviation it is also possible to determine whether 
or not the readings follow a Poisson law, which is 
the case when the equipment is Operating satis- 
factorily. 

As can be seen, the use of Latin squares for the 
statistical arrangement of measurements is an 
extremely powerful tool, by means of which a 
number of separate effects may be determined or 
eliminated. It has proved to be extremely efficient 
in the laboratory, not only in the actual calibration 
of standards, but in obtaining useful information 
relative to the long-range operating characteristics 
of the Geiger counters and associated electronic 
equipment. The general method as outlined should 
be applicable to a large variety of other types of 
measurements, especially where time effects or 
other disturbing factors are a source of trouble 

The results of 333 such sets of data as shown in 
tables 1 and 2 have been tabulated and are presented 
in figure 3, together with their theoretical distribu- 
tion. As can be seen from the figure, the experi- 
mental distribution closely approximates the theo- 
retical one. There is one more point that can be 
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made only after the accumulation of such a large and theoretical curves illustrates one further po 

amount of data as this. The small peaks at the that the Latin square arrangement has been ¢ 

right of figure 3 indicate a slight tendency for the pletely successful in eliminating extraneous, ; 

observed counts to be influenced by effects other Poisson fluctuations from the measurements. 

than statistical fluctuations. However, it must be 

remembered that the data include instances where 

either counters or scaling circuits have gone bad 

and have evidenced a need for repair or replacement, The author thanks W. J. Youden and J 

so that this slight excess is to be expected. It is Cameron, of the Statistical Engineering Section 

believed that this is the first time that such a large their suggestions as to the use of the Latin sq ui 

collection of Chi square values has been obtained and for many valuable hours of discussion: 

from the ratio of the error variance of Latin squares Margaret Selgin and Lucy Cavallo, of the Radio- 

to the expected (Poisson) variance activity Section, who performed the calculations 
It is well known that individual measurements of summarized in figure 3 

radioactive disintegrations follow a Poisson law, 

but the close agreement between the experimental WasnHincton, July 26, 1950 





A Problem in Precision Cam Design 
By Joseph Blum 


rhis paper proposes an analytic method for determining the profile of a cam required it 
a device that demands high precision. The method is applicable whenever the equation for 
the profile can be obtained by making a less stringent assumption, name lv, that the follower 
makes contact at a fixed point of its extremity The true profile can then be determined 
1umerically by considering it as the envelope of a family of curves 


The design of an X-ray spectrograph often requires In figure 1, a, O represents the axis about whic! 
that a constant rate of oscillation be maintained in the cam rotates, A the axis about which the crysta 
the crystal holder. The crystal rotates through an holder rotates, and # the center of the ball at 


angle 8B (the angle of travel), reverses its direction end of the follower. The follower ACB makes 
and then proceeds with the same angular velocity. contact with the cam and serves to rotate the crysta 
This type of motion may be obtained from a cam holder. The distances AO and AB are equal; denot 
rotating with constant angular velocity, #. An _ this common length by R. Let 7; and r: be, respe 
analytic method for computing the profile of the tively, the least and greatest radii measured fron 
cam is presented in this paper which extends, in a 0 to the point B; let r be the variable radius, OP 
certain sense, the application of a paper written by corresponding to the angle a, which r makes with th: 
J. B. Friauf;' the notation therein will be preserved least radius r;. 

to a considerable extent Denote the time by ¢ and let t=0 represent th: 
time when the least radius coincides with the lin 
segment extending from 0 to A. Then the angi 
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Ficgure 1, 
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tween 7, and OA is given by of. Let @ denote 
angle OAB and @# the angle AOB. Unless otherwise 
ndicated, all angles are in radians. 

From figure 1, a, it is clear that 


6 


. r 
sin 
‘ ‘ , 
2 2h 


angle 6, 
when 


smallest 
angle @, 


that the 
the largest 


follows 
and 


mm which it 
when 


ro, or 


eurs / / 
A 


. r 
— Se 


It is also clear that 6 Differentiating eq | 


vields 


4dée 1 di 1 dr da 
R dt R da dt 


COs 


2 dt 


da 
%) “dt 


snce a wt 
Substituting this into eq 4 yields 
do) 
dt* 


dedi 


4déa 
CUS 


2 dt 


1 d? 
Rda | a 


| 
(gain, since ¢=,5(4—9), 


, and sub- 


do/dt 
stituting this into eq 5 results in 


4dée ldr 
COs 


1d@ 
2dt Rda (+ ) 


2 dt 


The follower rotates through the angle p during the 
same time interval that it takes the cam to rotate 
Therefore, as the respective 


through 7 radians 


angular velocities are constant, 


dé 
dt 8B d@ Bw 


’ r 
Ww T - dt T 


Substituting this into eq 6 gives 


0 1 d? T l 
2 Rda\gt2) 


cos (6/2 1/2R)y4h?— Substituting 


From eq I 
this into eq 7 and simplifying vields 
2dr da 
\ t/}?? a B 


where B,=71/8 
\fter integrating this equation there is obtained 


a 


2 sin + 
B; 


Cc"; 


1 r 
aR 


where C is the constant of integration. 


When a=0, r=r,; solving for C, results in C,=2 


sin ~'r,/2R=—@,. Equation 9 can be expressed in the 
form 
lia 


2P sin +9, >> 
> 2 (8B, \ 


10) 
The greatest radius, 7 when a=2+ 8/2, 
which may be vertified by substituting into eq 10 

If the ball that makes contact with the cam were 
sufficiently small, then eq 10 would give one part of 
the cam’s contour, and a similar equation would give 
the remaining part. Frequently, however, this con- 
dition is not met, and the method given in the ref- 
erence in footnote 1 must be extended 
high level of precision. Consider then the family of 
circles of radius Pp, whose centers lie on the curve 
given by eq 10 The envelope of this family con- 
sists of two distinct curves; the one that lies nearer 
the center of curvature of eq 10 coincides with the 
contour of the cam. This envelope can be obtained 
by a numerical determination of the C-discriminant 

Let the curve of eq 10 be given in Cartesian co- 
ordinates by 


occurs 


to secure a 


1] 


Let (2,y) be an arbitrary point on this curve. The 
equation of the circle with radius p and center (2,y 
is given by 


p. 


Substituting y=/(2) into eq 12, there is obtained a 


one-parameter family of circles 
X—rz}?-+ ; p 
Differentiating eq 
0 14 


The C-discriminant is obtained by solving eq 13 and 
14 simultaneously for a discrete set of values of the 
parameter z. The method outlined below vields the 
numerical solution: 

1. For a selected set of values a;; determine the 
values of 7 by using eq 10 Denote these points 
by P,. 

2. Let 7; be the angle of inclination of the tangent 
line at P,; let y,; be the angle between the radius r, 
and the tangent line. Then, since 7,;=a;+ y¥,;, 


+tany 
tan a, tan y, 


tana 


' where r/ is the derivative dr/da 
This leads to 


Now tan y=) 
evaluated at P,. 


l/e 
tan y,; = 26, tan =( ; +4 ): 


?E. L. Ince, Ordinary differential equations, p. 85 to 
& Co., New York, N. Y., 1944). 
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These relations enable one to compute the value of 
f’ (x). 

3. The solution of eq 13 and 14 is easily obtained 
in the form 


Y pf'(a 
. r 
\! f(a 
(15 
) y¥it e. 
vil+it (ry) 


The points (Y,, Y,) are on the envelope, as required. 


To obtain the other part of the cam, consider that 
phase of the motion in which the angle @ decreases 
(fig. 1, b For this phase, 


dé Bw 
dt T 


Substituting this into eq 6 results in 
6 ldr/fl_ j 
cos . ( ) (16) 
2 hi dy 2 8 
where y denotes the angle that r makes with the 
greatest radius r,(see fig 1,b). Substituting cos (@/2) 


(1/2R) ,4k?—r* into eq (16) and simplifying gives 


2dr dy ~ 
’ (17) 
V 4/2? 7 B» 
where B,=1/2—-x/B 
When eq 17 is integrated there results 
2 sin-' — v4¢ (18) 
2k 8B 


where C, is the constant of integration. When y=0, 
r=?r,; solving for C, yields C,—2 sin~'r,/2R=8, 
Equation 18 can be expressed in the form 


+0, ): 


Asa check, when y= 2—8/2,r=2R sin }4 (@.— B) 
2R sin 4% @=—r,. The envelope of the family of 


l “ 
r=2R sin = ( (19) 








ABER 


\ 








Figure 2 


circles with centers on the curve (19) and radius p 
leads to the remainder of the contour of the cam. 

To complete the analysis, consider figure 2. The 
are ABC is the locus of eq 10, and are A’RB’C’ is 
its envelope. Similarly, the are DEF is the locus of 
eq 19; arc D’ RE’ F’ is its envelope and may be found 
numerically by following the procedure outlined 
above for eq 10. The high point of the cam is 
located numerically by solving for the intersection 
of the two envelopes. The low portion of the cam 
does not coincide with the envelope, as is evident 
from the figure. Let 7 be the point of intersection 
of the curves 10 and 19. The circular arc, PQ, with 
center at 7 and radius p constitutes the low portion 
of the cam. These considerations enable one to 
complete the numerical solution of the cam design. 

In summary, it should be noted that the essential 
idea in this approach is to first find the equation of 
the contour, assuming that the follower makes con- 
tact at a fixed point of its extremity; then, by using 
the true form of the follower, obtain the actual cam 
contour by determining the envelope of a family of 
curves. It is evident that the analytic method ap- 
plied in this particular case may be used in other 
problems of precise cam design. 


WasHiInGtTon, August 4, 1950. 
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Real Roots of Real Dirichlet L-Series ' 


By J. Barkley Rosser 


In the theory of the distribution of primes in arithmetic series, in assigning bounds in 
the three prime theorem, and in studying the class number of quadratic fields, a knowledge 


of the location of the real zeros of L(s, x) is of value. 
If proved, this result would be of value in each of the 
By a certain computational procedure the conjecture has already 
In the present paper, this earlier computational 


are no positive real zeros for any k. 
fields mentioned above. 
been verified for each individual k<67. 


procedure was tried for each k<227 and failed for k 
procedure is given in the present paper, but still the case k= 163 remained difficult. 


A long standing conjecture is that there 


163. An improved computational 


Finally, 


a new formula for L(s, x) was discovered that made it possible to treat many values of k 


simultaneously. 
adequately. 


I. Introduction 


g Let x be a real, non- 
principal character (mod &) and let 


Let k be a positive integer. 
| 


_—- a 


=. x(n 
Lis, x) > x ). 
n=) I 


By a computation using the methods of [1] it was 
shown that if 2<k<227 and k#163, then L(s,x) 
has no positive real zeros. This computation was 
laid out by G. Gourrich. Computation on IBM 
equipment was furnished by Miss L. Cutler and 
E. Rea under the direction of E. C. Yowell. Compu- 
tation on desk computers was furnished by Miss. L. 
Forthal and W. Paine under the direction of G 
Blanch. 

For k= 163, the method of [1] definitely failed. It 
seems likely that by making a careful refinement of 
the estimates of [1] by means of an extensive compu- 
tation, one could handle the case k= 163. However, 
this did not seem a very profitable undertaking, and 
so a further study of L(s,x) for positive real s was 
made, and various alternative methods were devised. 
Some of these seem clearly superior to the method 
of [1]. One such superior method is a generalization 
of the method of Chowla (see [3]). Using this method 
and a table of characters (mod k) prepared by Miss 
L. Cutler and E. Rea, it was a fairly quick matter for 
G. Gourrich to check that if k<227 and k#163, 
then L(s,x) has no positive real zeros. Even by this 
method, the case k=163 remained very difficult. 
Probably the method will handle the case k=163, 
but it seemed clear that even by this method the case 

=163 would require very extensive computations, 
and it seemed worthwhile to devise still further meth- 
This was done, and a method was finally found 
by which the case k= 163 can be handled rather easily, 
with only a minor computation. 

In the meantime, Chowla and Selberg (see [7]) 
have announced still another method for treating 
the case k= 163. 

One may conclude 
lished that if k<22 
real zeros. 


ods. 


hat it is now quite firmly estab- 


t 
7, then L(s, x) has no positive 


By means of this formula, the difficult case of k 


163 was finally treated 


II. Generalization of Chowla’s Method 


Throughout this section, we lay down the follow- 
ing conventions. x is a real primitive character. 
K shall be a positive integer, ¢(n) shall be a function 
of positive integers which is periodic with period A, 
L(s, ¢) shall denote the function that is got by analytic 
continuation from 


so OW) 
hema ; 
n=1 ft 


so that for R(s 


= P(N) 
— 


s 


Lis, >) 


n=1 % 


and f(2,@) shall denote the function that is got by 
analytic continuation from 


‘ . 
S* d(nje~’ 


n=l 


so that for R(z)>0, 


f(z, %=>5 onde 


n=! 


We cite without proof various results, the proofs 
for which can easily be derived from the results of 
chapter XIII of [5]. 


Theorem 1. For s#1, 


h 
L(s,¢)=K ‘D5 o(n)s(s, 


n=l 


a 
K ) 


Corollary. If 


then L(s,.@) is analytic for all s. 


! The preparation of this paper was sponsored (in part) by the Office of Nava! Research 
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Theorem 2. For all z different from 2ami/K (m 


an integer 


A 
yy 
i Dili je 
¢ n=l 


Corollary If 


K 
| 
S* o(n)=0. 


n 


then f(z,@) is analytic for 27/A, and 


; K n 
(0, d)= 3 ‘ d&n){ 1 —=; }- 
© ds o(n)( K) 


n= 


Theorem 3. If R(s)>1, then 


x’ ' f(x, ddr. 





I'(s) Lis, 6) 
Corollary. If 
A 
>> o(n) 0), 
then for R(s) >0, 


I'(s) Ls, o) : fir, dda 


Since 


, 
3 x(n) 
nn 


0), 


the results stated in the three corollaries hold if we 
replace @ by x and K by k. In particular, if we can 
prove that f(z, x)>0 for x>0, it will follow from the 
corollary to theorem 3 that L(s,x) has no positive 
real zeros. This can be readily proved for many 
values of k. The method of proof is as follows. 
With Chowla, we define 


D,\ ri) on) 
n 
‘ 
Ori 3(N) p: @,(m) 


m=! 


We then have the following known result (see [4]) 
Theorem 4. For s>0, 
f(z, @)=(1—e~*) f(z, @,) 

Then, if there is an r for which x,(n) is nonnegative 
for n>1, we infer that f(r, x)>0 for 2>0, and hence 

that Z (s, x) has no positive real zeros. 

For many k’s, we can prove by a brief computa- 
tion that there is an r such that x,(n)>0 for n>1. 
For a typical case, consider k=53 and let x be the 
real nonprincipal character (mod 53). In this case 


x(n), xi(m), and x(n) are sometimes negative, bu 


x;,(n) is nonnegative for all positive n. The com 
putations on which this statement is based ar 
given in table 1 (at the end of the paper). Th 
method of computation of table 1 is particular! 
simple, since from the definition of x,,, we have 
Xroi(h+ 1)— x-41(2) +x-(n+ 1) 
So for r>1 and n>2, x,(n) is the sum of the num 
bers immediately above it and immediately to th 
left of it in table 1. 
We have 
3)=x,(53)+ >> 


— 
m= 54 


~" 
~ 


xin x(m) 


However x,(53)=0 (see table 1) and x is periodi 
with period 53. So 


n 


, 
x(n + 53) S x(m) 


m= 1 


xi (4). 


SO x;(7) Is periodic with period 53. Similarly, xo" 
is periodic with period 53. However x,(53)=742 


(see table 1) so that 


x;(n + 53)=—=742- 


Xa(7) 
As x;(n)>0 for 1<n<53 (see table 1), it follows 
that x;(n)>0 for n> 1. 

In table 2, we have listed those k’s<227 for which 
there is a primitive x and for which we could find 
an r such that x,(n)>0 for n>1. Opposite each 
is listed the least value of r for which x,(n)>0 for 
n>1. Opposite 8 in table 2 is given the r corre- 
sponding to the character x(1)=1, x(3) l, 
x (5) 1, x(7)=1, and opposite 8* is given the 
corresponding to the character x(1)=1, x(3)=1, 
x (5) 1, x(7) 1. The corresponding x’s (mod 
24) are indicated by entries 24 and 24* in table 2 
Similarly for 40 and 40*, 56 and 56%, ete. 

It will be noted that so far we are using exactly 
Chowla’s method (see [3]), though our justification 
for the method is different from Chowla’s. As noted 
by Heilbronn in [4], there exist values of & such that 
no x,(n) is nonnegative for every n>1. In fact one 
can prove that k=163 is such a k; for actual com 
putation for k=163 discloses that f(log (10/7), x) is 
negative, so that by theorem 4 there cannot be any 
completely nonnegative x,. 

Our efforts to find an r for the cases k=43, 67, 88 
123, 148, 173, 187, 188, and 197 were sufficiently 
unrewarded that we suspect that for these values 
of k also there is nor. At any rate we devised an 
improvement of Chowla’s method to handle thes 
intractable k’s (except perhaps k= 163). 

Theorem 5. If 


F(s) 


r-'f(x)dz, 





then 


(1+ ar~*")F(s) ro") f(2)+af(rr)}dz 
i) 
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Theorem 6. If r is a positive integer and 6(n) is the 
coefficient of n~* in the Dirichlet series expansion of 
(1 +ar ‘) Lis, ®), 
then f 
I(x, 8) Mx, d&)+alf(rz, >) 
illustrate the use of these theorems for 
Note that for k=67, x(2)=x(3) 


We now 
the case k=67 
] So 


IC] x(P)p ie 


Te) 


———~ II (1 —x(p)p~*) 
Ip 7 
So the Dirichlet series expansion of 

(1+2-)L(s, x) 

consists of the Dirichlet series expansion of L(s,x) 

with all terms n~* removed for which n is divisible by 

2. Similarly, to get the Dirichlet series expansion of 

(1+: \1+3 *\L(s, x), 


for which n is divisible by 
Similarly for 


we remove all terms 7 
either 2 or 3 
*\(14+-37 "(1 5-*)L(s, x) 

We now verify by actual computation that 
y(n)>0 forn>25. Define d(n) to be the coefficient 
of n~* in the Dirichlet series expansion of 

(1+2-")L(s, x) 
Then 


\T'(s) Lis, x) ‘f(r, o)dz 


ir, 2) ©, X2) 7 T ry? f(2z, x2) 


\s each side of this is a power series In e~’, corre- 
ponding coefficients must be equal. So 


(2n) = xo(2n) + x(n) + x2(n 


o(2n+ 1)= x2(2n+ 1)4+ 2x2(n). 


Reealling that x(n) > 0 for n> 25, we see that surely 
— 
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¢(n)>0 for n>51. Actual computation of ¢,(n 
for 1<n<50 discloses that actually @(n)>0 for 
n>31. 

If we now define @(n) as the coefficient of n-* in 
the Dirichlet series expansion of 


3 *) Lis. ae 
then 


x'f(xr, Adz 


and 


f(r, @) f(x, o)+ f(32, o) 

This latter equation gives a relation between @, and 
¢@. Using this relation with the result ¢(n)>0 for 
n> 31 leads to the inference that 6,(n)>0 for n>95 
Actual computation of 6@,(n) for 1<n<94 discloses 
that actually 6.(n)>0 for n> 41. 

We proceed one step further, defining (n) as 
the coefficient of n~* in the Dirichlet series expansion 
of 

‘y1+5 *) Ls. x) 
Then 


*\P(s) Lis, x) f(s n)dax 


0 


Also we find that 9.(n)>0 for n>29. We now ascer- 

tain by actual computation that n,(n)>0 for n<28, 

and so conclude that ,(n)>0 for n>1. Then 

f(x,n)>0 for r>0 So 
(1+2-3 +3 *)\(1 +5 "I (s) Lis, x) >) 

0 for s>0. 

For other k’s we proceeded similarly, except that 
when x(—1)=1 it was necessary to work with xz, 
$3, 93, etc., instead of with xo, d, 4, etc.,/because the 
latter are periodic when x(—1)=1. 

In table 3 we have listed against * the combination 


for s>0 and so L(s,x) 


It l . P, Ls, x). 
which was used with that * and also the least value 
of r which sufficed with this combination. 

We note in passing that we can handle the case 
k—=43 with the combination 


(1+2-')L(s,x), 


but we need to take r=9, so that it is less laborious 
to take the combination 
(14--2-*)\ 3 *) Lis,x) 
for which r=2 suffices. 
We suspect that this method will work for * 
However for k=163, the combination 


163. 


(1 + po") L(s,x) 
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will require at least five factors, and likely more 


than five, and the 
to be very extensive 


computations involved 
So we sought other methods. 


Ill. Another Method 


We let x be a real primitive character. 


Define 
&(8) xr’ T (5) es) (1) 
&(8,x) (;) od (5) Ls. (2) 
if x(—1)=1 
E(s,x) (z) r(*3*) Lex > 
if x(—1 l 
We have the known results 
t(x) e(] (4) 
E(S, x) E(1 8.x) (2) 


Also known is 
Theorem 7. If k, and &, are any real constants, 
then for ki <a0<ky, 


and &(s, x) are bounded 
Theorem 8. There is a function f(z) with the fol- 
lowing properties: 
(1) f(z) is analytic for R(x) >0. 
(11) f(x) is positive for r >0 
(ill) f(x) is monotone decreasing for zr >v. 
(IV) For z>1, 


») 


“ED A J exp( ED 


‘V) There is a positive constant A such that for 


f(x)< exp/( 


r >v, 
| 2rzr\_ 
exp ( )s Ae: 
V tr l \ k 
(V1) For «a>e, 
E(s)E(s, x = fix\dr 
Proof. Define c, as the coefficient of n~* in the 
Dirichlet series expansion of ¢(s)L(s, x), so that 
foro >|, 





appear 





> ¢ 
. On 
¢(s) L(s, x) > > =. 

n=j/l 


Then ¢;=1 and ¢,>0. 


Case 1. x(—1)=—1 We have 


(4p) “Gen = ff ero i) 


J, J, emr*exo( 


Multiplying by e, 


rn(v+w)\dv dw 
k ) 


- v Ww 


. 
and summing gives for «><a, 
E(s)E(s, x) 
P ‘ 4 rni?—-w) dv dw 
(pw)*’* > Ce exp/ —_ )\ 
/0 0 ié l vk ,v Ww 
We now treat the integral on the right as a double 
integral over the first quadrant. Since the integrand 


is symmetric in v and w, we may replace the integral 
by twice the integral over the area in the first 


quadrant below the 45° ray thru the origin. In this 
we introduce 
Zr \ pw 
y 7 u 
as new variables of integration, getting 
E(s) ES, x) 
— P \ \ ; wnvy\) dy 
} x’ 'dr >» exp ) = , 
0 eo 2 (n=1 vk ) vy 47° 
We take 
, a ny \) dy 
f(x)=4 pate exp/( : ) | | 
e/ 2 (nel vk ) vy 47° 
and conclude (II) and (VI) of our theorem. We p 
y 22 +2 in eq 6 and get 
‘ > mn(2zr+-2) dz 
f(x) | ) Dy en exp/ ) 
Jo (n=l vk Yy2V4r+e2 
Then (1) and (IIT) of our theorem are evider 


Further, for zs >1, we write 














f(r)=4 exp(- 


<4 exp/ 


< exp( 7) AG exp( 7 
So (LV) holds. 


By eq 7, using only the first term of the series, we 
get 


f(x)> |, exp( a =) V2 ie t 
24 exp ( — AE) am 

. e\d 

> sateven(9fen(-): 


oO (V ) holds. 


In the relation 


K,\ Z exp { 2 cosh ) do 
see [5], p. 384, Example 40) put cosh ¢=1+t. There 
results 
: dt 
K (2) exp (— 2(1+4?)) . 
/0 yt? T 2t 


So for x >0 


» a i - 2 iz t t 
1K,(=""" ) | exp(— mnaz(l ') d 
vk JI yi yt >t 
Putting t= 2/27, we get 
— ’ 27-42) 
1Ko( === ) | exp ( WN LI ) d 
vk /0 vk yk \ tr 


™ by eq ft 


Case 2. x(—1) 1. By the duplication for- 


mula for the gamma function 


2vk(=7) rin 
\ ( ED (s)7 
2yk |, Z exp ( “7 ) lr 


2arzr é wn(2zr +2 
ff Shecemn(—eet®, 
vk /0 (n l vk 


wn(2 -+- 2) 


2rzr | ( 
[> exp/( + 
yk ) 0 S=1 vk 


(. 


) 


Multiplying b: 


E(s) E(s, x) 2k 
‘ 


We now take 


f(x) 


and conclude ( 
We can easily 


tHarj< 


and sO | lV } 


vives 


So we may ta 


Theorem 9. 
E(s)E(s, x 


Proof. 


Theorem l( 
/ 


Proof. 


Use 


==) dz 
vk Yaz y4r- 
=) dz 

yk Vy2 \4 T 


yc, and summing gives for «>a, 


2rnzr\) 
Jide 


‘ | | 
Ir Cy, @X » ( 
0 (t2% | vk 


. ae al a 
2yh 2, Cn exp ( ) 5) 
n=1 vk 


1), (11), (111), and (VI) of our theorem. 
show that for ari, 


») 


TE )) f(1) exp ( =) 


\ 


exp ( 


holds. Using the first term of eq 8 


2x 
f(x) >2 yk exp : 
(23) 
ke A 2yVk in (V). 
For all s 
&(1, x) ’ aa , 
) + (z*-'+-2-*) f(z)dz 
RUN } J1 " 


the method of Heilbronn (see [8]). 


). For R(x) >0, 


l ; 
( ) x f(z)+(2—-1)&(1, x). 
r 


Comparing (VI) of theorem 8 with theorem 


vives for o >O; 


"_ 
al 
0 


We 


can 


'f(x)dz 


rewrite 


eC, x) + | x 'f(x)dz. 
1) 1 ‘ 


SUA 


this as 


z*-'#(z)dz= (1, 


" 
20 | x 'dz-+ 








Replacing z by e~‘ in these integrals and using the 
uniqueness theorem for Laplace transforms, we infer 
that for 0<r<l 


So our theorem holds for 0<r<1. Then by an- 
alytic continuation, it holds for R (x) >0. 


If we let z—- in this result, we have zr f(r)->0 


and 
f(-)  f(1)>0 


So we have an alternative proof of the known result 
that &(1, x)>0. If we differentiate both sides and 
put r= 1, we infer: 

Theorem 11. 


E(1, x) f(1 2 f’(1) 


Another consequence of theorem 10 is that 


i202) 


as r—0. Hence in (VI) of theorem 8, we can take 


o,= 1. 


Theorem 12. For 0<s< 1, 


E(s)E(s, x) | is °+-3™ E(1,x)—azf(z)}dz. 
* 


Proof 


Use theorem 9 


It is clear from this that if z f(r)< &1, x) forl<z, 


then L(s,x)>0 for s>0. There are many k’s for 
which xf(r)<é1,x) for 1<z. For example, let 
x(— 1) land k<39. Then 
2x 
>] 
yi 


Hence by eq 8 


' 2mnzr 
ACs 2k SOc, exp ( ) 
n | yk 
. . Qarnr 
rf(x)=2ykD oer exp ( ) 
n=! yk 


However. with 2r/,k>1, 
V 


z exp (—"7"*) 


is decreasing for x>1 So rf(x) is a decreasing 
function for r>1. So we have only to prove 


(i): &(1, x) 





As 


i TT ON ATT SN Ree 


So, by theorem 11, it suffices to prove is 


f’(1)>f(1). 
However 
») ») 
> ,. 2En 2mn 
I’ VD=2 yk >> exp(- ) H 
n=] yA vk 
9 
~ = ii 
yk> ( exp( ) 
n=) yk 
2 
>t ( 1) 
If x(—1)= 1, one can carry out similar reasoning 


based on eq 7, since 


‘ Whe dz 
[~ xp( ."%2) 
0 k yiy4er 


\ 


is a decreasing function of x 

Thus we conclude by a very simple reasoning 
which does not even involve inspection of the values 
of x, that for k<39, Lis, x) has no positive real 
Zeros. 





For larger k’s, we would need to know the Values \ 
of x(n) for some of the smaller values of n. How- 
ever, usually a knowledge of the values of x(n) for 
n<»yk would be more than ample. 

Unfortunately, this method is not general. In f 
particular, it fails for k=163. Indeed (as we will . 
show in the next section) for k= 163, 

JO)> 1, x) 
so that it is impossible to have rs f(r)<&1, x) for 
r>1 So for k=163, more subtle methods are re- t} 
quired. 
IV. Treatment of k=163 
if 
Throughout this section, let k= 163. 
Then x(—1) 1. Also, the class number of ,—t 
is unity, so that 
&(1, x) \ k. 4 
Temporarily define 
g(r) &1, x) x f(x) (10 
Lemma 1. For l<2< yk/2n, g/x)<0. S, 
Proof. Using the first term of eq 8, we get 
2nz 
g(z)<&(1,x)—2yke exp (— ) 
vk . 
H 
2rz 
vk (1 22 exp (- )): 
vk 


») 


2rzr 
vb) 


f exp ( . 





ues 


, 
eal 


ues 
W- 


for 


In 


vill 


for 


re- 


is increasing for 1 <z: 











yk/2n, it suffices to prove 
9 


s<ew(-"") 


However, for k 163, 
9 


exp ( ==) 0.61, 


yh 


and our lemma is proved. 


Lemma 2. For (,k/2m)<2, g(x) >0 


Proof. For (yk/2x)<z, 
d 2rzr 
] ) exp ( )i <0 
( vk 


and for n>1 


d 2 iL 
ffs exp (—?™™4)<a 


Multiplying by 2¢,.,k and summing gives 


d 
da 


riir)-<O0 


for (yk/2r)<r. This establishes our lemma 


Lemma : There is a number a such that 
1) l<ia 
1) For l<r<a, g(r)<0 
II1) For a<r, g(r) >0. 
Proof. Obvious by lemma 1, lemma 2, and the fact 
that 


Lim g(r)=&1, x) >0 


Indeed, we have a > \ I 2r. but this fact isnot needed 
Lemma 4. For 1/2<s<l, 


sa monotone, nondecreasing function of s 
Proof By theorem 12, 


E(s)E(s, x ; I J Q\r) 
: di 
a a Jl a a I 
™) 
Oj E(S)E(S, x) ( 0 \2 tr "i gla j 
add 
ys f a a” *\ J) Os la’ a~*} 7 
However, 
O\r I ; (aa 1) (log r—log a 
Os la a ) (ar)(a*”-*-+-a 
ae a~ (log I og a 
lar) \a  -@ 


So for 1/2<s, we see that 


d (2*-!+27*) 
Os (a*~'+a7*) 
is negative if 1<z<a and positive if a<r. So by 
lemma 3, we see that for 1/2<s 
re) (az tt r *) g(r) 
Os fa‘ tt a *\ r 


is nonnegative for 1<z. Thus our lemma is proved. 
If we can prove that the function 


tis)e(s ) 
S\S/S\8, X 


a a" 


is positive at s=1/2, we can conclude that it is posi- 
tive for s>1/2._ As it is unchanged by replacing s 
by 1—s, we could conclude that it is positive for 
O0<s<l. So we have proved 

Lemma 5. If L(s, x) is positive for s=1/2, then 
it is positive for s>0. 

We have now reduced the problem to that of 
proving that Z(1/2,x)>0. This result is proved in 
7]. Alternatively, one can compute L(1/2,x) explic- 
itly by means of the formula 


This is not exactly an agreeable computation but 
by using a table of the incomplete gamma function, 
it is not unduly laborious. The computation was 
carried out by Miss S. Marks under the direction of 
G. Blanch, and the value 


L | a, x ) 0.06910 


was obtained 

\ shorter computation results if we use theorem 9. 
Putting s=1/2 in this, we see that it suffices to prove 
that 


$E(1, x 2] ] f(r)da 
is negative. That is, using eq 8 and eq 9, it suffices 


to prove 
2rna ; 
exp ( )de 
yh 


That is, it suffices to prove 


So: fn ot fe [en 
\ yk 


y2yh yf 
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where 


( exp ( *#) dt. 


1 ee & 


erfe (y) 


Now for k= 163, we have x(p 1 for every prime 
p with 2<p<37. So for n<40, ¢,=1 when n is a 
perfect square, and zero otherwise. So we need only 


show 
l l 
4 erfe » ad erfe al 
v2yk (24) . (+7) 
l T 
+— erfe {6 } 
o ( \ \ .) 


that is, we need only show 


+ 


0.1979> erfe (0.9921)+-— erfe (2(0.9921)) 


l 
+— erfe (3(0.9921))+. 
» 


Using a four-place table of erfc, we get 


erfe (0.9921)=0.1606 
I 
5 erfe (2(0.9921))=0.0118 
— erfe (3(0.9921 0.0005 
» 
l 
‘ erfe (4(0.0021 0.0000 


As 0.1979 >0.1729, we conclude that L(1/2, x) >0 


V. Miscellaneous Results 
If Lis,x 


is to be positive for 0<s, then from 


. 


I'(s)L(s, x) x*-'f(x,x)dz 


we see that 


xr 'f(z,x)dz 
0 


must be positive even when f(z,x) is negative for 
some xz (as when k=163). In an effort to prove this 
for various k’s, we undertook to find more about 
the behavior of f(z, x Nothing particularly useful 
for our purpose was discovered However miscel- 
laneous results were found, and we list them without 


proof. 


Theorem 13. For s¥1, 
il S se 
L(s, ¢)= -—.—. (—z)*"'f(z,. edz 
2 ri 





where the contour begins at + © on the real ay 

encircles the origin counterclockwise once, a 

returns to+ © on the real axis, and where the « 

tour does not encircle any of the points +2xmi/K 

(m a positive integer), and where|arg(—2)| <7 
Theorem 14 


K 
>> o(n)=0, 


n=l 
then 


f™ (0, @¢)=(— 1)"L(—m, ¢). 


Define 


me _ , aS a 
T,=T(m,K, ¢)= or) (1-Z) 


r=1 


Theorem 15. For M>1, 


( 1) 


VV 
"M > > \ 
K = 


Tu 1)* gy, (A)A* (1), 


where A* (1) denotes the N-th difference of x™ evalu- 
ated at r=1. That is 


A* (1%) =—( yes 1) N! (1+r)* 
re (N—r)!r! 
Theorem 16. If 
h 
> o(n)=0 
n= | 
then 
' "m 
(mi) K"™(,,, m+1,, ,(m+)na pm 
f'™(0.¢) 4 Re B,7 
(m+1)!1 2 2! 
(m+ 1)(m)(m— 1)(m— 2) wp / 
11 BLT, -3+ ( 





where B,, is the m-th Bernoulli number as defined on 
p. 125 of [5]. 


Theorem 17. For x positive and small 


*S‘A.(1—e-*)’, 


F(z, >) é 


where 


l d d d , 
o AL dz" Ma; rs ): ‘(45+ , ) fix, | 


Theorem 18. If x is a primitive real characte: 


and x(—1)=—1, then: 
9 
(1) For 7 
VE 2, (k 
f(z ) > 1)™ ) L(2+2m., 
. e > (3; x 











) > 9 ° 
. 2ah «nmi k . 
) . 7 \ (7) 1 
- For I: > $= k t( x) 5° x 
: . © fel ? 14 ( kz \ 
r 2ar/ 
Vk kz ®-) x(r l 2ami 
(2.x ‘ (4) For 2z+-— 
w 2@fe-1 ? ,; kz ) } 
2m? 
k M h 
m4 R-1 1(z, x) . p Ht "(5 ) La 2m, x 
ivi \ (! ) ) ro 2 ~~ x")? T m=0 “7 
w Sa 2a . om XI rs | 2m \ ‘a 
1) yi eo ) <1 xr) ] 
. 9 9 oma i+2M 
; 2rmi 2m \2a a ¥, 14 k ) 
> For = > - 
h re gi 
TABLE | Characters and sums of characters (mod 53 
j \ k ke <x") ] 
x 
x 2afo7 7? kz r x x x x 
| + 
2m 
] l l 
2 1 0 " 
» “ ] 
~ a Win 
} For x 4 ‘ ‘ 
h ' ; 
t l ( 1 ( 
7 l ( 
Vf 24 m ‘ l 0 
yr kz? 
/ > > 1) (= ) L 2 2m. X 10 | 4 
T m= ot 2 
11 | 1 
12 ] & . 
. > w+ l ! lt 2 
+ 1 Vu yi h ) Ss x\") l . 14 2 13 ‘ 
ao at 2+ 2M kee \3 ! le 6 
“ 14 ( < ) : ~ 
2ar 17 l : 4 103 
Is ‘ 20 i 
rere ‘ i9 ; 2 164 
Theorem 19. If x is a primitive real character 20 2 ‘ 19 
and x( 1) 1, then: 21 ! % 23 
22 l ( 5 ths 
23 I 4 30 
9 24 l +4 4 
. a 25 l 
1) For |z ; ’ 
k at th “un 
27 4 440 
a l 4 ‘74 
2y $5 an 
vk kz\?" . 0 ~ 
f x) > » ( 1) ( - ) Liit+2m x , , - 
wT m= 2n 4 s 
2 2 {2 aL 
3 ; 20 430 
, ‘ ‘ 25 wi 
9 . Ink 2rmi 20 is4 
(2) For !z'- ss 
k k oe ‘ 16 TOO 
37 ; 13 l 
”s l 2 ll 724 
a) ; 5 732 
f \ I < 1 x l 0 » ms - 
T | j kz = l 42 
+( ») -) 43 l 1 ( 42 
on! 4 0 42 
45 1 4] 
E Jeo \2m R-1 (r) 46 “ 
1. > > | (5-) ra +2m,x)— >> abil ( ‘s 0 ~o 
T m= <7 ( —e 8 uv ! 1 41 
“) i ] 42 
9 ° 42 
> aA = Wht 2 1 ‘ 742 
3) For oft 42 
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TABLE Values of k and corresponding values of 1 TABLE 3 Values of k and corresponding combinations ay ; 
values of r 
I 
r k r k r 
k Combination 
‘ 
109 2 15 2 
‘ l * 1! l 147 ! 43 i+2-° 1+3-*) Lis,x 
2 6° 2 2 14s 2 67 | (1+2-*)(14+3-) (145-*) L(s,x 
l 2 1! i 1an* 2 AS (1+3-*) L(s,x 4 
5 2 ' it 172 2 123 | (14+2-*)(14-5-*) L(s,x : 
as a 148 | (14+3>-*)(14+5-*) (14+7-*) L(y 
s 60 119 17 2 173 | (14+2-*)(14+3-*)(14+5-*) L(a,y ‘ | 
f 2 a 179 2 187 | (14+2-*)(143-*) (145-*) L(s,x 
» » ‘ ) a” 2 Is! 2 IRS 1+3-*) Lis,x 
I 2 is l a4 2 Iso 1 197 1+2-*) Lis. 2 
! “iy 2 2 Is4 2 
17 2 71 120 2 In4* 2 | 
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Forced Oscillations in Nonlinear Systems’ 
By Mary L. Cartwright 


his paper shows how the approximate form of the solutions of a certain nonlinear 
differential equation occurring in radio work may be obtained from certain general results 
and gives the proof of the general results in detail. The proof of the general statement 
depends on a type of method that can be applied with minor modifications to any equation 
of the type 
F+k f(r)rt+e(z kp(t 


where p(t has period 2#/\, and p(t)dt is bounded for all ¢, f(r)>1 for 2x >a, and 
0 
g(r)/z> lfor xr >a 
For some years Professor J. E. Litthkewood and I suggested that it corresponded to a physical systen 
have been working on nonlinear differential equa- investigated by him and van der Mark. Fo 


tions * of a type which occur in radio work and _ certain values of the parameters the physical systen 
elsewhere. One of the most interesting of these had two possible stable oscillations, one of period 


equations is 4ne/d and one of period (2n+1)2x/yk. As a matte! 
of fact in the case of (1), owing to the strictly 
g=k(1 — 2*)¢+-2—bkd cos (At+ a), (1) symmetrical nonlinear function 1—z?, the perio: 


tna/d does not occur, but for certain values of 6 the: 


spec > k large 0<b< 2/3 , nti . . . 
especially for & large and »< 2/3. Our attention 4-6 two stable oscillations of periods (2n 4-1)2x/X. 


was drawn to it by a remark of van der Pol,*? which 
a . It would take too long to give a complete pro 


This paper contains material presented in lecture form to the staff of the of this statement here, but I propose to show ho 


Institute for Numerical Analysis of the National Bureau of Standards on January . . . 

28, 1949. Miss Cartwright was a consultant at the INA at the time this lecture the approximate form of the solutions may b 
s delivered > . . . , oUy » Oc 

Tae Te TL. Cte ene 2. & tntiee 2.1 obtained from certain very general results, and giv 


ondon Math. Soe. 20, 180 . : ae 
189 (1945), and Ann. Math. 48, 472-494 (1947): also M. L. Cartwright. J. Inst the proof of the general results in detail. The pro 
Elec. Eng. (Radio Section) 95 (III &-# (1048, and Proc. Cambridge Phil. Soc. © 
45, 495 (1949 
B. van der Pol, Proc. Inst. Radio Eng. 22, 1051-1086 (1934 ‘ B. van der Pol and J. van der Mark, Nature 120, 353-364 (1927) 
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of the general statement depends on the type of 
method which we use throughout. It can be 
applied with minor modifications to any equation 
of the form 


+k f(r)e- g(r) k p(t), (2) 


. 
t 


where p(t) has period 27/\, and p(tdt is bounded 


JV 
for all t, f(x7)>1 for x)/>a, and g(xr)/x>1 for |x\>a, 
provided of course that f, g, p satisfy the usual con- 
ditions required for the existence and uniqueness of 
solutions, and are either independent of & or satisfy 
certain simple inequalities independent of k. 

We assume throughout that #>1, and usually 
that it is large. We first observe that the @ in (1) 
is merely inserted for convenience. It enables us to 
choose the origin of ¢ as we wish. We shall write 
r,.¢, for the values of 2 and # at t=0. 

There are two main weapons of attack besides the 
equation (1) itself. The first is the integrated equation 


ri—Z,-+ k( Z -r- y ; ro)+ |. rdt 
bk |sin (At+ @)—sin al; (3) 


the second is the energy equation which is obtained 
by multiplying (1) by 22 and integrating, 


t 
— 724-2k (x?— 1)#2*dt + x?— z? 
2bkr Z cos (At+ ae)di. (4) 


In some ways the integrated equation is more 
fundamental than (1), and that is why the X is 
inserted on the righthand side of (1). The energy 
equation may be rewritten with z=y, #,—y, in the 
form 


r* + y*—(xrot+ Yo 


2h) brxF cos (At+ a) —(2? 1)? dt. (4°) 


The righthand side of (4’) will be dominated by the 
term 2°z* when s is large unless ¢ is then very small, 
and it seems improbable that # can be small for 
most of the time that s is large. This suggests that 
r+y* decreases rapidly over any are for which z is 
large, and therefore that + is bounded for sufficiently 
large t. The general result which | propose to prove 
later is the following 


Theorem 1 If x r(t.z,.%,) 18 any solution of (1) 
for which r—2 2 i, when t=O. then 


J Bt < Bk. (5 


rhe re B ik a constant inde pe nde nt o} I and : jor I > # 


t>t,(x,. #,). 


Assuming for the moment the result of theorem 1, 
we may argue as follows: eq 3 can be rearranged in 
the form 


>) vo *7 
= r=b6 sin (At+a)+C a rdt, (3°) 
5) k k. 0 
where (C is a constant depending on z,,%,,6 and a, 
and in virtue of (5) the last term is O(1/k) for large 
values of k and 0<t<2z/dX. Hence there are two 
extreme possibilities: either # is comparable with & 
in magnitude, or x is given approximately by the 
equation z= Y, where 
FLX) ; XN 6 sin (At+ a)+¢ (6) 


) 
» 


Both these possibilities may occur for arbitrarily 
large values of ¢t, and also of course there are transi- 
tions from one to the other. For fixed ¢ and 6 and 
C' such that (—2/3)<6 sin (f+ a) +C< 2/3, eq 6 has 
three TOOots, 


X\< 1 X.<.1< X,;=X;¢,C), 
as may be seen from figure 1. If b>2/3 
b sin (At+a)+C runs outside the interval |[—2/3, 2/3 


for some ¢ in each period; we suppose that 0<b< 2/: 
because this is the most interesting range. It may 
be observed from figure 1 that if Y, and .Y; increase 
with ¢, .Y. decreases and vice-versa, so the solutions 
over O<t<2r/A are approximately as shown in 
figure 2, provided that # <B. 

If we integrate over a second period 24, \<t<4-7 
A, we have (3’) with C” in place of C where 


C"=¢ zat=C 
k Jon) 


Nt. )dt-4 O( 13 ); 


Hence if « <B, the solution is given approximately 
by X,(¢,C") over the next riod. The difference 
(”’—C is actually of the order as ¢/k, which 
occurs in the error term i «3"), but it can be shown 
that in the circumstances considered ¢ is of the form 
7,(1,C)+ Ok Putting this and r=.X,(t,C’) +O(k 

in (3°), and remembering that .Y, and 7, are periodic, 
we see than the nonperiodic error is 











1( l 
VtOvdt O( ) 
kJ | 
F x? 
Feox 
xX, Xe 
-2 : [oemon) 
-! xX~ 
= 
3 
Figure 1, 


515 





























Figure 2 


and so the second wave will certainly be lower than 
the first in the case of Y, by an amount of order 1/f, 
and the solution will perform descending waves as 
shown in figure 3. The second wave will be higher 
than the first in the case of .Y,, and higher or lower 
for .Y,, unless 


(so that C’=(C' approximately). This can only 
occur for C=0, and then we have a solution with 
period 27/X approximately. There is in fact a 
strictly periodic solution with period 22/\ which is 
for reasons which we shall give presently, unstable. 

Returning to the case in which |! is comparable 
with & in magnitude, we suppose first that <> dk, 
where 6>0 for a time of length D/k. Then x changes 
by at least Dé, and this is excluded by (5) if Di >2B. 
Similarly ¢ is not less than —4éé for an interval of 
time greater than )/k, for if it is 2 runs outside the 
strip 2\< B. 

As regards the transitions, if 1—2z? >6>0 and |z 
is large, the second term in (1) is much more im- 
portant than z or 6b ky\ cos(M+a), and so # has 
the same sign as # in s <1—6 and the opposite if 
r >1+4, (provided that # is large and & is large). 
Hence any small but significant deviation from X, 
or .X; is at once corrected, whereas a small but signi- 
cant deviation from .Y, causes # to increase rapidly. 
In the latter case # cannot decrease again until z ap- 
proaches either the .Y, or \, curve which corresponds. 

Finally near 2 1 the term 64 cos (At+a@) dom- 
inates except when M4 1/27. The dips 
of YY, towards r=1 and crests of Y, toward z | 
do in fact occur in these intervals, and a complete 
theory of the behavior of solutions near these points 
is very complicated. However, in spite of that we 


ais near tirn 


can now form a fairly good general picture of | 
behavior of solutions. They settle above s=1 

a long descent in waves of the form r=.X, (t,C) a 
proximately, dipping a little lower each time un 
they reach the neighborhood of z=1. There th 
have three alternatives: (1) another wave zs 

just above z=1, (2) an unstable wave r= NX, ji 
below z=1, or (3) a rush down to the correspondi; 
X, which, as may be seen from figure 1, I 
r 2. The stable oscillations naturally follo 
the first or third alternative, and have period 
(2n+-1)2x/X or (2n—1)2x/d according as they hay: 
n+1/2 or n—1/2 waves above z=1. Separating 
them there are many types of unstable motion fo 
lowing an .Y, curve which begins and ends on z= 1 fo 
part or the whole of the way (or a similar .Y, cury 
after a long ascent near zs 1). From an, curv: 
they may pull up sharply to XY, or rush down to .\ 
at any stage. All types starting near r=1 finally 
plunge down to the .Y, curve, and perform a cor 
responding long ascent to r 1, the whole pheno- 
menon for r<0 corresponding to that for z>0. All 
this can be rigorously established, and moreover the 
relative positions of solutions approximating to 
NX3(¢,C,) and .Y;(¢,C,) remain the same throughout a 
long descent unless (, and (, differ by something 


Is Iie 


which is extremely small for large k, such as « 

We now return to the proof of theorem 1. Tix 
preceding analysis depended to a large extent on th 
fact that 1—2* changes sign twice; in the work which 
follows the significant point is that the coefficient 
of z is positive for large x and that the function g(a 
in (2) has the sign of x. Physicists may conside1 
it intuitive that a system with a restoring force and 
positive damping for large x should have bounded 
solutions, and I hope to show that this is also intui- 
tive by mathematical commonsense. Incidentally 
the lemmas give a good deal of other information 
about the solutions. I shall as usual refer to f as 
the time, and say that ¢, is before or after ¢, according 
as t,<t, or t; >t. The constants B are independent 
of & and ¢, and are not necessarily the same in each 
place unless a suffix is attached, and, as we said 
earlier k>1. 

We first ensure that all solutions eventually com: 
fairly near r=0. 

Lemma 1. A solution of (1) cannot have absolut 
value greater than 3” for all large 1. 

Suppose that this is not so, and that r>3 
Integrating from ¢t, where r=2, 7=2;, 


fort 


we have 


Slliia 


t—?,4 k(~ 


C4 i ; 
r + r,)+ rdt=bk| sia (At+-a) 
» t, 


and so 
rdt 


a) 


3 Ct t)< 


r- kx ( l “)+0 1) 


as t-» ©, where the canstant implied depends on 

The lefthand side tengs to ©, and so r->— ~, b 

this implies that z->— ©, which contradicts 27> 
3° for arbitrarily large ¢, and similai 


Hence r<3" 
z> —3” for some large ¢. 
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We next establish that as long as |x| is not too large 


cannot increase too much. 
on an are PQ, then 


the value of 
Lemma 2. /f x <3” 
ty|<.|¢p| + Bk, 
and, if the are lasts a time longe r than 4 k ~4, rg 
More ge nerally if ix < b. on an are PQ, then 


fp + Bb )k 


If 2 <3" and <> on an are, it can only last a 
time at most 2.3°/k<4/k. For if it lasts longer, the 
solution travels a distance more than 2.3 and there- 
fore cannot remain in “q~ strip. The same is true if 
i k on an are in\r <: 

Let P; be the last ot before Q at which |% <fk, 
or P itself, whichever is the latest. Suppose first 
that zp, >0, then ¢>0 on the are P,Q. For, if not, 
z=0 at some point of P,Q and P, is not the last point 

which |# <k. By the integrated equation, since 


rdt< 


result. 
If the time 


the second 


where B,>1 depends on 6, which gives the 
Similarly if tp,< O, we have Zp, ex Bik. 
than 4k, ip, <k, and we 
form of the result. The 
the same method 

The 


are QR outside the strip jx <3 


have 
result for |z'< 6, 


is greater 
follows by 
next two lemmas show that the height of an 
time taken 
to describe it are bounded by numbers depending on 


, and the 


tg; in other words the velocity with which it emerges 
from the strip 

Lemma 3. /f the are QR lies above r=3" and 
he gins at Qona 3°. the greatest he ight h satisfie r 


I 
Q 4 kh, 


2k 


When r=A=3"%+h,.2=—0. 
point Q@ to the point /7 at 


Integrating from the 


which r=hA, we have 


l)da 


rdt+ bk{sin (At+ a) 


} 
Sin a) 


since ws 


VR, 
O0<ig—2kh, 
and so for h. 
If OR is an are above x 


. then the time t taken 
prov “led that Ig I 


and the result follows for A 
Lemma 4 
at Y on I 3 


less than B Lo, 


. } ‘ . 
»>* OF quining 


to describe QR is 


that #g>k, and that FP is on r=3" so 


rr—3". On QR it is easy to see that 


Suppose 


that r¢ 


*f 


Ie Bk rdt 


Big 


because zg >k 


0 Ir Biot 


and so t< Byte 
The reduction in energy over any sufficiently high 
arc is established in lemma 5 which gives effect to 
the remarks we made about (4’). It is much the 
most difficult part of the proof. 
Lemma 5. Jf QR is an are above x 
at Q on r=3" and ending at R on x==3", for 
B,>1 there exists a By> B,>1 suc *h that if Ig 


;™ begin ning 
a give n 
BA, 
2<#—4B ki 


Rr™*@ 


The energy eq 4 for the are QP is 


bd 


*t 
2k (72?—1)#*dt+-2bkxr f sin (At+ ajdt 


0 


"¢ 
lerdt>k rdt 


2 -+2bkx J, idt 
- Bkt( fe dt) 


by Cauchy's inequality. So by lemma 4 and (8 


to < 2J- 


provided that J > Bikig, where B; depends on B, and 
ig>k. We may obviously choose B}>48,, and we 
have the result except for the which 
J < Bikig 


Suppose 


cause in 


that J< / from 


2/2+3" 


now i keg, and integrate 


Q@ until t=; we or z whichever comes 


. nn ] 
first. Then #>>5%g>0. Hence z is increasing on 


this are and 


; 


sin At+-a rdt sin a dy: a 2 
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Hence 


#42 2k | (2?—1)adt 


> —2J—kBB:-—B 


> — 2k Btig—kB> —? # 


—_ } y 
‘ovided that ¢g>kB 2 bi 

provided that tg >A B,. It follows that z*> jie and 

l ~wiax ek 
therefore #> <7, so that x reaches 2B}+-3” first. 
Now since 2?—1>2 on the are 

°28,+3 
J>2k it dr>kig- 283, 


and we have a contradiction. 
the lemma is true. 

It remains to combine the results of lemmas 2 
and 5. 

Lemma 6. Jf Q is a point on xr=3” such that 
Ig »>Byk>O, then the solution returns to the strip 


Hence the result of 
































r\ <3" at R and emerges again at S with | zs tq 
Byk 
By lemma 2 with 2? and S in place of P and Q 
ig |?<\dg!?+2B klip! + Bl 
2 
| 
t 
x 9 oa 
-| 
-2 
Figure 3 
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kb» sin (At+a)dr—B 





Using lemma 5 in the form stated and also in t] 
weaker form \zg!<\ie!, we have 
¥s\?< \%g\*,—4 Bk ige\+2Bkiie 


+ Bik 
(%o|— B,k)’, 
From which the result follows. For #g>B,k> B,/ 
We now have the result stated intheorem1. Fo 
by lemma 1 the solution must enter the stri 
r<3", and by the second part of lemma 2 z# < B,/ 
if it stays there. If on the other hand it emerges 
at Y@ with tg <B,k, the height (or depth) of th 
subsequent are outside the strip is less than B,+- B, 
For by symmetry all the lemmas for —>3”* have 
strictly corresponding forms for z<—3”". But by 
lemma 6 if |tg >B,k, the solution emerges the next 
time at S with (%s)< |#g' — Bk, so that it must either 
stay in the strip or emerge eventually with # < BA 
and then by the second part of lemma 2, 2 <B,+ 8 
for all subsequent ¢ 
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